Import Modes from FEST3D File

With this option,impressed forward travelling modes in all waveguides of a multi-port network can be imported from a .fim file.



Figure 1. The AW - Waveguide port dialog for importing from a FEST3D file.

Parameters:

Reference point for FEST3D model (named point)
A named point indicating the translation of the imported model in the Feko coordinate system. Note that this point must have been previously defined with a DP card.
Rotation around the X axis
This specifies the rotation of the imported model around the X axis in degrees.
Rotation around the Y axis
This specifies the rotation of the imported model around the Y axis in degrees.
Rotation around the Z axis
This specifies the rotation of the imported model around the Z axis in degrees.
File name
The name of the .fim file.

In order to model a waveguide port excitation by an impressed mode, the cross section of the waveguide at the port location must be meshed into metallic triangles with a unique label. The propagation direction is given by the unit vector u 3 , see the small graphics in the AW card panel above.

In general, specific meshing rules exist in Feko relating the triangular patch size to the wavelength. When meshing the cross section of a waveguide to define a waveguide port, the mesh size must be small enough to capture the field distribution of the highest mode ( m , n ) which is included in the expansion. Feko checks this automatically and gives a warning for coarse meshes or an error if the mesh size is too large. One must then either refine the mesh just at the port or reduce the maximum mode indices used in the expansion.

The following restrictions apply when using a waveguide port excitation:
  • Waveguide ports are available for models containing metallic objects (wires and surfaces and wire/surface junctions, including PO) and dielectrics (solved using the SEP, FEM) or dielectric coatings and thin dielectric sheets.
  • Special Green’s functions may not be used in conjunction with waveguide port feeds
  • When using waveguide ports, then UTD is not allowed in the same model. Note, however, that in Feko it is possible (using the AR or AP cards) to decompose a model (say a horn antenna in front of a reflector) into different sub-problems. If this approach is followed, then waveguide feeds may be used. See Example_35 in the Scripting Examples guide for an illustration of this decomposition technique.

The reflection coefficient at each waveguide port (S11) is always calculated and available for display in POSTFEKO on an S-parameter graph, even when no S-parameter calculation has been requested. Requesting S-parameters with the SP card is supported for waveguide ports. Multiple ports (active and/or passive) can be present in the model. S-parameters are directly based on the waveguide impedance of the specific mode under consideration. The reference impedance as specified at the SP card is not used for waveguide ports.

Examples for the application of waveguide feeds can be found in the Examples Guide (E-2) and in the Script Examples, example_08 and example_34.

In order to rule out any possible doubts and ambiguities regarding the waveguide mode definitions, we give here the explicit expressions of the modes as used in Feko. This implementation follows closely the conventions in S. Ramo, J. R. Whinnery, and T. van Duzer, Fields and Waves in Communication Electronics, John Wiley & Sons, Inc., 3rd ed., 1994.

Rectangular waveguide expressions

Local Cartesian coordinates ( u 1 , u 2 , u 3 ) are assumed. The factor e jωt±jβ u 3 is omitted for brevity, where β is the complex modal propagation coefficient. In the expressions below, let a be the dimension of the waveguide port in u 1 (distance between S1 and S2), and b the dimension of the waveguide port in u 2 (distance between S1 and S3). The modal cutoff coefficient is the same for TM- and TE-modes, and is given by,

(1) β c = ( ( mπ a ) 2 + ( nπ b ) 2 ) 1/2
For transverse magnetic (TM) modes the axial magnetic field component vanishes, and the axial electric field component for the TM m,n mode is expressed as,
(2) E u 3 TM ( u 1 , u 2 )=Asin( mπ u 1 a )sin( nπ u 2 b )
for m = 1,2... , and n = 1,2... A is the complex amplitude of the impressed mode in V m . The remaining field components are,
(3) E u 1 TM ( u 1 , u 2 )= Z TM H u 2 TM ( u 1 , u 2 )= jβ β c 2 ( mπ a )Acos( mπ u 1 a )sin( nπ u 2 b )
(4) E u 2 TM ( u 1 , u 2 )= Z TM H u 1 TM ( u 1 , u 2 )= jβ β c 2 ( nπ b )Asin( mπ u 1 a )cos( nπ u 2 b )
(5) Z TM = β ωε

For transverse electric (TE) modes the axial electric field component vanishes, and the axial magnetic field component for the T E m,n mode is expressed as,

(6) H u 3 TE ( u 1 , u 2 )=Acos( mπ u 1 a )cos( nπ u 2 b )
for m=0,1,2,..., and n=0,1,2,..., (but not both m and n zero). A is the complex amplitude of the impressed mode A m . The remaining field components are,
(7) E u 1 TE ( u 1 , u 2 )= Z TE H u 2 TE ( u 1 , u 2 )= jβ β c 2 ( nπ b ) Z TE Acos( mπ u 1 a )sin( nπ u 2 b )
(8) E u 2 TE ( u 1 , u 2 )= Z TE H u 1 TE ( u 1 , u 2 )= jβ β c 2 ( mπ a ) Z TE Asin( mπ u 1 a )cos( nπ u 2 b )
(9) Z TE = ωμ β

Circular waveguide expression

Local cylindrical coordinates (r,φ,z) are assumed with the Z axis on the waveguide axis (S1-S3). The factor e jωt±jβz is omitted for brevity, where β is the complex modal propagation coefficient. The expressions below are valid for the fields inside a circular waveguide, for example, ra , where a is the radius of the waveguide port. J m (x) is the m th order Bessel function of the first kind, and J ' m (x) denotes the derivative with respect to the argument.

For transverse magnetic (TM) modes the axial magnetic field component vanishes, and the axial electric field component for the TM m,n mode is expressed as,

(10) E z TM (r,φ)=A J m ( β c r)cos[ m(φ φ 0 ) ]

for m=0,1,2,..., and n = 1,2... A is the complex amplitude of the impressed mode in V m and φ 0 is the rotation angle. The modal cutoff coefficient β c is the n th zero of J m ( β c a)=0 . The remaining field components are,

(11) E r TM (r,φ)= Z TM H φ TM (r,φ)= jβ β c A J m ( β c r)cos[ m(φ φ 0 ) ]
(12) E φ TM (r,φ)= Z TM H r TM (r,φ)= jβm β c 2 r A J m ( β c r)sin[ m(φ φ 0 ) ]
(13) Z TM = β ωε

For transverse electric (TE) modes the axial electric field omponent vanishes, and the axial magnetic field component for the T E m,n mode is expressed as,

(14) H z TE (r,φ)=A J m ( β c r)cos[ m(φ φ 0 ) ]

for m=0,1,2,..., and n = 1,2... A is the complex amplitude of the impressed mode in A m and φ 0 is the rotation angle. The modal cutoff coefficient β c is the n th zero of J ' m ( β c a)=0 . The remaining field components are,

(15) E r TE (r,φ)= Z TE H φ TE (r,φ)= jβm β c 2 r Z TE A J m ( β c r)sin[ m(φ φ 0 ) ]
(16) E φ TE (r,φ)= Z TE H r TE (r,φ)= jβ β c Z TE A J m ( β c r)cos[ m(φ φ 0 ) ]
(17) Z TE = ωμ β

Coaxial waveguide expressions

Local cylindrical coordinates (r,φ,z) are assumed with the Z axis on the waveguide axis (S1-S3). The factor e jωt±jβz is omitted for brevity, where β is the complex modal propagation coefficient. The expressions below are valid for the fields inside a coaxial waveguide, for example, for r i r r o , where r o is the radius of the outer conductor and r i is the radius of the inner conductor of the coaxial waveguide port. J m (x) and N m (x) are m th order Bessel function of the first and second kind, respectively and J ' m (x) and N ' m (x) denote the derivatives with respect to the argument.

The fundamental mode in a coaxial waveguide is a TEM wave and propagates with β=ω με .The axial electric and magnetic field components are zero for a TEM-mode, and the transverse field components have a static field distribution,

(18) E r TEM (r,φ)= Z TEM H φ TEM (r,φ)= A r
(19) Z TEM = μ ε

A is the complex amplitude of the impressed mode in V

For transverse magnetic (TM) modes the axial magnetic field component vanishes, and the axial electric field component for the TM m,n mode is expressed as,

(20) E z TM (r,φ)=A( J m ( β c r) J m ( β c r o ) + N m ( β c r) N m ( β c r o ) )cos[ m(φ φ 0 ) ]

for m=0,1,2,..., and n = 1,2... A is the complex amplitude of the impressed mode in V m and φ 0 is the rotation angle. The modal cutoff coefficient β c is the n th root of the transcendental characteristic equation, N m ( β c r i ) J m ( β c r o )= N m ( β c r o ) J m ( β c r i ) , enforcing the boundary condition that E z must be zero at r i and r o . The remaining field components are,

(21) E r TM (r,φ) = Z TM H φ TM (r,φ) = jβ β c A( J m ( β c r) J m ( β c r o ) + N m ( β c r) N m ( β c r o ) )cos[ m(φ φ 0 ) ]
(22) E φ TM (r,φ) = Z TM H r TM (r,φ) = jβm β c 2 r A( J m ( β c r) J m ( β c r o ) + N m ( β c r) N m ( β c r o ) )sin[ m(φ φ 0 ) ]
(23) Z TM = β ωε

For transverse electric (TE) modes the axial electric field component vanishes, and the axial magnetic field component for the T E m,n mode is expressed as,

(24) H z TE (r,φ)=A( J m ( β c r) J m ( β c r o ) + N m ( β c r) N m ( β c r o ) )cos[ m(φ φ 0 ) ]

for m=0,1,2,..., and n = 1,2... A is the complex amplitude of the impressed mode in A m and φ 0 is the rotation angle. The modal cutoff coefficient β c is the n th root of the transcendental characteristic equation, N m ( β c r i ) J m ( β c r o )= N m ( β c r o ) J m ( β c r i ) , enforcing the boundary condition that the derivative of H z normal to the conductors must be zero at the inner and outer radii. The remaining field components are,

(25) E r TE (r,φ) = Z TE H φ TE (r,φ) = jβm β c 2 r Z TE A( J m ( β c r) J m ( β c r o ) + N m ( β c r) N m ( β c r o ) )sin[ m(φ φ 0 ) ]
(26) E φ TE (r,φ) = Z TE H r TE (r,φ) = jβ β c Z TE A( J m ( β c r) J m ( β c r o ) + N m ( β c r) N m ( β c r o ) )cos[ m(φ φ 0 ) ]
(27) Z TE = ωμ β