COMSOL 5.0 RF Tutorials
COMSOL 5.0 RF Tutorials
October 2014
COMSOL 5.0
Contact Information
Visit the Contact COMSOL page at www.comsol.com/contact to submit general
inquiries, contact Technical Support, or search for an address and phone number. You can
also visit the Worldwide Sales Offices page at www.comsol.com/contact/offices for
address and contact information.
If you need to contact Support, an online request form is located at the COMSOL Access
page at www.comsol.com/support/case.
Other useful links include:
Support Center: www.comsol.com/support
Product Download: www.comsol.com/product-download
Product Updates: www.comsol.com/support/updates
Discussion Forum: www.comsol.com/community
Events: www.comsol.com/events
COMSOL Video Gallery: www.comsol.com/video
Support Knowledge Base: www.comsol.com/support/knowledgebase
Part number: CM021002
Branch-Line Coupler
Introduction
A branch line coupler, also known as a quadrature (90) hybrid, is a four-port network
device with one input port, two output ports, with a 90 phase difference between
them, and one isolated port. Due to its symmetry, any port can be used as the input
port.
Isolated
Z0/2
Z0
Output
Z0
(4)
(3)
/4
Z0
Z0
/4
Z0
Z0
(1)
(2)
Input
Output
Z0/2
Model Definition
The form of the branch line coupler is shown schematically in Figure 1. The layout
design is based upon Ref. 1, and is tuned to operate at 3 GHz. The design is realized
as microstrip lines patterned onto a 0.060 inch dielectric substrate. The microstrip
lines are modeled as perfect electric conductor (PEC) surfaces, and another PEC
surface on the bottom of the dielectric substrate acts as a ground plane. The entire
modeling domain is bounded by PEC boundaries that represent the device packaging.
The four ports are modeled as small rectangular faces that bridge the gap between the
1 |
BRANCH-LINE COUPLER
PEC face that represents the ground plane, and the PEC faces that represent the
microstrip line at each port.
Lumped port 2
PEC enclosure
Lumped port 3
Air
Lumped port 1
60 mil substrate
Lumped port 4
Figure 2: The model of the branch line coupler. Some exterior faces are removed for
visualization.
The model is shown in Figure 2. A small air domain bounded by a PEC surface around
the device is also modeled. The model is meshed using a tetrahedral mesh. A good rule
of thumb is to use approximately five elements per wavelength in each material.
2 |
BRANCH-LINE COUPLER
Figure 3: The frequency response of the branch line coupler shows good input matching
(S11) and isolation (S41) around 3 GHz. The coupled signal at the two output ports
(S21 and S31) is about -3 dB at 3 GHz.
Figure 4: The phases on the two output ports show approximately 90-degree shift at 3 GHz.
3 |
BRANCH-LINE COUPLER
Reference
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
4 |
BRANCH-LINE COUPLER
Expression
Value
Description
thickness
60[mil]
0.0015240 m
Substrate thickness
l_s
40[mm]
0.040000 m
Length, substrate
w_line2
5[mm]
0.0050000 m
Width, line 2
l_line2
13[mm]
0.013000 m
Length, line 2
l_line1
(l_s-l_line2)/
2
0.013500 m
Length, line 1
w_line1
3.2[mm]
0.0032000 m
Width, line 1
w_line3
3[mm]
0.0030000 m
Width, line 3
l_line3
13.6[mm]
0.013600 m
Length, line 3
f_min
1[GHz]
1.0000E9 Hz
Minimum frequency in
sweep
f_max
5[GHz]
5.0000E9 Hz
Maximum frequency in
sweep
lda_min
c_const/f_max
0.059958 m
h_max
0.2*lda_min
0.011992 m
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 2*w_line1+l_line3.
4 In the Height text field, type l_s.
5 Locate the Position section. From the Base list, choose Center.
6 Click the Build Selected button.
5 |
BRANCH-LINE COUPLER
Rectangle 2 (r2)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type w_line2*2+l_line3.
4 In the Height text field, type l_line2.
5 Locate the Position section. From the Base list, choose Center.
6 On the Work plane toolbar, click Build All.
Rectangle 3 (r3)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type l_line3.
4 In the Height text field, type l_line2.
5 Locate the Position section. In the xw text field, type -l_line3/2.
6 In the yw text field, type l_line2/2+w_line3.
7 Click the Build Selected button.
8 Click the Zoom Extents button on the Graphics toolbar.
Array 1 (arr1)
1 On the Work plane toolbar, click Transforms and choose Array.
2 Select the object r3 only.
3 In the Settings window for Array, locate the Size section.
4 From the Array type list, choose Linear.
5 In the Size text field, type 3.
6 Locate the Displacement section. In the yw text field, type -l_line2-w_line3.
7 On the Work plane toolbar, click Build All.
Difference 1 (dif1)
1 On the Work plane toolbar, click Booleans and Partitions and choose Difference.
2 Select the objects r2 and r1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the three rectangles belonging to the array object (arr1).
6 Clear the Keep interior boundaries check box.
6 |
BRANCH-LINE COUPLER
Extrude 1 (ext1)
1 On the Geometry toolbar, click Extrude.
2 In the Settings window for Extrude, locate the Distances from Plane section.
3 In the table, enter the following settings:
Distances (mm)
thickness
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
7 |
BRANCH-LINE COUPLER
Union 1 (uni1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Union.
2 In the Settings window for Union, locate the Union section.
3 Clear the Keep interior boundaries check box.
4 Select the objects blk1 and ext1 only.
5 Locate the Selections of Resulting Entities section. Select the Create selections check
box.
6 Click the Build All Objects button.
7 Right-click Component 1 (comp1)>Geometry 1>Union 1 (uni1) and choose Rename.
8 In the Rename Union dialog box, type Substrate in the New label text field.
9 Click OK.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type l_s.
4 In the Depth text field, type l_s+l_s/8.
5 In the Height text field, type thickness*5.
6 Locate the Position section. From the Base list, choose Center.
7 In the z text field, type thickness*5/2.
8 Click the Build All Objects button.
9 Right-click Component 1 (comp1)>Geometry 1>Block 2 (blk2) and choose Rename.
10 In the Rename Block dialog box, type Package in the New label text field.
8 |
BRANCH-LINE COUPLER
11 Click OK.
DEFINITIONS
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, locate the Input Entities section.
3 From the Geometric entity level list, choose Boundary.
4 Select Boundary 13 only.
5 Right-click Component 1 (comp1)>Definitions>Explicit 1 and choose Rename.
6 In the Rename Explicit dialog box, type Microstrip Line in the New label text field.
7 Click OK.
View 1
Hide three boundaries to get a better view of the interior parts when reviewing the
mesh.
1 On the 3D view toolbar, click Hide Geometric Entities.
9 |
BRANCH-LINE COUPLER
2 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
3 From the Geometric entity level list, choose Boundary.
4 Select Boundaries 1, 2, and 4 only.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
Now set up the physics. The default boundary condition is perfect electric conductor,
which is applied to all exterior boundaries. Apply this condition also to the interior
boundaries of the microstrip lines.
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic
Waves, Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
section.
3 From the Selection list, choose Microstrip Line.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 24 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Wave excitation at this port list, choose On.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 25 only.
Lumped Port 3
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
10 |
BRANCH-LINE COUPLER
Lumped Port 4
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 14 only.
Lumped ports are assigned at each end of the microstrip lines. Wave excitation is on only
at the first port.
MATERIALS
Assign material properties to the model. First, apply air to all domains.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
11 |
BRANCH-LINE COUPLER
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Substrate.
4 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
3.38
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
12 |
BRANCH-LINE COUPLER
STUDY 1
13 |
BRANCH-LINE COUPLER
The input power is evenly split between the two output ports.
S-Parameter (emw)
1 In the Model Builder window, under Results click S-Parameter (emw).
2 In the Settings window for 1D Plot Group, click to expand the Title section.
3 From the Title type list, choose None.
4 Click to expand the Legend section. From the Position list, choose Lower right.
1D Plot Group 3
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Title section.
3 From the Title type list, choose None.
4 On the 1D plot group toolbar, click Global.
5 In the Settings window for Global, locate the y-Axis Data section.
14 |
BRANCH-LINE COUPLER
Unit
Description
arg(emw.S21)
deg
S-parameter, 21 component,
phase
arg(emw.S31)
deg
S-parameter, 31 component,
phase
Derived Values
1 On the Results toolbar, click Global Evaluation.
2 In the Settings window for Global Evaluation, locate the Data section.
3 From the Parameter selection (freq) list, choose From list.
4 In the Parameter values (freq (Hz)) list, select 3.0000E9.
5 Locate the Expression section. In the Expression text field, type
arg(emw.S21)-arg(emw.S31).
15 |
BRANCH-LINE COUPLER
16 |
BRANCH-LINE COUPLER
Model Definition
This example considers three geometries:
a rectangular cavity of dimensions 0.9 in-by-0.9 in-by-0.4 in;
a cylindrical cavity of radius 0.48 in and height 0.4 in; and
a spherical cavity of radius 1.35 cm.
The cavity walls are assumed to be a good conductor, such as copper, with an electric
conductivity of 5.7107 S/m, and relative permeability and permittivity of unity. The
interior of the cavity is assumed to be vacuum, with zero electric conductivity, and unit
permeability and permittivity. The analytic solutions to these three cases are given in
Ref. 1.
The lossy walls of the cavity are represented via the impedance boundary condition.
This boundary condition accounts for the frequency dependent losses on the walls of
a cavity due to the non-zero electric conductivity, which makes the eigenvalue problem
nonlinear. When solving any eigenvalue problem, it is necessary to provide a frequency
around which to search for modes. In addition, when solving a nonlinear eigenvalue
problem, it is also necessary to provide a frequency at which to initially evaluate the
frequency-dependent surface losses. Although the guesses for these frequencies do not
need to be very close, solution time is less the closer they are.
It is usually possible to estimate the resonant frequency of interest, and to use this as
an initial guess. It is also possible to quickly estimate the resonant frequency by
building a second model that uses the perfect electrical conductor (PEC) boundary
condition instead of the impedance boundary condition. A model that uses only PEC
boundaries will result in a linear eigenvalue problem, and is less computationally
intensive to solve. Such a model only requires a rough guess at the frequency of the
1 |
mode, and does not require a frequency at which to evaluate the surface losses.
Therefore, it is often convenient to also solve a version of a model without losses.
Q-FACTOR AND RESONANT FREQUENCY IN CAVITY STRUCTURES
(1)
There are two dominant modes for a cylindrical cavity. One dominant mode of the
cylindrical cavity is TE111 when the ratio between the height and radius is more than
2.03. The other dominant mode is TM010 when the ratio is less than 2.03. For this
case, the Q-factor and resonant frequency are given as
1
2,4049 2
1,2025
Q TM 010 = --------------------------, f TM 010 = ------------------ ------------------
a
a
2
R s 1 + ---
(2)
(3)
r
---------- and is the
2
These two analytical approaches are compared with the Q-factor obtained from
Eigenfrequency analysis.
2 |
used by the finite element method better approximate the analytic solutions, which are
described by sinusoidal functions for the rectangular cavity and Bessel functions for the
cylindrical and spherical cavities. This difference between the numerical results and the
analytic solution is discretization error, and is always reduced with mesh refinement.
TABLE 1: RESULTS FOR THE TE101 MODE OF A RECTANGULAR CAVITY
MAXIMUM
MESH SIZE
Q-FACTOR
(ANALYTIC=7770)
h_max
9.706
7039
h_max/2
9.283
7687
h_max/4
9.273
7765
h_max/8
9.273
7770
Q-FACTOR
(ANALYTIC=8065)
h_max
9.458
7891
h_max/2
9.419
8004
h_max/4
9.411
8056
h_max/8
9.411
8065
Q-FACTOR
(ANALYTIC=14594)
h_max
9.752
14121
h_max/2
9.723
14430
h_max/4
9.701
14616
h_max/8
9.697
14641
Note that convergence with respect to the mesh is fastest for the rectangular cavity and
slowest for the spherical cavity. This is because the isoparametric finite-element mesh
represents curved surfaces approximately, via second order polynomials by default.
This introduces some small geometric discretization error that is always reduced with
mesh refinement. Although it is possible to use different element orders, the default
second-order curl element (also known as a vector or Nedelec element) is the best
compromise between accuracy and memory requirements. Because memory
requirements for three-dimensional models increase exponentially with increasing
element order, and increasing number of elements, there is strong motivation to use
as coarse a mesh as reasonable. Figure 1 shows the fields within the cavities, as well as
the surface currents and surface losses.
3 |
Figure 1: Arrow plots of electric and magnetic fields. Slice plot of electric field.
4 |
5 |
Reference
1. C.A. Balanis, Advanced Engineering Electromagnetics, John Wiley & Sons, 1989.
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
6 |
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type a_r.
4 In the Depth text field, type a_r.
5 In the Height text field, type b_r.
7 |
Now set up the physics. Override the default perfect electric conductor condition on
the exterior boundaries by an impedance condition.
Assign material properties on the model. First, apply vacuum to all domains.
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
8 |
Name
Value
Unit
Property group
Relative permittivity
epsilonr
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 From the Selection list, choose All boundaries.
5 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
sigma_wall
S/m
Basic
Add variables for Q-factor calculation and visualization. For this Q-factor calculation,
add two integration coupling operators: one for volume and the other for surface
integration.
Integration 1 (intop1)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_v in the Operator name text field.
3 Locate the Source Selection section. From the Selection list, choose All domains.
9 |
Integration 2 (intop2)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_s in the Operator name text field.
3 Locate the Source Selection section. From the Geometric entity level list, choose
Boundary.
4 From the Selection list, choose All boundaries.
Variables 1
1 On the Definitions toolbar, click Local Variables.
2 In the Settings window for Variables, locate the Variables section.
3 Click Load from File.
4 Browse to the models Model Library folder and double-click the file
cavity_resonators_model1_variables.txt.
The emw. prefix is for the Electromagnetic Waves, Frequency Domain interface in
the first model. Wav and Qsh are Energy density time average and Surface losses,
respectively. Qfactor included in this text file shows up in orange indicating an
unknown variable. It will be known after solving the model.
MESH 1
The maximum mesh size is one dimension of the cavity scaled inversely by d_f, a
discretization factor defined in Parameters. The discretization factor is also used as a
parametric sweep variable to see the effect of the mesh refinement.
Size
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
type h_max_r/d_f.
5 In the Maximum element growth rate text field, type 2.
6 In the Curvature factor text field, type 1.
7 In the Resolution of narrow regions text field, type 0.1.
8 Click the Build All button.
10 |
STUDY 1
Provide the number of modes and a frequency around which to search for modes.
Step 1: Eigenfrequency
1 In the Model Builder window, expand the Study 1 node, then click Step 1:
Eigenfrequency.
2 In the Settings window for Eigenfrequency, locate the Study Settings section.
3 In the Desired number of eigenfrequencies text field, type 1.
4 In the Search for eigenfrequencies around text field, type 9e9.
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
d_f
1 2 4 8
Parameter unit
11 |
5 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
with that shown in Figure 1, top. The exact numbers that you get may differ slightly.
Add a surface plot of the surface losses and an arrow plot of the surface current
(Figure 2, top).
3D Plot Group 2
1 On the Model toolbar, click Add Plot Group and choose 3D Plot Group.
2 In the Model Builder window, under Results right-click 3D Plot Group 2 and choose
Surface.
3 In the Settings window for Surface, click Replace Expression in the upper-right corner
12 |
ROOT
1 On the Model toolbar, click Add Physics to open the Add Physics window.
2 Go to the Add Physics window.
3 In the Add physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
4 Click Add to Component in the window toolbar.
5 On the Model toolbar, click Add Physics to close the Add Physics window.
ADD STUDY
1 On the Model toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Custom Studies>Empty
Study.
You will copy the settings from the existing study later on.
4 Click Add Study in the window toolbar.
5 On the Model toolbar, click Add Study to close the Add Study window.
GEOMETRY 2
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type a_c.
4 In the Height text field, type height_c.
13 |
Set up the second physics interface. The steps are same as for the first model.
Assign material properties to the second model. Apply vacuum to all domains.
Material 3 (mat3)
1 In the Model Builder window, under Component 2 (comp2) right-click Materials and
14 |
Name
Value
Unit
Property group
Relative permittivity
epsilonr
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Material 4 (mat4)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 From the Selection list, choose All boundaries.
5 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
sigma_wall
S/m
Basic
Add variables and two integration coupling operators. The purpose of these is same as
in the first model.
Integration 3 (intop3)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_v in the Operator name text field.
3 Locate the Source Selection section. From the Selection list, choose All domains.
15 |
Integration 4 (intop4)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_s in the Operator name text field.
3 Locate the Source Selection section. From the Geometric entity level list, choose
Boundary.
4 From the Selection list, choose All boundaries.
Variables 2a
1 On the Definitions toolbar, click Local Variables.
2 In the Settings window for Variables, locate the Variables section.
3 Click Load from File.
4 Browse to the models Model Library folder and double-click the file
cavity_resonators_model2_variables.txt.
The emw2. prefix refers to the Electromagnetic Waves, Frequency Domain interface
for the second model.
MESH 2
Apply the same logic in the mesh set up as you have done in the first model.
Size
1 In the Model Builder window, under Component 2 (comp2) right-click Mesh 2 and
type h_max_c/d_f.
5 In the Maximum element growth rate text field, type 2.
6 In the Curvature factor text field, type 1.
7 In the Resolution of narrow regions text field, type 0.1.
8 Click the Build All button.
STUDY 2
1 Select both Study 1> Step 1: Eigenfrequency 1 and Study 1> Parametric Sweep 1 using
16 |
Step 1: Eigenfrequency 1
1 In the Model Builder window, expand the Study 2 node, then click Step 1:
Eigenfrequency 1.
2 In the Settings window for Eigenfrequency, locate the Physics and Variables Selection
section.
3 In the table, enter the following settings:
Physics interface
Solve for
Discretization
physics
physics
17 |
Again, add a surface plot of the surface losses and an arrow plot of the surface current
(Figure 2, middle).
3D Plot Group 4
1 On the Model toolbar, click Add Plot Group and choose 3D Plot Group.
2 In the Settings window for 3D Plot Group, locate the Data section.
3 From the Data set list, choose Study 2/Parametric Solutions 2 (6).
4 Right-click Results>3D Plot Group 4 and choose Surface.
5 In the Settings window for Surface, click Replace Expression in the upper-right corner
18 |
12 Click OK.
19 |
ROOT
1 On the Model toolbar, click Add Physics to open the Add Physics window.
2 Go to the Add Physics window.
3 In the Add physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
4 Click Add to Component in the window toolbar.
5 On the Model toolbar, click Add Physics to close the Add Physics window.
ADD STUDY
1 On the Model toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Custom Studies>Empty
Study.
4 Click Add Study in the window toolbar.
5 On the Model toolbar, click Add Study to close the Add Study window.
GEOMETRY 3
Sphere 1 (sph1)
1 On the Geometry toolbar, click Sphere.
2 In the Settings window for Sphere, locate the Size section.
3 In the Radius text field, type a_s.
20 |
Material 5 (mat5)
1 In the Model Builder window, under Component 3 (comp3) right-click Materials and
21 |
Name
Value
Unit
Property group
Relative permittivity
epsilonr
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Material 6 (mat6)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 From the Selection list, choose All boundaries.
5 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
sigma_wall
S/m
Basic
Integration 5 (intop5)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_v in the Operator name text field.
3 Locate the Source Selection section. From the Selection list, choose All domains.
22 |
Integration 6 (intop6)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_s in the Operator name text field.
3 Locate the Source Selection section. From the Geometric entity level list, choose
Boundary.
4 From the Selection list, choose All boundaries.
Variables 3a
1 On the Definitions toolbar, click Local Variables.
2 In the Settings window for Variables, locate the Variables section.
3 Click Load from File.
4 Browse to the models Model Library folder and double-click the file
cavity_resonators_model3_variables.txt.
The emw3. prefix in the imported table is for the physics interface, Electromagnetic
Waves, Frequency Domain, in the third model.
MESH 3
Size
1 In the Model Builder window, under Component 3 (comp3) right-click Mesh 3 and
type h_max_s/d_f.
5 In the Maximum element growth rate text field, type 2.
6 In the Curvature factor text field, type 1.
7 In the Resolution of narrow regions text field, type 0.1.
8 Click the Build All button.
STUDY 3
1 Select both Study 2> Step 1: Eigenfrequency 1, Study 2>Parametric Sweep 1 using
Step 1: Eigenfrequency 1
1 In the Model Builder window, under Study 3 click Step 1: Eigenfrequency 1.
23 |
2 In the Settings window for Eigenfrequency, locate the Physics and Variables Selection
section.
3 In the table, enter the following settings:
Physics interface
Solve for
Discretization
physics
physics
physics
3D Plot Group 6
1 On the Model toolbar, click Add Plot Group and choose 3D Plot Group.
24 |
2 In the Settings window for 3D Plot Group, locate the Data section.
3 From the Data set list, choose Study 3/Parametric Solutions 3 (12).
4 Right-click Results>3D Plot Group 6 and choose Surface.
5 In the Settings window for Surface, click Replace Expression in the upper-right corner
25 |
Derived Values
Finish by evaluating the Q-factor and resonant frequency. Compare them with those
values in Table 1, Table 2 and Table 3.
1 On the Results toolbar, click Global Evaluation.
2 In the Settings window for Global Evaluation, locate the Data section.
3 From the Data set list, choose Study 1/Parametric Solutions 1.
4 From the Eigenfrequency selection list, choose First.
5 From the Table columns list, choose Inner solutions.
6 Click Replace Expression in the upper-right corner of the Expression section. From
26 |
27 |
28 |
Connecting a 3D Electromagnetic
Wave Model to an Electrical Circuit
Introduction
A model built with the RF Module can be connected to an electrical circuit equivalent,
if there is some structure outside of the model space that you wish to approximate as
a circuit equivalent. An example is shown in Figure 1, the 3D model of a coaxial cable
is connected to a voltage source, in series with a matched impedance, and sees a load,
also of matched impedance.
Zsource
FE model of a coaxial cable
VAC
Zload
Model Definition
The geometry in this model is a short section of a air-filled coaxial transmission line,
shown schematically in Figure 1. A 3D modeling space is used to model the coaxial
cable. The walls of the coax are treated as perfect electric conductors. This is
appropriate when the skin depth, and the losses in the conductors, are insignificant.
At one end of the coaxial cable, Lumped Port boundary condition is used to connect
the model to nodes 0 and 2 of the Electrical Circuit. A Voltage Source between circuit
nodes 0 and 1 excites the system, and a Resistor representing the source impedance is
added between nodes 1 and 2. Node 0 is specified as the Ground Node by default,
1 |
C O N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O A N E L E CT R I C A L CI RC U I T
which fixes the absolute voltage. The connection from the Electrical Circuit model to
the Electromagnetic Waves interface is via the External I Vs. U features.
At the other end of the coaxial cable, another Lumped Port boundary condition is used
to connect the model to nodes 3 and 0 of the Electrical Circuit. A Resistor which works
as a matched load is added between nodes 3 and 0. At any non-zero frequency, the
absolute voltage has no well-defined meaning, voltage only has a meaning as the path
integral of electric field between two points, so any arbitrary point in the model can be
chosen to have zero voltage. If you are working with a purely RF model, without an
electrical circuit, it is not even possible to fix the absolute voltage. However, when
using the Electrical Circuit interface, it requires that the absolute voltage be fixed at
one node (Node 0) in the model.
When solving such a model, some changes to the solver settings may be needed, since
the default solver suggestions may not be appropriate. The appropriate way to solve
this problem is with a Fully Coupled solver, using the default Iterative solver. Follow
the Modeling Instructions for the solver settings.
2 |
CO N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O AN E L E C T R I CA L CI RC U I T
Figure 2: Electric field norm (multislices) and magnetic field, and power flow (green,
blue arrows) inside the coaxial cable.
The fields and power flow plot shows the TEM wave propagation inside the coaxial
cable, which is excited by the Electrical Circuit interface.
3 |
C O N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O A N E L E CT R I C A L CI RC U I T
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
4 |
CO N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O AN E L E C T R I CA L CI RC U I T
Expression
Value
Description
r_coax
1[mm]
0.0010000 m
R_coax
2[mm]
0.0020000 m
L_coax
20[mm]
0.020000 m
f0
10[GHz]
1.0000E10 Hz
Frequency
lda0
c_const/f0
0.029979 m
Wavelength, air
h_max
0.2*lda0
0.0059958 m
Z_coax
Z0_const/
(2*pi)*log(R_coax/
r_coax)
41.560
Analytical impedance
Here, c_const and Z0_const are predefined COMSOL constants for the speed of
the light and the wave impedance in vacuum, respectively.
GEOMETRY 1
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type R_coax.
4 In the Height text field, type L_coax.
5 Locate the Axis section. From the Axis type list, choose x-axis.
6 Right-click Component 1 (comp1)>Geometry 1>Cylinder 1 (cyl1) and choose Rename.
7 In the Rename Cylinder dialog box, type Coax outer in the New label text field.
8 Click OK.
Cylinder 2 (cyl2)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type r_coax.
4 In the Height text field, type L_coax.
5 Locate the Axis section. From the Axis type list, choose x-axis.
5 |
C O N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O A N E L E CT R I C A L CI RC U I T
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object cyl1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the object cyl2 only.
6 Click the Build All Objects button.
6 |
CO N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O AN E L E C T R I CA L CI RC U I T
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 1 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Type of lumped port list, choose Coaxial.
5 From the Terminal type list, choose Circuit.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
7 |
C O N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O A N E L E CT R I C A L CI RC U I T
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Type of lumped port list, choose Coaxial.
5 From the Terminal type list, choose Circuit.
ELECTRICAL CIRCUIT (CIR)
Voltage Source 1
1 On the Physics toolbar, click Voltage Source.
2 In the Settings window for Voltage Source, locate the Node Connections section.
3 In the table, enter the following settings:
Label
Node names
4 Locate the Device Parameters section. From the Source type list, choose AC-source.
Resistor 1
1 On the Physics toolbar, click Resistor.
2 In the Settings window for Resistor, locate the Node Connections section.
8 |
CO N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O AN E L E C T R I CA L CI RC U I T
Node names
4 Locate the Device Parameters section. In the R text field, type Z_coax.
Resistor 2
1 On the Physics toolbar, click Resistor.
2 In the Settings window for Resistor, locate the Node Connections section.
3 In the table, enter the following settings:
Label
Node names
4 Locate the Device Parameters section. In the R text field, type Z_coax.
External I Vs. U 1
1 On the Physics toolbar, click External I Vs. U.
2 In the Settings window for External I Vs. U, locate the Node Connections section.
3 In the table, enter the following settings:
Label
Node names
4 Locate the External Device section. From the V list, choose Lumped port voltage
(emw/lport1).
External I Vs. U 2
1 On the Physics toolbar, click External I Vs. U.
2 In the Settings window for External I Vs. U, locate the Node Connections section.
3 In the table, enter the following settings:
Label
Node names
9 |
C O N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O A N E L E CT R I C A L CI RC U I T
4 Locate the External Device section. From the V list, choose Lumped port voltage
(emw/lport2).
MATERIALS
Next, assign material properties on the model. Specify air for the coaxial cable.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
MESH 1
10 |
C O N N E C T I N G A 3 D E L E C T RO M AG N E T I C WAVE M O D E L T O A N E L E CT R I C A L C I RC U IT
STUDY 1
Solution 1
1 On the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 1 node.
3 Right-click Stationary Solver 1 and choose Fully Coupled.
4 On the Study toolbar, click Compute.
RESULTS
11 |
C O N N E C T I N G A 3 D E L E C T RO M A G N E T I C WAVE M O D E L T O A N E L E CT R I C A L CI RC U I T
1 In the Model Builder window, under Results right-click Electric Field (emw) and
12 |
C O N N E C T I N G A 3 D E L E C T RO M AG N E T I C WAVE M O D E L T O A N E L E CT R I C A L C I RC U IT
Outer conductor
Dielectric
Center
conductor
Figure 1: Cross section of a coaxial cable. The arrows visualize the magnetic field.
Model Definition
Because a coax operates in TEM modewith the electric and magnetic fields normal
to the direction of propagation along the cablemodeling a 2D cross section suffices
to compute the fields and the impedance. For this model, assume perfect conductors
1 |
and a lossless dielectric with relative permittivity r = 2.4. The inner and outer radii are
0.5 mm and 3.43 mm, respectively.
The characteristic impedance, Z0 = V/I, of a transmission line relates the voltage
between the conductors to the current through the line. Although the model does not
involve computing the potential field, the voltage of the TEM waveguide can be
evaluated as a line integral of the electric field between the conductors:
ri
V = V i V o = E dr
(1)
ro
Similarly, the current is obtained as a line integral of the magnetic field along the
boundary of either conductor or any closed contour, C, bisecting the space between
the conductors:
I =
C H dr
(2)
The voltage and current in the direction out of the plane are positive for integration
paths oriented as in Figure 2.
Figure 2: The impedance of a coaxial cable can be found from the voltage, V, and current,
I, which are computed via line integrals as shown.
The value of Z0 obtained in this way, should be compared with the analytic result
ro
1 0
Z 0,analytic = ------ --------- log ----- 74.5
2 r
ri
0
2 |
Figure 3: Electric field magnitude (surface) and magnetic field (arrows) inside the
coaxial cable.
The impedance computed with the default mesh is Z0 = 74.65 . When the mesh is
refined, the result will approach the analytic value of 74.5 .
3 |
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
r_i
0.5[mm]
5.0000E-4 m
r_o
3.43[mm]
0.0034300 m
eps_r
2.4
2.4000
Relative
dielectric
constant
Z0_analytic
(Z0_const/
(2*pi*sqrt(eps_r)))*
log(r_o/r_i)
74.531
Characteristic
impedance,
analytic
4 |
GEOMETRY 1
Create the geometry using a single circle node with the radius of the outer conductor
and an extra layer representing the inner conductor.
Circle 1 (c1)
1 On the Geometry toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Object Type section.
3 From the Type list, choose Curve.
4 Locate the Size and Shape section. In the Radius text field, type r_o.
5 Click to expand the Layers section. In the table, enter the following settings:
Layer name
Thickness (m)
Layer 1
r_o-r_i
An advantage of using layers is that you automatically get a radial line to use for
computing the voltage as a line integral of the electric field.
7 Click the Zoom Extents button on the Graphics toolbar.
5 |
MATERIALS
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
Name
Value
Unit
Property group
Relative permittivity
epsilonr
eps_r
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Add a variable for the characteristic impedance computed as the voltage between the
conductors divided by the current through the cable. For this purpose, you need two
integration coupling operators.
Integration 1 (intop1)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_rad in the Operator name text
field.
3 Locate the Source Selection section. From the Geometric entity level list, choose
Boundary.
6 |
4 Select Boundary 1 only. Any of the four interior boundaries that connect the two
Integration 2 (intop2)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, type int_circ in the Operator name text
field.
3 Locate the Source Selection section. From the Geometric entity level list, choose
Boundary.
7 |
Now define the variable for the characteristic impedance computed from the
simulation.
Variables 1
1 On the Definitions toolbar, click Local Variables.
2 In the Settings window for Variables, locate the Variables section.
3 In the table, enter the following settings:
Name
Expression
Unit
Description
int_rad(-emw.Ex*t1x-emw.Ey*t1y)
Voltage
-int_circ(emw.Hx*t1x+emw.Hy*t1y)
Current
Z0_model
V/I
Characteristic
impedance
Here, t1x and t1y are the tangential vector components along the integration
boundaries ('1' refers to the boundary dimension). Shortly, you will determine the
tangential vector directions along the boundaries using an arrow plot of t1. The
signs in the definitions above are chosen such that V = V1 - V0 (see Equation 1) and
to have a positive current value correspond to a current in the positive z direction.
8 |
The emw. prefix gives the correct physics-interface scope for the electric and
magnetic field vector components.
ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN (EMW)
Keep the default physics settings, which include perfect electric conductor conditions
for the outer boundaries.
MESH 1
STUDY 1
9 |
RESULTS
10 |
A comparison with Equation 1 reveals that the line integral for the voltage computes
the potential difference V1 - V0. When computing the line integral for the current,
the clockwise orientation of the integration contour would mean that a positive
current is directed in the negative z direction, that is, into the modeling plane. The
minus sign added in the definition of I reverses this direction.
14 Right-click Results>Electric Field (emw)>Arrow Line 1 and choose Disable to retrieve
Derived Values
1 On the Results toolbar, click Global Evaluation.
2 In the Settings window for Global Evaluation, click Replace Expression in the
1 Go to the Table window. The result, roughly 74.65 , is within 0.2% of the analytic
value, 74.53 .
11 |
12 |
Model Definition
The model setup, schematically depicted in Figure 1, is a short section of an air-filled
coaxial transmission line. The symmetry of the structure allows for a 2D axisymmetric
model geometry.
Zref
Zref
V0(t)
Model space
1 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
At the other end of the coax, consider, in turn, three different boundary conditions:
1 perfect electric conductor (PEC)to simulate the short condition;
2 perfect magnetic conductor (PMC)to simulate an open condition; and
3 lumped portto simulate a matched load.
On the walls of the coax, apply a PEC boundary condition; this condition is
appropriate when both skin depth and losses in the conductors are very small.
Use a triangular mesh with the maximum element size chosen such that there are at
least two elements in the radial direction and at least eight elements per wavelength.
The only changes required to the default solver settings are to tighten the relative
tolerance from the default value, and to adjust the timespan and output time steps. The
internal time steps taken by the solver are auto-selected based on the specified relative
tolerance.
2 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
PEC
PMC
Matched load
Figure 2: Radial component of electric field at the input port versus time for three
different termination conditions: short (blue), open (green), and matched load (red).
Modeling Instructions
From the File menu, choose New.
NEW
3 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Time Dependent.
6 Click Done.
DEFINITIONS
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
r_coax
1[mm]
0.001000 m
R_coax
2[mm]
0.002000 m
L_coax
40[mm]
0.04000 m
Length of coax
core into cavity
20[GHz]
2.000E10 Hz
Pulse frequency
c_const/f
0.01499 m
Wavelength, free
space
1/f
5.000E-11 s
Period
h_max
min(L/8,(R_coax-r_coax)/2)
5.000E-4 m
Maximum element
size
Next, define the excitation, V0(t), in terms of a Gaussian pulse and a sine function.
text field.
3 Locate the Parameters section. In the Location text field, type 2*T.
4 In the Standard deviation text field, type T/2.
Analytic 1 (an1)
1 On the Model toolbar, click Functions and choose Global>Analytic.
2 In the Settings window for Analytic, type V0 in the Function name text field.
4 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
Lower limit
Upper limit
0.2[ns]
8 Click to collapse the Plot parameters section. Click the Plot button.
GEOMETRY 1
An elongated rectangle offset from the symmetry axis represents the straight coaxial
cable.
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type R_coax-r_coax.
4 In the Height text field, type L_coax.
5 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
Set up a point probe for plotting the electric field component Er while solving. You
will also use this plot to reproduce Figure 2.
1 On the Definitions toolbar, click Probes and choose Domain Point Probe.
2 In the Settings window for Domain Point Probe, locate the Point Selection section.
3 In row Coordinates, set r to r_coax.
4 Select the Snap to closest boundary check box.
5 In the Model Builder window, expand the Component 1 (comp1)>Definitions>Domain
Point Probe 1 node, then click Point Probe Expression 1 (ppb1).
6 In the Settings window for Point Probe Expression, locate the Expression section.
7 Click temw.Er - Electric field, r component in the upper-right corner of the section.
E L E C T RO M A G N E T I C WAVE S , TR A N S I E N T ( T E M W )
Now set up the physics. Begin by defining the Lumped port input condition.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 2 only (the bottom boundary).
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Wave excitation at this port list, choose On.
5 In the V0 text field, type V0(t).
6 Locate the Settings section. In the Zref text field, type (Z0_const/2/
pi)*log(R_coax/r_coax).
Finally, define a lumped port condition to use for the matched load case.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
6 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
1 On the Model toolbar, click Add Material to open the Add Material window.
ADD MATERIAL
1 On the Model toolbar, click Add Material to close the Add Material window
MESH 1
Size
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
type h_max.
7 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
STUDY 1
To study the short termination case first, disable the PMC and lumped port
conditions so that the default PEC condition is activated on the termination
boundary.
6 Locate the Physics and Variables Selection section. Select the Modify physics tree and
variables for study step check box.
7 In the Physics and Variables Selection section, under Component 1>Electromagnetic
Waves, Transient: Ctrl-click first Perfect Magnetic Conductor 1 and then Lumped Port
2 so both are selected. Then right-click and choose Disable (or use the Disable
button).
8 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
STUDY 1
Electric Field
Click on the Probe Plot 1 tab to place it in focus.
Electric Field
The default surface plot shows the electric field in the coax at the end of the simulation
interval. Because the transient has died out, the solution you see is only noise. Modify
the time to get a more interesting plot.
1 In the Model Builder window, under Results click Electric Field.
2 In the Settings window for 2D Plot Group, locate the Data section.
3 From the Time (s) list, choose 1.500E-10.
9 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
10 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
2 In the Settings window for Study, locate the Study Settings section.
3 Clear the Generate default plots check box.
4 On the Model toolbar, click Compute.
RESULTS
11 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
E L E C T RO M A G N E T I C WAVE S , TR A N S I E N T ( T E M W )
Lumped Port 2
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Transient (temw) right-click Lumped Port 2 and choose Enable. Note that you do not
need to disable the PMC condition because it is overridden by the lumped port.
STUDY 1
12 |
TR A N S I E N T M O D E L I N G O F A C O A X I A L C A B L E
Conical Antenna
Introduction
Conical antennas are useful for many applications due to their broadband
characteristics and relative simplicity. This example includes an analysis of the antenna
impedance and the radiation pattern as functions of the frequency for a monoconical
antenna with a finite ground plane and a 50 coaxial feed. The rotational symmetry
makes it possible to model this in axially symmetric 2D. When modeling in 2D, you
can use a dense mesh, giving an excellent accuracy for a wide range of frequencies.
Model Definition
The antenna geometry consists of a 0.2 m tall metallic cone with a top angle of
90 degrees on a finite ground plane of a 0.282 m radius. The coaxial feed has a central
conductor of 1.5 mm radius and an outer conductor (screen) of 4.916 mm radius
separated by a teflon dielectric of relative permittivity of 2.07. The central conductor
of the coaxial cable is connected to the cone, and the screen is connected to the ground
plane.
Figure 1: The geometry of the antenna. The central conductor of the coaxial cable is
connected to the metallic cone, and the cable screen is connected to the finite ground plane.
The model takes advantage of the rotational symmetry of the problem, which allows
modeling in 2D using cylindrical coordinates. You can then use a very fine mesh to
achieve an excellent accuracy.
1 |
CONICAL ANTENNA
DOMAIN EQUATIONS
BOUNDARY CONDITIONS
2 |
CONICAL ANTENNA
obtain maximum transmission into free space. This is quite well fulfilled in the high
frequency range.
3 |
CONICAL ANTENNA
Figure 3: Polar plot of the antenna radiation pattern in the near field versus the elevation
angle for 200 MHz, 800 MHz, and 1.5 GHz. The scale is logarithmic.
Figure 4: Polar plot of the antenna radiation pattern for the far field versus the elevation
angle for 200 MHz, 800 MHz, and 1.5 GHz. This plot is normalized differently but has
a shape resembling the near field.
4 |
CONICAL ANTENNA
As the frequency increases the antenna impedance gets closer to 50 , which means
that a voltage generator connected to the input of the antenna should have an output
impedance of 50 .
Modeling Instructions
From the File menu, choose New.
NEW
Import 1 (imp1)
1 On the Model toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 Click Browse.
4 Browse to the models Model Library folder and double-click the file
conical_antenna.mphbin.
5 |
CONICAL ANTENNA
5 Click Import.
The imported geometry is effectively a half circle with the metal areas removed. You
model the electromagnetic waves in the air and the dielectric material inside the
coaxial cable. There is no need to include the metal as a domain in the model
because the fields in it are essentially zero except on its surface.
DEFINITIONS
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
Z_tl
50[ohm]
50.000
Characteristic
transmission line
impedance
Variables 1
1 On the Model toolbar, click Variables and choose Local Variables.
2 In the Settings window for Variables, locate the Variables section.
6 |
CONICAL ANTENNA
Expression
Z_tl*(1+emw.S11)/
(1-emw.S11)
Unit
Description
Antenna impedance
Define the following selections in order to get easy access to some frequently used
domains and boundaries.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 1 and choose Rename.
3 In the Rename Explicit dialog box, type Air in the New label text field.
4 Click OK.
5 Select Domain 1 only.
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 2 and choose Rename.
3 In the Rename Explicit dialog box, type Dielectric in the New label text field.
4 Click OK.
The dielectric domain is inside the coaxial cable just below the cone. It is easier to
select it if you zoom in a little.
5 Select Domain 2 only.
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 3 and choose Rename.
3 In the Rename Explicit dialog box, type Outer Air Boundaries in the New label
text field.
4 Click OK.
5 In the Settings window for Explicit, locate the Input Entities section.
6 From the Geometric entity level list, choose Boundary.
7 |
CONICAL ANTENNA
With all selections and expressions now defined, it is time to set up the materials and
the physics of the model.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
Air (mat1)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Air (mat1).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Air.
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Dielectric.
4 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
2.07
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
Set up Boundary 6, at the bottom of the coaxial cable, to be a Port. You can zoom
in on this part of the geometry to easier find and select this boundary.
2 Select Boundary 6 only.
3 In the Settings window for Port, locate the Port Properties section.
8 |
CONICAL ANTENNA
The Scattering boundary condition is a simple way of letting the waves undergo only
minor artificial reflections as they leave the computational domain through the
exterior boundaries. To minimize these reflections, but at a greater computational
cost, you can use Perfectly Matched Layers.
Far-Field Domain 1
1 On the Physics toolbar, click Domains and choose Far-Field Domain.
2 Select Domain 1 only.
Far-Field Calculation 1
1 In the Model Builder window, expand the Far-Field Domain 1 node, then click Far-Field
Calculation 1.
2 In the Settings window for Far-Field Calculation, locate the Boundary Selection
section.
3 From the Selection list, choose Outer Air Boundaries.
Adding a Far Field Calculation feature does not affect the physics of the model, but
makes it possible to study the far field generated by the antenna. Select the
boundaries to use for this computation so that, in the physical (3D) geometry, they
surround all sources and reflecting objects. The outer air boundaries are a
convenient choice.
9 |
CONICAL ANTENNA
MESH 1
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
Size
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Size.
2 In the Settings window for Size, locate the Element Size section.
3 Click the Custom button.
4 Locate the Element Size Parameters section. In the Maximum element size text field,
type 25[mm].
This global maximum element size makes sure that the mesh everywhere resolves
the wavelength.
Size 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Free
Triangular 1 and choose Size.
2 In the Settings window for Size, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Domain.
4 From the Selection list, choose Dielectric.
5 Locate the Element Size section. Click the Custom button.
6 Locate the Element Size Parameters section. Select the Maximum element size check
box.
7 In the associated text field, type 0.5[mm].
You now have an especially fine mesh inside the coaxial cable, where the wave is
produced.
Size 2
1 Right-click Free Triangular 1 and choose Size.
2 In the Settings window for Size, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 Select Boundaries 4 and 8 only.
10 |
CONICAL ANTENNA
box.
7 In the associated text field, type 2.5[mm].
STUDY 1
The frequency range you just entered runs from 200 MHz to 1.5 GHz in steps of
25 MHz.
4 On the Model toolbar, click Compute.
11 |
CONICAL ANTENNA
RESULTS
12 |
CONICAL ANTENNA
1D Plot Group 5
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 On the 1D plot group toolbar, click Global.
3 In the Settings window for Global, locate the y-Axis Data section.
4 In the table, enter the following settings:
Expression
Unit
Description
real(Z)
Resistance
imag(Z)
Reactance
5 Click to expand the Coloring and style section. Locate the Coloring and Style section.
Find the Line style subsection. From the Line list, choose Cycle.
6 On the 1D plot group toolbar, click Plot.
You have now reproduced Figure 2. Next, visualize the near-field and far-field
radiation patterns using polar plots.
Select a few of the frequencies from the list of parameter values. Showing the
radiation pattern for all of them would take a bit of time and lead to a cluttered plot.
2 In the Settings window for Polar Plot Group, locate the Data section.
3 From the Parameter selection (freq) list, choose From list.
4 In the Parameter values (freq) list, choose 2.0000E8, 8.0000E8, and 1.5000E9.
13 |
CONICAL ANTENNA
9 Locate the r-Axis Data section. In the Expression text field, type
10*log10(emw.nPoav).
The variable emw.nPoav represents the outgoing power flow through the
boundaries where it is evaluated. The expression you just entered gives you the same
in a logarithmic scale.
10 Select the Description check box.
11 In the associated text field, type Near-field radiation pattern.
12 Locate the Angle Data section. From the Parameter list, choose Expression.
13 In the Expression text field, type atan2(z,r).
14 Click to expand the Coloring and style section. Locate the Coloring and Style section.
Find the Line style subsection. From the Line list, choose Cycle.
15 Click to collapse the Coloring and style section. Click to expand the Legends section.
18 Click to collapse the Legends section. Click the Zoom Extents button on the Graphics
toolbar.
Your near-field radiation plot should look like that in Figure 3.
Finally, visualize the far-field radiation pattern.
corner of the r-axis data section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Far field>emw.normdBEfar - Far-field
norm, dB.
4 Locate the r-Axis Data section. In the Description text field, type Far-field
radiation pattern.
14 |
CONICAL ANTENNA
15 |
CONICAL ANTENNA
16 |
CONICAL ANTENNA
Metallic housing
Polytetrafluoroethylene
(PTFE) ring
Cylindrical cavity
Dielectric
substrate
Figure 1: An evanescent mode cavity filter. The signal fed from a microstrip line is slot
coupled into the cylindrical cavity loaded with a PTFE ring.
Model Definition
The resonant frequency of the empty cylindrical waveguide cavity TE111 mode can be
calculated from the equation
p' nm 2
c
l 2
f nml = ---------------------- ------------ + -----
a
d
2 r r
1 |
where a and d are the radius and length of the cylinder, respectively, and p'nm is the
mth root of the Bessel function J'n(x). The TE111 mode is the dominant TE mode of
the cylindrical cavity resonator, and for a cavity of 25 mm radius and 100 mm height
this resonance is at 3.823 GHz. The starting point of this model was a computation
(not presented here) of the TE111 mode resonant frequency of an empty cylindrical
cavity and a subsequent verification of agreement with the analytic solution.
This basic model was then modified by the addition of a metal box at either end
representing a housing. Inside is a dielectric substrate and a microstrip line which is slot
coupled into the cavity. This represents the input and output of the device.
The slots are located on the center of the cavity ends to induce symmetric fields and
they are also parallel to each other to couple the injected fields maximally. The size of
the slots are tuned to provide a better matching to the reference characteristic
impedance assigned on ports. The model uses lumped ports to excite the structure.
The end of each microstrip line over the slots is shorted to couple the fields from the
microstrip lines through the slots and vice versa. The cavity is partially filled with a ring
of PTFE, r = 2.1, which causes the resonant frequency to shift down.
2 |
Figure 2: The frequency response of the filter shows bandpass filter characteristics. The
center frequency is lower than the dominant mode resonant frequency of the metallic
cavity.
Figure 3: The dielectric tube inside the cavity does not distort the electric field distribution
at resonance significantly.
3 |
Reference
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
4 |
Expression
Value
Description
60[mil]
0.0015240 m
Substrate thickness
l_slot
42[mm]
0.042000 m
Slot length
w_slot
3[mm]
0.0030000 m
Slot width
f_min
3.45[GHz]
3.4500E9 Hz
Minimum frequency in
sweep
f_max
3.61[GHz]
3.6100E9 Hz
Maximum frequency in
sweep
lda_min
c_const/f_max
0.083045 m
h_max
0.2*lda_min
0.016609 m
Here 'mil' refers to the unit milliinches, that is 1 mil = 0.0254 mm.
GEOMETRY 1
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type 25.
4 In the Height text field, type 100.
5 Click the Build Selected button.
6 Right-click Component 1 (comp1)>Geometry 1>Cylinder 1 (cyl1) and choose Rename.
7 In the Rename Cylinder dialog box, type Cavity in the New label text field.
8 Click OK.
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
5 |
Create a substrate.
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 60.
4 In the Depth text field, type 60.
5 In the Height text field, type d.
6 Locate the Position section. From the Base list, choose Center.
7 In the z text field, type -d/2.
8 Click the Build Selected button.
9 Right-click Component 1 (comp1)>Geometry 1>Block 1 (blk1) and choose Rename.
10 In the Rename Block dialog box, type Bottom_plate in the New label text field.
11 Click OK.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 3.2.
4 In the Depth text field, type 25.
5 In the Height text field, type d.
6 Locate the Position section. From the Base list, choose Center.
7 In the y text field, type 25/2-w_slot/2.
8 In the z text field, type -d/2.
9 Click the Build Selected button.
10 Right-click Component 1 (comp1)>Geometry 1>Block 2 (blk2) and choose Rename.
11 In the Rename Block dialog box, type Bottom_feed in the New label text field.
6 |
12 Click OK.
Block 3 (blk3)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 60.
4 In the Depth text field, type 60.
5 In the Height text field, type 20.
6 Locate the Position section. From the Base list, choose Center.
7 In the z text field, type -10.
8 Click the Build Selected button.
9 Click the Zoom Extents button on the Graphics toolbar.
10 Right-click Component 1 (comp1)>Geometry 1>Block 3 (blk3) and choose Rename.
11 In the Rename Block dialog box, type Housing in the New label text field.
12 Click OK.
13 Click the Wireframe Rendering button on the Graphics toolbar to see the interior.
7 |
Rotate 1 (rot1)
1 On the Geometry toolbar, click Transforms and choose Rotate.
2 Select all objects except for the cylinder, that is, blk3, wp1, blk1, and blk2.
3 In the Settings window for Rotate, locate the Rotation Angle section.
4 In the Rotation text field, type 0,180.
5 Locate the Point on Axis of Rotation section. In the z text field, type 50.
6 Locate the Axis of Rotation section. From the Axis type list, choose Cartesian.
7 In the x text field, type 1.
8 In the z text field, type 0.
9 Click the Build Selected button.
10 Click the Zoom Extents button on the Graphics toolbar.
Cylinder 2 (cyl2)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type 25.
4 In the Height text field, type 20.
5 Locate the Position section. In the z text field, type 40.
Cylinder 3 (cyl3)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type 15.
4 In the Height text field, type 20.
5 Locate the Position section. In the z text field, type 40.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object cyl2 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the object cyl3 only.
8 |
The default boundary condition is perfect electric conductor, which applies to all
exterior boundaries. Assign a perfect electric conductor condition to the remaining
boundaries of the cavity.
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic
Waves, Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
9 |
You can do this most easily by copying the text '21, 28, 35, and 42', clicking in the
selection box, and then pressing Ctrl+V, or by using the Paste Selection dialog box.
10 |
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 44 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Wave excitation at this port list, choose On.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 34 only.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
11 |
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Select Domains 2, 3, 7, and 8 only.
3 In the Settings window for Material, locate the Material Contents section.
4 In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
3.38
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Material 3 (mat3)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Select Domain 6 only.
3 In the Settings window for Material, locate the Material Contents section.
4 In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
2.1
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
12 |
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
STUDY 1
13 |
S-Parameter (emw)
1 In the Model Builder window, under Results click S-Parameter (emw).
2 In the Settings window for 1D Plot Group, click to expand the Legend section.
3 From the Position list, choose Lower right.
14 |
Di po le Ant e nna
Introduction
The dipole antenna is one of the most straightforward antenna configurations. It can
be realized with two thin metallic rods that have a sinusoidal voltage difference applied
between them. The length of the rods is chosen such that they are quarter wavelength
elements at the operating frequency. Such an antenna has a well-known torus-like
radiation pattern.
Conductive radiator
PML
/4
Lumped port
Figure 1: A dipole antenna. The model consists of two cylindrical arms of conductive
material with a voltage source in between. A region of free space bounded by a perfectly
matched layer (PML) surrounds the antenna.
Model Definition
The model of the antenna consists of two cylinders representing each of the dipole
arms. The free space wavelength at the antennas operating frequency is 4 m. Thus,
each of the antenna arms is 1 m long and aligned with the z-axis. The arm radius is
chosen to be 0.05 m. In the limit as the radius approaches zero, this antenna will
approach the analytic solution.
1 |
DIPOLE ANTENNA
A small cylindrical gap of size 0.01 m between the antenna arms represents the voltage
source. The power supply and feed structure are not modeled explicitly, and it is
assumed that a uniform voltage difference is applied across these faces. This source
induces electromagnetic fields and surface currents on the adjacent conductive faces.
The dipole arm surfaces are modeled using the Impedance Boundary Condition,
which is appropriate for conductive surfaces that have dimensions much larger than the
skin depth. This boundary condition introduces a finite conductivity at the surface as
well as resistive losses.
The air domain around the antenna is modeled as sphere of free space of radius 2 m,
which is approximately the boundary between the near-field and the far-field. This
sphere of air is truncated with a perfectly matched layer (PML) that acts as an absorber
of outgoing radiation. The far-field pattern is computed on the boundary between the
air and the PML domains.
The mesh is manually adjusted such that there are five elements per free space
wavelength and that the boundaries of the antenna are meshed more finely. The PML
is swept with a total of five elements along the radial direction.
2 |
DIPOLE ANTENNA
Figure 2: A slice plot of the electric field magnitude around the antenna.
Figure 3: The polar plot of the far field pattern in the xy-plane is isotropic.
3 |
DIPOLE ANTENNA
Figure 4: A 3D visualization of the far-field pattern of the dipole shows the expected
torus-shaped pattern.
Modeling Instructions
From the File menu, choose New.
NEW
4 |
DIPOLE ANTENNA
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
lambda0
4[m]
4.0000 m
Operating wavelength
arm_lengt
h
lambda0/4
1.0000 m
r_antenna
arm_length/20
0.050000 m
gap_size
arm_length/100
0.010000 m
GEOMETRY 1
Create a sphere with a layer. The outer layer presents the PML.
Sphere 1 (sph1)
1 On the Geometry toolbar, click Sphere.
2 In the Settings window for Sphere, locate the Size section.
3 In the Radius text field, type 2.5*arm_length.
4 Click to expand the Layers section. In the table, enter the following settings:
Layer name
Thickness (m)
Layer 1
0.5*arm_length
Then, add a cylinder with layers. The top and bottom parts are the antenna radiators.
A small gap between the antenna radiators is for the voltage source.
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
5 |
DIPOLE ANTENNA
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type r_antenna.
4 In the Height text field, type 2*arm_length+gap_size.
5 Locate the Position section. In the z text field, type -(arm_length+gap_size/2).
6 Click to expand the Layers section. In the table, enter the following settings:
Layer name
Thickness (m)
Layer 1
arm_length
The domain inside the antenna radiators is not part of the model analysis.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object sph1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the object cyl1 only.
6 |
DIPOLE ANTENNA
DEFINITIONS
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, locate the Input Entities section.
3 From the Geometric entity level list, choose Boundary.
4 Click Paste Selection.
5 In the Paste Selection dialog box, type 13-15, 18-20, 28, 30, 39, 41 in the
Selection text field.
7 |
DIPOLE ANTENNA
6 Click OK.
8 |
DIPOLE ANTENNA
View 1
Suppress some domains and boundaries. This helps to see the interior parts when
setting up the physics and reviewing the mesh.
1 On the 3D view toolbar, click Hide Geometric Entities.
2 Select Domains 1 and 2 only.
9 |
DIPOLE ANTENNA
Set up the physics for the model. Add an Impedance Boundary Condition that
overrides the default PEC boundary condition on the antenna radiator surface.
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic
Waves, Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
10 |
DIPOLE ANTENNA
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Click the Zoom In button on the Graphics toolbar in a couple of times to see the small
4 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
5 From the Type of lumped port list, choose User defined.
6 In the hport text field, type gap_size.
7 In the wport text field, type 2*pi*r_antenna.
8 Specify the ah vector as
0
11 |
DIPOLE ANTENNA
Far-Field Calculation 1
1 In the Model Builder window, expand the Far-Field Domain 1 node, then click Far-Field
Calculation 1.
2 In the Settings window for Far-Field Calculation, locate the Boundary Selection
section.
3 Click Clear Selection.
4 Click Paste Selection.
5 In the Paste Selection dialog box, type 9-12, 26-27, 34, 37 in the Selection text
field.
6 Click OK.
MATERIALS
Assign air as the material for all domains and override the antenna radiator surface with
copper.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
ADD MATERIAL
Copper (mat2)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Copper
(mat2).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 From the Selection list, choose Antenna.
5 On the Model toolbar, click Add Material to close the Add Material window.
12 |
DIPOLE ANTENNA
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
STUDY 1
13 |
DIPOLE ANTENNA
The results show the E-field norm distribution on the antenna radiators. It is plotted
in Figure 2.
3D Plot Group 3
1 In the Model Builder window, under Results click 3D Plot Group 3.
2 Click the Zoom Extents button on the Graphics toolbar.
Derived Values
Finish the result analysis by evaluating the port impedance.
1 On the Results toolbar, click Global Evaluation.
2 In the Settings window for Global Evaluation, click Replace Expression in the
14 |
DIPOLE ANTENNA
Fresnel Equations
Introduction
A plane electromagnetic wave propagating through free space is incident at an angle
upon an infinite dielectric medium. This model computes the reflection and
transmission coefficients and compares the results to the Fresnel equations.
Model Definition
A plane wave propagating through free space (n = 1) as shown in Figure 1 is incident
upon an infinite dielectric medium (n = 1.5) and is partially reflected and partially
transmitted. If the electric field is p-polarizedthat is, if the electric field vector is in
the same plane as the Poynting vector and the surface normalthen there will be no
reflections at an incident angle of roughly 56, known as the Brewster angle.
Reflected
Transmitted
Unit cell
n1
n2
Incident
Figure 1: A plane wave propagating through free space incident upon an infinite
dielectric medium.
Although, by assumption, space extends to infinity in all directions, it is sufficient to
model a small unit cell, as shown in Figure 1; a Floquet-periodic boundary condition
applies on the top and bottom unit-cell boundaries because the solution is periodic
along the interface. This model uses a 3D unit cell, and applies perfect electric
conductor and perfect magnetic conductor boundary conditions as appropriate to
1 |
FRESNEL EQUATIONS
model out-of-plane symmetry. The angle of incidence ranges between 090 for both
polarizations.
For comparison, Ref. 1 and Ref. 2 provide analytic expressions for the reflectance and
transmittance. Reflection and transmission coefficients for s-polarization and
p-polarization are defined respectively as
n 1 cos incident n 2 cos transmitted
r s = -----------------------------------------------------------------------------------n 1 cos incident + n 2 cos transmitted
(1)
2n 1 cos incident
t s = -----------------------------------------------------------------------------------n 1 cos incident + n 2 cos transmitted
(2)
(3)
2n 1 cos incident
t p = -----------------------------------------------------------------------------------n 1 cos transmitted + n 2 cos incident
(4)
n 2 cos transmitted 2
T = ------------------------------------------ t
n 1 cos incident
(5)
(6)
2 |
FRESNEL EQUATIONS
(7)
Figure 2: Electric field, Ey (slice) and power flow (arrows) for TE incidence at 70 inside
the unit cell.
3 |
FRESNEL EQUATIONS
Figure 3: Magnetic field, Hy (slice) and power flow (arrows) for TM incidence at 70
inside the unit cell.
Note that the sum of reflectance and transmittance in Figure 4 and Figure 5 equals 1,
showing conservation of power. Figure 5 also shows that the reflectance around 56
the Brewster angle in the TM caseis close to zero.
4 |
FRESNEL EQUATIONS
Figure 4: The reflectance and transmittance for TE incidence agree well with the analytic
solutions.
Figure 5: The reflectance and transmittance for TM incidence agree well with the analytic
solutions. The Brewster angle is also observed at the expected location.
5 |
FRESNEL EQUATIONS
References
1. C.A. Balanis, Advanced Engineering Electromagnetics, Wiley, 1989.
2. B.E.A. Saleh and M.C. Teich, Fundamentals of Photonics, Wiley, 1991.
Modeling Instructions
From the File menu, choose New.
NEW
Define some parameters that are useful when setting up the mesh and the study.
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
6 |
FRESNEL EQUATIONS
Expression
Value
Description
n_air
1.0000
n_slab
1.5
1.5000
lda0
1[m]
1.0000 m
Wavelength
f0
c_const/lda0
2.9979E8 1/s
Frequency
alpha
70[deg]
1.2217 rad
Angle of incidence
beta
asin(n_air*sin(al
pha)/n_slab)
0.67701 rad
Refraction angle
h_max
lda0/6
0.16667 m
alpha_bre
wster
atan(n_slab/
n_air)
0.98279 rad
r_s
(n_air*cos(alpha)
-n_slab*cos(beta)
)/
(n_air*cos(alpha)
+n_slab*cos(beta)
)
-0.54735
Reflection coefficient,
TE
r_p
(n_slab*cos(alpha
)-n_air*cos(beta)
)/
(n_air*cos(beta)+
n_slab*cos(alpha)
)
-0.20613
Reflection coefficient,
TM
t_s
(2*n_air*cos(alph
a))/
(n_air*cos(alpha)
+n_slab*cos(beta)
)
0.45265
Transmission
coefficient, TE
t_p
(2*n_air*cos(alph
a))/
(n_air*cos(beta)+
n_slab*cos(alpha)
)
0.52925
Transmission
coefficient, TM
The angle of incidence is updated while running the parametric sweep. The
refraction (transmitted) angle is defined by Snell's law with the updated angle of
incidence. The Brewster angle exists only for TM incidence, p-polarization, and
parallel polarization.
7 |
FRESNEL EQUATIONS
Variables 1
1 On the Model toolbar, click Variables and choose Local Variables.
2 In the Settings window for Variables, locate the Variables section.
3 In the table, enter the following settings:
Name
Expression
Unit
Description
ka
emw.k0
rad/m
kax
ka*sin(alpha)
rad/m
kay
kaz
ka*cos(alpha)
rad/m
kb
n_slab*emw.k0
rad/m
kbx
kb*sin(beta)
rad/m
kby
kbz
kb*cos(beta)
rad/m
GEOMETRY 1
First, create a block composed of two domains. Use layers to split the block.
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 0.2.
4 In the Depth text field, type 0.2.
5 In the Height text field, type 0.8.
6 Click to expand the Layers section. In the table, enter the following settings:
Layer name
Thickness (m)
Layer 1
0.4
8 |
FRESNEL EQUATIONS
Set up the physics based on the direction of propagation and the E-field polarization.
First, assume a TE-polarized wave which is equivalent to s-polarization and
perpendicular polarization. Ex and Ez are zero while Ey is dominant.
1 In the Model Builder window, expand the Component 1 (comp1)>Electromagnetic
Waves, Frequency Domain (emw) node, then click Electromagnetic Waves, Frequency
Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
9 |
FRESNEL EQUATIONS
exp(-i*kax*x)[V/m]
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 3 only.
3 In the Settings window for Port, locate the Port Mode Settings section.
4 Specify the E0 vector as
0
exp(-i*kbx*x)[V/m]
The bottom surface is an observation port. The S21-parameter from Port 1 and Port
2 provides the transmission characteristics.
The E-field polarization has Ey only and the boundaries are always either parallel or
perpendicular to the E-field polarization. Apply periodic boundary conditions on the
boundaries parallel to the E-field except those you already assigned to the ports.
Periodic Condition 1
1 On the Physics toolbar, click Boundaries and choose Periodic Condition.
2 Select Boundaries 1, 4, 10, and 11 only.
3 In the Settings window for Periodic Condition, locate the Periodicity Settings
section.
4 From the Type of periodicity list, choose Floquet periodicity.
10 |
FRESNEL EQUATIONS
11 |
FRESNEL EQUATIONS
Now set up the material properties based on refractive index. The top half is filled with
air.
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
12 |
Property
Name
Value
Unit
Property group
Refractive index
n_air
Refractive index
FRESNEL EQUATIONS
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Select Domain 1 only.
3 In the Settings window for Material, locate the Material Contents section.
4 In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Refractive index
n_slab
Refractive index
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
13 |
FRESNEL EQUATIONS
STUDY 1
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
alpha
range(0,2[deg],90[deg])
Parameter unit
field, type 0.
4 Find the z-planes subsection. In the Planes text field, type 0.
5 Locate the Coloring and Style section. From the Color table list, choose Wave.
6 In the Model Builder window, right-click Electric Field (emw) and choose Arrow
Volume.
14 |
FRESNEL EQUATIONS
7 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
that in Figure 2.
Add a 1D plot to see the reflection and transmission versus the angle of incidence.
1D Plot Group 2
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Plot Settings section.
3 Select the x-axis label check box.
4 In the associated text field, type Angle of Incidence.
5 Select the y-axis label check box.
6 In the associated text field, type Reflectance and Transmittance.
7 Click to expand the Legend section. From the Position list, choose Upper left.
8 On the 1D plot group toolbar, click Global.
9 In the Settings window for Global, locate the y-Axis Data section.
10 In the table, enter the following settings:
Expression
Unit
Description
abs(emw.S11)^2
Reflectance
abs(emw.S21)^2
Transmittance
11 Click to expand the Coloring and style section. Locate the Coloring and Style section.
Find the Line style subsection. From the Line list, choose None.
12 Find the Line markers subsection. From the Marker list, choose Cycle.
13 On the 1D plot group toolbar, click Global.
15 |
FRESNEL EQUATIONS
14 In the Settings window for Global, locate the y-Axis Data section.
15 In the table, enter the following settings:
Expression
Unit
Description
abs(r_s)^2
Reflectance, analytic
n_slab*cos(beta)/
(n_air*cos(alpha))*abs(t_s)^2
Transmittance, analytic
The remaining instructions are for the case of TM-polarized wave, p-polarization, and
parallel polarization. In this case, Ey is zero while Ex and Ez characterize the wave. In
other words, Hy is dominant while Hx and Hz are effectless. Thus, the H-field is
perpendicular to the plane of incidence and it is convenient to solve the model for the
H-field.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
Port 1
1 In the Settings window for Port, locate the Port Mode Settings section.
2 From the Input quantity list, choose Magnetic field.
3 Specify the H0 vector as
0
exp(-i*kax*x)[A/m]
Port 2
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw) click Port 2.
2 In the Settings window for Port, locate the Port Mode Settings section.
3 From the Input quantity list, choose Magnetic field.
16 |
FRESNEL EQUATIONS
exp(-i*kbx*x)[A/m]
Solution 1
1 In the Model Builder window, expand the Study 1>Solver Configurations node.
2 Right-click Solution 1 and choose Solution>Copy.
RESULTS
17 |
FRESNEL EQUATIONS
STUDY 1
click Multislice 1.
2 In the Settings window for Multislice, click Replace Expression in the upper-right
Unit
Description
abs(r_p)^2
Reflectance, analytic
n_slab*cos(beta)/
(n_air*cos(alpha))*abs(t_p)^2
Transmittance, analytic
4 On the 1D plot group toolbar, click Plot. The plot should look like Figure 5. The
Brewster angle is observed around 56 degrees, which is close to the analytic value.
18 |
FRESNEL EQUATIONS
H-Bend Waveguide 2D
Introduction
This example is a 2D version of H-Bend Waveguide 3D, which shows how to model a
bending rectangular waveguide for microwaves. For a general introduction, see the
model H-Bend Waveguide 3D.
The dimensions of the waveguide and the frequency range used in this model are such
that TE10 is the single propagating mode. In this mode, if the bend is in the xy-plane,
the electric field only has a z-component, which furthermore is independent of the
z-coordinate. This makes it possible to set up and solve the model in a 2D geometry.
Model Definition
The considered geometry is an xy-plane cross-section of the 3D geometry, as seen in
Figure 1. This figure also sums up the material and boundary settings, which are the
same as in the 3D model. The main advantage with setting up the model in 2D is that
it solves much faster and uses less memory. As a consequence, this version of the model
1 |
H- B E N D WAVE G UI D E 2 D
does not stress the need to adapt the mesh to the wavelength, but simply lets you apply
a mesh that is more than fine enough.
n=1.44
Ez=cos((/2)(y/a))
2 |
H- B E N D WAVE GU I D E 2 D
3 |
H- B E N D WAVE G UI D E 2 D
Modeling Instructions
From the File menu, choose New.
NEW
4 |
H- B E N D WAVE GU I D E 2 D
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GEOMETRY 1
Circle 1 (c1)
1 On the Geometry toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type 0.08.
Circle 2 (c2)
1 On the Geometry toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type 0.04.
Square 1 (sq1)
1 On the Geometry toolbar, click Primitives and choose Square.
2 In the Settings window for Square, locate the Size section.
3 In the Side length text field, type 0.08.
4 Locate the Position section. In the y text field, type -0.08.
Compose 1 (co1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Compose.
2 Click in the Graphics window and then press Ctrl+A to select all objects.
3 In the Settings window for Compose, locate the Compose section.
4 In the Set formula text field, type sq1*(c1-c2).
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 0.04.
4 In the Height text field, type 0.1.
5 Locate the Position section. In the x text field, type 0.04.
Rectangle 2 (r2)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
5 |
H- B E N D WAVE G UI D E 2 D
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
Name
Value
Unit
Property group
Refractive index
Refractive index
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Right-click Material 2 (mat2) and choose Rename.
3 In the Rename Material dialog box, type Silica Glass in the New label text field.
4 Click OK.
5 Select Domain 2 only.
6 In the Settings window for Material, click to expand the Material properties section.
7 Locate the Material Properties section. In the Material properties tree, select
Electromagnetic Models>Refractive Index>Refractive index (n).
6 |
H- B E N D WAVE GU I D E 2 D
Name
Value
Unit
Property group
Refractive index
1.44
Refractive index
With TE waves, only the z-component of the electric field needs to be solved for
1 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Components section.
2 From the Electric field components solved for list, choose Out-of-plane vector.
3 Locate the Physics-Controlled Mesh section. Select the Enable check box.
The default boundary condition is perfect electric conductor, which is fine for all
exterior boundaries except the ports. The software automatically imposes continuity
on interior boundaries.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 1 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
5 From the Wave excitation at this port list, choose On.
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
7 |
H- B E N D WAVE G UI D E 2 D
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
This gives you a range of frequencies from 4 GHz to 5.2 GHz, with a pitch of 0.025
GHz.
4 On the Model toolbar, click Compute.
RESULTS
S-Parameter (emw)
1 In the Model Builder window, expand the Results>S-Parameter (emw) node, then click
Global 1.
2 In the Settings window for Global, click to expand the Coloring and style section.
3 Locate the Coloring and Style section. Find the Line markers subsection. From the
Marker list, choose Cycle.
4 On the 1D plot group toolbar, click Plot.
8 |
H- B E N D WAVE GU I D E 2 D
H-Bend Waveguide 3D
Introduction
This example shows how to model a rectangular waveguide for microwaves. A single
hollow waveguide can conduct two kinds of electromagnetic waves: transversal
magnetic (TM) or transversal electric (TE) waves. This model examines a TE wave,
one that has no electric field component in the direction of propagation. More
specifically, for this model you select the frequency and waveguide dimension so that
TE10 is the single propagating mode. In that mode the electric field has only one
nonzero componenta sinusoidal with two nodes, one at each of the walls of the
waveguide. This makes it possible to set up and solve the model in 2D, which is done
in a separate version; see H-Bend Waveguide 2D.
One important design aspect is how to shape a waveguide to go around a corner
without incurring unnecessary losses in signal power. Unlike in wires, these losses
usually do not result from ohmic resistance but instead arise from unwanted
reflections. You can minimize these reflections by keeping the bend smooth with a
large enough radius. In the range of operation the transmission characteristics (the
ability of the waveguide to transmit the signal) must be reasonably uniform for
avoiding signal distortions.
With air as the inside medium of the waveguide, the transmission is nearly perfect
throughout the range of operation. In this model, to make the simulation and the
results more interesting, the bend will be filled with Silica glass, a dielectric medium.
The model also shows how to systematically compute and export all S-parameters to a
Touchstone file.
Model Definition
This example illustrates how to create a model that computes the electromagnetic
fields and transmission characteristics of a 90 bend for a given radius. This type of
waveguide bends changes the direction of the H field components and leaves the
direction of the E field unchanged. The waveguide is therefore called an H-bend. The
H-bend design used in this example is well-proven in real-world applications and you
can buy similar waveguide bends online from a number of manufacturers. This
particular bend performs optimally in the ideal case of perfectly conducting walls.
1 |
H- B E N D WAVE G UI D E 3 D
The waveguide walls are typically plated with a very good conductor, such as silver. In
this example the walls are considered to be made of a perfect conductor, which means
that the tangential component of the electric field is zero, or that n E = 0 on the
boundaries. This boundary condition is referred to as a perfect electric conductor
(PEC) boundary condition.
The geometry is as follows:
The waveguide is considered to continue indefinitely before and after the bend. This
means that the input wave needs to have the form of a wave that has been traveling
through a straight waveguide. The shape of such a wave is determined by the boundary
conditions of Maxwells equations on the sides of the metallic boundaries, that is, the
PEC boundary condition. If polarized according to a TE10 mode, the shape is known
analytically to be E = (0, 0, sin( (a y)/(2 a))) cos(t) given that the entrance
boundary is centered around the y = 0 axis, and that the width of the waveguide, in
the y direction, is 2a.
The model is set up using the time-harmonic electromagnetic waves physics interface.
This means that only the phasor component of the field is modeled. The incident field
then has the form E = (0, 0, E0z) = (0, 0, sin( (a y)/(2 a))), and is considered as part
2 |
H- B E N D WAVE GU I D E 3 D
PEC
Ez=cos((/2)(y/a))
The width of the waveguide is chosen so that it has a cutoff frequency of 3.7 GHz.
This makes the waveguide operational up to 7.5 GHz. At higher frequencies other
modes than the TE10 appear, causing a dirty signal. The input wave then splits into
several modes that are hard to control without having large power losses. Below the
cutoff frequency, no waves can propagate through the waveguide. This is an intrinsic
property of microwave waveguides.
The cutoff frequency of different modes in a straight waveguide is given by the relation
c m 2
n 2
( c ) mn = --- ----- + ---
2 a
b
where m and n are the mode numbers (m = 1, n = 0 for the TE10 mode), a and b are
the lengths of the sides of the waveguide cross-section, and c is the speed of light.
For this waveguide, a = 2b and b = 2 cm.
The first few cutoff frequencies are (c)10 = 3.7 GHz, (c)01 = 7.5 GHz,
(c)11 = 8.4 GHz. The frequencies used in this model are from 4.0 GHz to 5.2 GHz,
and hence entirely within the single-mode range.
3 |
H- B E N D WAVE G UI D E 3 D
On the input boundary, the Port boundary condition lets you choose which mode to
send in. Any reflected waves having the same shape will be transmitted back through
this same boundary. The output boundary also uses a Port condition, but without field
excitation, to specify the shape of the wave that it lets pass through. Using port
boundary conditions means that you automatically gain access to postprocessing
variables for the S-parameters.
Figure 1: The z-component of the electric field for a frequency of 5.1 GHz.
The S-parameters are shown as functions of the frequency in Figure 2.
4 |
H- B E N D WAVE GU I D E 3 D
Modeling Instructions
From the File menu, choose New.
NEW
5 |
H- B E N D WAVE G UI D E 3 D
Circle 1 (c1)
1 On the Work plane toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type 0.08.
Circle 2 (c2)
1 On the Work plane toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type 0.04.
Square 1 (sq1)
1 On the Work plane toolbar, click Primitives and choose Square.
2 In the Settings window for Square, locate the Size section.
3 In the Side length text field, type 0.08.
4 Locate the Position section. In the yw text field, type -0.08.
Compose 1 (co1)
1 On the Work plane toolbar, click Booleans and Partitions and choose Compose.
2 Click in the Graphics window and then press Ctrl+A to select all objects.
3 In the Settings window for Compose, locate the Compose section.
4 In the Set formula text field, type sq1*(c1-c2).
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 0.04.
4 In the Height text field, type 0.1.
6 |
H- B E N D WAVE GU I D E 3 D
Rectangle 2 (r2)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 0.1.
4 In the Height text field, type 0.04.
5 Locate the Position section. In the xw text field, type -0.1.
6 In the yw text field, type -0.08.
7 On the Work plane toolbar, click Build All.
You may click Close on the current Work Plane toolbar to access the Geometry toolbar.
Extrude 1 (ext1)
1 On the Geometry toolbar, click Extrude.
2 In the Settings window for Extrude, locate the Distances from Plane section.
3 In the table, enter the following settings:
Distances (m)
0.02
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
7 |
H- B E N D WAVE G UI D E 3 D
9 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Refractive index
Refractive index
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Right-click Material 2 (mat2) and choose Rename.
3 In the Rename Material dialog box, type Silica Glass in the New label text field.
4 Click OK.
5 Select Domain 2 only.
6 In the Settings window for Material, click to expand the Material properties section.
7 Locate the Material Properties section. In the Material properties tree, select
Electromagnetic Models>Refractive Index>Refractive index (n).
8 Click Add to Material.
9 Locate the Material Contents section. In the table, enter the following settings:
Property
Refractive index
Name
Value
Unit
Property group
1.44
Refractive index
1 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
2 Select the Enable check box.
8 |
H- B E N D WAVE GU I D E 3 D
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 1 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
5 From the Wave excitation at this port list, choose On.
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 15 only.
3 In the Settings window for Port, locate the Port Properties section.
4 In the Port name text field, type 2.
5 From the Type of port list, choose Rectangular.
The default boundary condition is perfect electric conductor, which is fine for all
exterior boundaries except the ports. The software automatically imposes continuity
on interior boundaries.
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
The range command you just entered means that you will get solutions for
frequencies from 4 GHz to 5.2 GHz, with a pitch of 25 MHz. The solution process
should only take a few minutes, but if you want to speed it up, you can increase the
pitch.
4 On the Model toolbar, click Compute.
9 |
H- B E N D WAVE G UI D E 3 D
RESULTS
S-Parameter (emw)
1 In the Model Builder window, expand the Results>S-Parameter (emw) node, then click
Global 1.
2 In the Settings window for Global, click to expand the Coloring and style section.
3 Locate the Coloring and Style section. Find the Line markers subsection. From the
Marker list, choose Cycle.
4 On the 1D plot group toolbar, click Plot.
The result, which should look like Figure 2, shows that the transmission varies
throughout the frequency range. Note in particular that S21 has two deep dips,
corresponding to almost perfect transmission. This is the result of resonances in the
bend. To confirm this, try looking at the field distribution for the frequency where
the upper peak is located, 5.1 GHz.
10 |
H-BEND WAVEGUIDE 3D
The standing wave pattern still remains in the bend, but at this frequency it is almost
completely gone in the input section.
For an alternative view, you can plot the instantaneous value of the electric field
inside the waveguide. Only the z-component will be substantially non-zero. For a
better view, add also deformation. Replace the Multislice with a single horizontal
slice plot.
5 In the Model Builder window, expand the Electric Field (emw) node.
6 Right-click Multislice 1 and choose Delete.
7 Right-click Electric Field (emw) and choose Slice.
8 In the Settings window for Slice, click Replace Expression in the upper-right corner
11 |
H- B E N D WAVE G UI D E 3 D
The Wave color table looks its best using a symmetric range. You can also play with
a deformed shape plot to make the waves appear more clearly.
13 Click to expand the Range section. Locate the Coloring and Style section. Select the
Symmetrize color range check box.
14 Right-click Results>Electric Field (emw)>Slice 1 and choose Deformation.
15 In the Settings window for Deformation, click Replace Expression in the upper-right
The remaining instructions show you how to systematically solve with one port active
at a time, and save the results in the Touchstone format.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
Click the Browse button and select a file to which you want to export the results. If
the file does not exist, it will be created.
DEFINITIONS
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
PortName
1.000
STUDY 1
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
12 |
H-BEND WAVEGUIDE 3D
Description
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
PortName
1 2
Parameter unit
The parameter is the same as the name suggested by the port sweep. The parameter
values should be the same as your port numbers.
5 On the Study toolbar, click Compute.
RESULTS
S-Parameter (emw)
As you can see, after performing the parametric sweep over the ports, the S-parameter
plot you created previously is empty. To restore the plot, you need to change the data
set and specify the inner parameter - that is, the frequency - as the quantity to display
along the horizontal axis.
1 In the Model Builder window, under Results>S-Parameter (emw) click Global 1.
2 In the Settings window for Global, locate the Data section.
3 From the Data set list, choose Study 1/Parametric Solutions 1.
To verify the reciprocity of the waveguide, you can add the S-parameters S12dB and
S22dB to the Expressions table and change the parameter selection for PortName:
4 Click Add Expression in the upper-right corner of the y-axis data section. From the
13 |
H- B E N D WAVE G UI D E 3 D
14 |
H-BEND WAVEGUIDE 3D
1 |
constructed from three rectangular waveguide sections joining at 120 and with a
ferrite post inserted at the center of the joint.
Port 3
Ferrite post
Port 2
Port 1
Figure 1: The post is magnetized by a static H0 bias field along its axis. The bias field is
supplied by external permanent magnets which are not explicitly modeled in this model.
IMPEDANCE MATCHING
An important step in the design of any microwave device is to match its input
impedance for a given operating frequency. Impedance matching is equivalent to
minimizing the reflections back to the inport. The parameters that need to be
determined are the size of the ferrite post and the width of the wider waveguide section
surrounding the ferrite. In this tutorial, these are varied in order to minimize the
reflectance. The scattering parameters (S-parameters) used as measures of the
reflectance and transmittance of the circulator are automatically computed.
The nominal frequency for the design of the device is chosen as 3 GHz. The circulator
can be expected to perform reasonably well in a narrow frequency band around
3 GHz, and so a frequency range of 2.8 3.2 GHz is studied. It is desired that the
device operates in single mode. Thus a rectangular waveguide cross section of 6.67 cm
by 3.33 cm is selected to set the cut-off frequency for the fundamental TE10 mode to
2.25 GHz. The cut-off frequencies for the two nearest higher modes, the TE20 and
TE01 modes, are both at 4.5 GHz, leaving a reasonable safety margin.
2 |
Model Definition
One of the rectangular ports is excited by the fundamental TE10 mode. At the ports,
the boundaries are transparent to the TE10 mode. The following equation applies to
the electric field vector E inside the circulator:
1
2
j
( r E ) k 0 r --------- E = 0
where r denotes the relative permeability tensor, is the angular frequency, is the
conductivity tensor, 0 is the permittivity of vacuum, r is the relative permittivity
tensor, and k0 is the free space wave number. In this particular model, the conductivity
is zero everywhere. Losses in the ferrite are introduced as complex-valued permittivity
and permeability tensors. The magnetic permeability is of key importance as it is the
anisotropy of this parameter that is responsible for the nonreciprocal behavior of the
circulator. For simplicity, the rather complicated material expressions are predefined in
a text file that is imported into the model. The expressions are also included in the next
section for reference.
THE LOSSY FERRITE MATERIAL MODEL
Complete treatises on the theory of magnetic properties of ferrites can be found in Ref.
1 and Ref. 2. The model assumes that the static magnetic bias field, H0, is much
stronger than the alternating magnetic field of the microwaves, so the quoted
expressions are a linearization for a small-signal analysis around this operating point.
Under these assumptions, and including losses, the anisotropic permeability of a ferrite
magnetized in the positive z direction is given by:
j 0
[ ] = j 0
0 0 0
where
= j 0 xy
= 0 ( 1 + xx )
and the unique elements of the magnetic susceptibility tensor are given by:
3 |
m ( 20 + 2 ( 1 + 2 ) )
0 m ( 20 2 ) + 0 m 2 2
- j ------------------------------------------------------------------------------------ xx = ------------------------------------------------------------------------------------2
2
( 20 2 ( 1 + 2 ) ) + 420 2 2
( 20 2 ( 1 + 2 ) ) + 420 2 2
2 0 m 2
m ( 20 2 ( 1 + 2 ) )
- + j ------------------------------------------------------------------------------------ xy = ------------------------------------------------------------------------------------2
2
( 20 2 ( 1 + 2 ) ) + 4 20 2 2
( 20 2 ( 1 + 2 ) ) + 4 20 2 2
where
0 = 0 H 0
m = 0 M s
0 H
= -----------------2
Here 0 denotes the permeability of free space; is the angular frequency of the
microwave field; 0 is the precession resonance frequency (Larmor frequency) of a
spinning electron in the applied magnetic bias field, H0; m is the electron Larmor
frequency at the saturation magnetization of the ferrite, Ms; and is the gyromagnetic
ratio of the electron. For a lossless ferrite ( = 0), the permeability becomes infinite at
= 0. In a lossy ferrite ( 0), this resonance becomes finite and is broadened. The
loss factor, , is related to the line width, H, of the susceptibility curve near the
resonance as given by the last expression above. The material data,
Ms = 5.41104 A/m, r = 14.5
with an effective loss tangent of 210-4 and H = 3.18103 A/m, are taken for
aluminum garnet from Ref. 2. The applied bias field is set to H0 = 7.96103 A/m. The
electron gyromagnetic ratio taken from Ref. 2 is 1.7591011 C/kg.
4 |
5 |
The plot of the S-parameter from the parametric sweep of sc_ferrite indicates a
minimum for a scale factor of 0.518.
6 |
The plot of the S-parameter from the parametric sweep of sc_chamfer indicates a
minimum for a scale factor of about 3.0.
7 |
At the center frequency most of the standing waves are gone with the optimized values
of sc_ferrite and sc_chamfer.
Figure 5: Electric field norm plot with the optimized sc_ferrite and sc_chamfer values.
8 |
9 |
From the below plot, it should be possible to identify the model at first glance so it has
to display the geometry and some characteristic simulation results.
Reference
1. R.E. Collin, Foundations for Microwave Engineering, 2nd ed., IEEE Press/
Wiley-Interscience, 2000.
2. D.M. Pozar, Microwave Engineering, 3rd ed., John Wiley & Sons Inc, 2004.
Modeling Instructions
From the File menu, choose Open.
Browse to the models Model Library folder and double-click the file
lossy_circulator_3d_geom.mph.
10 |
GEOMETRY 1
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
On the Model toolbar, click Add Material to close the Add Material window.
Air (mat1)
In the Electromagnetic Waves interface, the ferrite is entered as a separate, user-defined
equation model referring to the global variables defined above.
ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN (EMW)
murxy
murxz
muryx
muryy
muryz
murzx
murzy
murzz
11 |
9 Locate the Conduction Current section. From the list, choose User defined. One
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 1 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
5 From the Wave excitation at this port list, choose On.
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 18 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
Port 3
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 19 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
The mesh needs to resolve the local wavelength and, for lossy domains, the skin
depth. The skin depth in the ferrite is large so the main concern is to resolve the
local wavelength. This is done by providing maximum mesh sizes per domain. The
rule of thumb is to use a maximum element size that is one fifth of the local
wavelength (at the maximum frequency) or smaller.
MESH 1
Free Tetrahedral 1
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Free Tetrahedral.
Size 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Free
Tetrahedral 1 and choose Size.
2 In the Settings window for Size, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Domain.
12 |
box.
7 In the associated text field, type 1.5e-2.
Size 2
1 Right-click Free Tetrahedral 1 and choose Size.
2 In the Settings window for Size, locate the Element Size section.
3 Click the Custom button.
4 Locate the Geometric Entity Selection section. From the Geometric entity level list,
choose Domain.
5 Select Domain 2 only.
6 Locate the Element Size Parameters section. Select the Maximum element size check
box.
7 In the associated text field, type 4.5e-3.
8 In the Model Builder window, right-click Mesh 1 and choose Build All.
13 |
STUDY 1
The electric field norm gives a good indication on where the main power is flowing
and where there are standing waves due to reflections from the impedance mismatch
at the center. See Figure 2.
The remaining work is to vary the two design parameters in order to minimize
reflections at the nominal frequency. To do this, perform parametric sweeps over the
design parameters (scale factors). To avoid accumulating a lot of data while solving,
throw away the solution and log only the S-parameter representing reflection in a
table. For this purpose, add a global probe to the model.
DEFINITIONS
1 On the Definitions toolbar, click Probes and choose Global Variable Probe.
2 In the Settings window for Global Variable Probe, click Replace Expression in the
upper-right corner of the Expression section. From the menu, choose Component 1
(comp1)>Electromagnetic Waves, Frequency Domain>Ports>S-parameter,
dB>emw.S11dB - S-parameter, dB, 11 component.
STUDY 1
Modify the study in order to vary the scale factor determining the size of the ferrite
post. The study type is still Frequency Domain.
The parametric sweep over the scale factor is added as an extension to the frequency
domain study.
14 |
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
sc_ferrite
range(0.5,3e-3,0.53)
Parameter unit
5 Locate the Output While Solving section. Select the Accumulated probe table check
box.
6 Find the Memory settings for jobs subsection. From the Keep solutions in memory list,
In Accumulated Probe Table 1, delete column number 2 from the left with heading
freq.
2 Click Table Graph in the window toolbar.
RESULTS
1D Plot Group 3
Compare with the plot shown Figure 3. The plot of the S-parameter indicates a
minimum for a scale factor of 0.518, so freeze the parameter at this value and add a
new study for varying the next scale factor.
1 In the Model Builder window, expand the Results>Tables node.
15 |
DEFINITIONS
Parameters
1 In the Model Builder window, expand the Global>Definitions node, then click
Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
sc_ferrit
e
0.518
0.5180
STUDY 1
Parametric Sweep
1 In the Model Builder window, under Study 1 click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 In the table, enter the following settings:
Parameter name
sc_chamfer
range(2.8,0.04,3.2)
Parameter unit
In Accumulated Probe Table 1, delete column number 2 from the left with heading
freq.
2 Click Table Graph in the window toolbar.
16 |
RESULTS
1D Plot Group 3
See Figure 4. The plot of the S-parameter indicates a minimum for a scale factor of
about 3.0, so leave the parameter at this value and add a study for the frequency
response.
ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN (EMW)
1 In the Model Builder window, expand the Results>Tables node, then click Component
1 (comp1)>Electromagnetic Waves, Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Port Sweep Settings section.
3 Select the Activate port sweep check box.
DEFINITIONS
Parameters
1 In the Model Builder window, under Global>Definitions click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
PortName
1.000
Port name
ADD STUDY
1 On the Study toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Preset Studies>Frequency
Domain.
4 Click Add Study in the window toolbar.
5 On the Study toolbar, click Add Study to close the Add Study window.
STUDY 2
17 |
In the reproduced Figure 5 most of the standing waves are gone at the center
frequency.
Finally plot all the S-parameters as a function of frequency.
S-Parameter (emw)
1 In the Model Builder window, expand the S-Parameter (emw) node, then click Global
1.
2 In the Settings window for Global, locate the y-Axis Data section.
The table should now only contain emw.S11dB, emw.S21dB and emw.S31dB.
6 On the 1D plot group toolbar, click Plot.
18 |
Parametric Sweep
The parametric sweep is used to control which port is excited. It overrides the settings
on individual port features and drives one port at a time using 1 W of input power.
1 In the Model Builder window, under Study 1 click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 In the table, enter the following settings:
Parameter name
PortName
1 2 3
Parameter unit
4 Locate the Output While Solving section. Find the Memory settings for jobs
Derived Values
1 On the Results toolbar, click More Derived Values and choose Other>Global Matrix
Evaluation.
2 In the Settings window for Global Matrix Evaluation, locate the Data section.
19 |
View 3
1 In the Model Builder window, under Component 1 (comp1) right-click Definitions and
choose View.
2 In the Settings window for View, locate the View section.
3 Clear the Show grid check box.
RESULTS
20 |
8 In the Model Builder window, expand the Electric Field (emw) 2 node, then click Slice
1.
9 In the Settings window for Slice, locate the Plane Data section.
10 From the Plane list, choose xy-planes.
11 In the Planes text field, type 1.
12 Right-click Results>Electric Field (emw) 2>Slice 1 and choose Deformation.
13 In the Settings window for Deformation, locate the Expression section.
14 Click emw.Ex,emw.Ey,emw.Ez - Electric field in the upper-right corner of the section.
The port excitation can now be selected on the plot group. For the model
thumbnail, select the second port.
24 In the Model Builder window, click Electric Field (emw) 2.
25 In the Settings window for 3D Plot Group, locate the Data section.
26 From the Parameter value (PortName) list, choose 2.000.
27 On the 3D plot group toolbar, click Plot.
21 |
22 |
1 |
2 |
r(Xg, Yg)
Figure 1: A convex metamaterial lens. Both the shape and the dielectric distribution are
defined on a rectangular domain, and mapped into the deformed state.
1 |
Model Definition
Consider a 2D model geometry as shown in Figure 2. A square air domain, bounded
by a perfectly matched layer (PML) on all sides, encloses a rectangular region in which
the metamaterial lens is defined.
PML
Figure 2: The modeling domain consists of the metamaterial lens in an air domain, and
a surrounding PML. A Gaussian beam is incident from the left.
Model a Gaussian beam entering the domain from the left side, via a surface current
excitation at an interior boundary. The surface current, Js0,can also be thought of as
a displacement current excitation. The waist of the beam is at the boundary, so the
excitation at this boundary can be specified as
y 2
J s0 = exp -------
w 0
where w0 is the waist size. The excitation is at the boundary between a domain of free
space and the PML, and excites a wave propagating in both directionsinto the PML
and into the modeling domain. The wave propagating into the PML is completely
absorbed, and the wave propagating into the domain is diffracted by the lens.
2 |
Both the shape and the dielectric distribution of the metamaterial lens are defined with
respect to the original Cartesian coordinate system, as shown in Figure 1. The true
shape of the lens is described by the relationship
1
--- X g ( 2 Y g 2 )
x = F x ( X g, Y g ) = 2
y
F y ( X g, Y g )
1 2
Y g 1 + --- x
2
where Xg, Yg are the Cartesian coordinates of the undeformed frame.
The dielectric distribution is defined on the original Cartesian domain as:
1 2 2
r = 1 + --- Y g
2
The above expression introduces a variation in the dielectric in the y-coordinate of the
undeformed lens. On the deformed lens, the dielectric varies in both directions.
The Deformed Geometry functionality is used to define the mapping of the dielectric
from the initially rectangular domain onto the desired shape. The deformation and the
dielectric distribution within the lens domain is completely specified via the above
functions.
3 |
Figure 3: The norm of the electric field shows the Gaussian beam diffracted by the
metamaterial lens.
4 |
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
5 |
Expression
Value
Description
f0
3[GHz]
3.000E9 Hz
Operating frequency
lda0
c_const/f0
0.09993 m
w0
lda0*4
0.3997 m
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
GEOMETRY 1
First, create a square for the entire model domain. Add a layer on each side of the
square.
Square 1 (sq1)
1 On the Geometry toolbar, click Primitives and choose Square.
2 In the Settings window for Square, locate the Size section.
3 In the Side length text field, type 3.
4 Locate the Position section. From the Base list, choose Center.
5 Click to expand the Layers section. In the table, enter the following settings:
Layer name
Thickness (m)
Layer 1
lda0
Rectangle 1 (r1)
Add a rectangle for the lens.
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Height text field, type 2.
4 Locate the Position section. From the Base list, choose Center.
6 |
DEFINITIONS
Add a selection for the lens domain which will be recalled frequently while setting up
the model properties.
Explicit 1
1 On the Definitions toolbar, click Explicit.
7 |
Next, add a set of variables for the shape and the dielectric distribution of the
metamaterial lens.
Variables 1a
1 On the Definitions toolbar, click Local Variables.
2 In the Settings window for Variables, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Domain.
4 From the Selection list, choose Lens.
5 Locate the Variables section. In the table, enter the following settings:
8 |
Name
Expression
Unit
Description
xp
0.5[m]*Xg[1/m]*(2-(Yg[1/m])^2)
Mapping of Xg -> x
yp
Yg*(1+(0.5*(xp[1/m])^2))
Mapping of Yg -> y
erp
(1+0.5*(Yg[1/m])^2)^2
Dielectric
distribution
Set up Deformed Geometry. You need to specify Free Deformation, Prescribed Mesh
Displacement and Prescribed Deformation.
1 In the Model Builder window, under Component 1 (comp1) click Deformed Geometry
(dg).
2 In the Settings window for Deformed Geometry, locate the Frame Settings section.
3 From the Geometry shape order list, choose 1.
Free Deformation 1
1 On the Physics toolbar, click Domains and choose Free Deformation.
9 |
10 |
3 In the Settings window for Prescribed Mesh Displacement, locate the Prescribed
Mesh Displacement section.
4 In the dx text field, type xp-Xg.
5 In the dy text field, type yp-Yg.
Prescribed Deformation 1
1 On the Physics toolbar, click Domains and choose Prescribed Deformation.
2 In the Settings window for Prescribed Deformation, locate the Domain Selection
section.
3 From the Selection list, choose Lens.
4 Locate the Prescribed Mesh Displacement section. In the dx text-field array, type
xp-Xg on the first row.
perform the calculation for the out-of-plane component. The in-plane components
are both zero.
section.
3 From the Selection list, choose Lens.
4 Locate the Electric Displacement Field section. From the r list, choose User defined.
default value 1.
11 |
6 Locate the Conduction Current section. From the list, choose User defined. Leave
Surface Current 1
1 On the Physics toolbar, click Boundaries and choose Surface Current.
2 Select Boundary 10 only.
3 In the Settings window for Surface Current, locate the Surface Current section.
4 Specify the Js0 vector as
0
exp(-(y/w0)^2)
MATERIALS
Set all domain with vacuum. The lens domain material properties are explicitly
configured by Wave Equation, Electric 2 in Electromagnetic Waves, Frequency
Domain.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
12 |
1 On the Model toolbar, click Add Material to close the Add Material window.
MESH 1
Size
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
type lda0/10.
5 In the Minimum element size text field, type 0.0012.
6 In the Model Builder window, right-click Mesh 1 and choose Build All. You may zoom
The model is analyzed with two study steps. First, make sure that Stationary study step
is solved only for Deformed Geometry.
Step 1: Stationary
1 In the Model Builder window, under Study 1 click Step 1: Stationary.
2 In the Settings window for Stationary, locate the Physics and Variables Selection
section.
3 In the table, enter the following settings:
Physics interface
Solve for
Discretization
physics
Add a Frequency Domain study step and set as solved only for Electromagnetic Waves,
Frequency Domain.
13 |
2 In the Settings window for Frequency Domain, locate the Study Settings section.
3 In the Frequencies text field, type f0.
4 Locate the Physics and Variables Selection section. In the table, enter the following
settings:
Physics interface
Solve for
Discretization
Deformed Geometry
physics
Add a filled contour plot describing the dielectric distribution over the lens.
2D Plot Group 2
1 On the Model toolbar, click Add Plot Group and choose 2D Plot Group.
2 In the Model Builder window, under Results right-click 2D Plot Group 2 and choose
Contour.
3 In the Settings window for Contour, click Replace Expression in the upper-right
14 |
5 Locate the Coloring and Style section. From the Contour type list, choose Filled.
6 From the Color table list, choose GrayScale.
7 Select the Reverse color table check box. The plot for the dielectric distribution is
shown in Figure 4.
15 |
16 |
Microwave Oven
Introduction
This is a model of the heating process in a microwave oven. The distributed heat source
is computed in a stationary, frequency-domain electromagnetic analysis followed by a
transient heat transfer simulation showing how the heat redistributes in the food.
Model Definition
The microwave oven is a metallic box connected to a 500 W, 2.45 GHz microwave
source via a rectangular waveguide operating in the TE10 mode. Near the bottom of
the oven there is a cylindrical glass plate with a spherical potato placed on top of it. A
part of the potato is cut away for mechanical stability, which also facilitates the creation
of a finite element mesh in the region where it is in contact with the plate. Symmetry
is utilized by simulating only half of the problem. The symmetry cut is applied
vertically through the oven, waveguide, potato, and plate. Figure 1 below shows the
reduced geometry.
1 |
MICROWAVE OVEN
The model uses copper for the walls of the oven and the waveguide. Although resistive
metals losses are expected to be small, the impedance boundary condition on these
walls ensures that they get accounted for. For more information on this boundary
condition, see the section Impedance Boundary Condition in the RF Module Users
Guide. The symmetry cut has mirror symmetry for the electric field and is represented
by the boundary condition n H = 0.
The rectangular port is excited by a transverse electric (TE) wave, which is a wave that
has no electric field component in the direction of propagation. At an excitation
frequency of 2.45 GHz, the TE10 mode is the only propagating mode through the
rectangular waveguide. The cutoff frequencies for the different modes are given
analytically from the relation
c m 2
n 2
( c ) mn = --- ----- + ---
2 a
b
where m and n are the mode numbers and c denotes the speed of light. For the TE10
mode, m = 1 and n = 0. With the dimensions of the rectangular cross section
(a = 7.8 cm and b = 1.8 cm), the TE10 mode is the only propagating mode for
frequencies between 1.92 GHz and 3.84 GHz.
The port condition requires a propagation constant , which at the frequency is given
by the expression
2
2 2
= ------ c
c
With the stipulated excitation at the rectangular port, the following equation is solved
for the electric field vector E inside the waveguide and oven:
1
2
j
( r E ) k 0 r --------- E = 0
where r denotes the relative permeability, j the imaginary unit, the conductivity,
the angular frequency, r the relative permittivity, and 0 the permittivity of free space.
The model uses material parameters for air: = 0 and r = r = 1. In the potato the
same parameters are used except for the permittivity which is set to r = 65 20j where
the imaginary part accounts for dielectric losses. The glass plate has = 0, r = 1 and
r = 2.55.
2 |
MICROWAVE OVEN
3 |
MICROWAVE OVEN
Figure 3: Temperature in the center of the potato during the first 5 seconds of heating.
4 |
MICROWAVE OVEN
Modeling Instructions
From the File menu, choose New.
NEW
Add a Frequency Domain study type for the Electromagnetic Waves, Frequency
Domain interface.
5 In the Select study tree, select Custom Studies>Preset Studies for Some Physics
Interfaces>Frequency Domain.
6 Click Done.
STUDY 1
Add a Time Dependent study for the Heat Transfer in Solids interface.
The Frequency Domain study is only used for the Electromagnetics physics
interface, whereas the time-dependent study is only applicable for the heat transfer
physics interface. Notice that the electromagnetic heat source will be computed
first, and then used in the time-dependent heat transfer study step.
2 In the Settings window for Time Dependent, locate the Physics and Variables
Selection section.
5 |
MICROWAVE OVEN
Solve for
Discretization
physics
Solve for
Discretization
physics
DEFINITIONS
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
6 |
Name
Expression
Value
Description
wo
267[mm]
0.2670 m
Oven width
do
270[mm]
0.2700 m
Oven depth
ho
188[mm]
0.1880 m
Oven height
wg
50[mm]
0.05000 m
Waveguide width
dg
78[mm]
0.07800 m
Waveguide depth
hg
18[mm]
0.01800 m
Waveguide height
rp
113.5[mm]
0.1135 m
hp
6[mm]
0.006000 m
bp
15[mm]
0.01500 m
rpot
31.5[mm]
0.03150 m
Potato radius
T0
8[degC]
281.2 K
Initial potato
temperature
MICROWAVE OVEN
GEOMETRY 1
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type wo.
4 In the Depth text field, type do/2.
5 In the Height text field, type ho.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type wg.
4 In the Depth text field, type dg/2.
5 In the Height text field, type hg.
6 Locate the Position section. In the x text field, type wo.
7 In the z text field, type ho-hg.
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type rp.
4 In the Height text field, type hp.
5 Locate the Position section. In the x text field, type wo/2.
6 In the z text field, type bp.
Sphere 1 (sph1)
1 On the Geometry toolbar, click Sphere.
2 In the Settings window for Sphere, locate the Size section.
3 In the Radius text field, type rpot.
4 Locate the Position section. In the x text field, type wo/2.
5 In the z text field, type rpot+bp.
6 Click the Build All Objects button.
The sphere you have created for the potato now overlaps the glass plate. This in itself
is not a problem, but where the sphere touches the bottom of the glass plate, you
7 |
MICROWAVE OVEN
risk getting very thin mesh elements. To avoid this problem, you will delete the part
of the sphere that overlaps the cylinder. To retain the cylinder after this operation,
begin by making a copy of it.
Copy 1 (copy1)
1 On the Geometry toolbar, click Transforms and choose Copy.
2 Select the object cyl1 only.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object sph1 only to add it to the Objects to add list.
Finally, make a geometric operation to keep only the part of the potato and the plate
that overlaps the half oven.
Compose 1 (co1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Compose.
2 Click in the Graphics window and then press Ctrl+A to select all objects.
3 In the Settings window for Compose, locate the Compose section.
4 In the Set formula text field, type (blk1+blk2)*(dif1+copy1).
5 Select the Keep input objects check box.
6 Click the Build All Objects button.
section.
3 From the Geometric entity level list, choose Object.
4 Select the objects dif1 and copy1 only.
5 Click the Build All Objects button.
6 Click the Wireframe Rendering button on the Graphics toolbar.
8 |
MICROWAVE OVEN
DEFINITIONS
Create the following selections definitions in order to make Domain and Boundary
selections easier as you walk through these model instructions. Note that if you have
problems finding certain numbers, you can always choose Selection List.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 1 and choose Rename.
3 In the Rename Explicit dialog box, type Potato in the New label text field.
4 Click OK.
5 Select Domain 3 only.
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 2 and choose Rename.
3 In the Rename Explicit dialog box, type Plate in the New label text field.
4 Click OK.
5 Select Domain 2 only.
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 3 and choose Rename.
3 In the Rename Explicit dialog box, type Air in the New label text field.
4 Click OK.
5 Select Domains 1 and 4 only.
Explicit 4
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 4 and choose Rename.
3 In the Rename Explicit dialog box, type No Heat Transfer in the New label text
field.
4 Click OK.
9 |
MICROWAVE OVEN
Explicit 5
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 5 and choose Rename.
3 In the Rename Explicit dialog box, type Port Boundary in the New label text field.
4 Click OK.
5 In the Settings window for Explicit, locate the Input Entities section.
6 From the Geometric entity level list, choose Boundary.
7 Select Boundary 23 only.
Explicit 6
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 6 and choose Rename.
3 In the Rename Explicit dialog box, type Symmetry Boundaries in the New label text
field.
4 Click OK.
5 In the Settings window for Explicit, locate the Input Entities section.
6 From the Geometric entity level list, choose Boundary.
7 Select Boundaries 2, 7, 10, and 19 only.
Explicit 7
1 On the Definitions toolbar, click Explicit.
2 In the Model Builder window, under Component 1 (comp1)>Definitions right-click
Explicit 7 and choose Rename.
3 In the Rename Explicit dialog box, type Metal Boundaries in the New label text
field.
4 Click OK.
5 In the Settings window for Explicit, locate the Input Entities section.
6 From the Geometric entity level list, choose Boundary.
7 Select Boundaries 1, 35, 17, and 2022 only.
MATERIALS
Next, define the materials. Air and Copper are already in the Material Library.
10 |
MICROWAVE OVEN
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
Air (mat1)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Air (mat1).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Air.
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Right-click Material 2 (mat2) and choose Rename.
3 In the Rename Material dialog box, type Potato in the New label text field.
4 Click OK.
5 In the Settings window for Material, locate the Geometric Entity Selection section.
6 From the Selection list, choose Potato.
7 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
65-20*
j
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Thermal conductivity
0.55
W/
(mK)
Basic
Density
rho
1050
kg/m
Basic
Cp
3.64e3
J/
(kgK)
Basic
Material 3 (mat3)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Right-click Material 3 (mat3) and choose Rename.
3 In the Rename Material dialog box, type Glass in the New label text field.
11 |
MICROWAVE OVEN
4 Click OK.
5 In the Settings window for Material, locate the Geometric Entity Selection section.
6 From the Selection list, choose Plate.
7 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
2.55
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
You do not need to define the listed thermal properties, as the glass plate will not
be in the thermal part of the model.
ADD MATERIAL
1 On the Model toolbar, click Add Material to close the Add Material window.
Copper (mat4)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Copper
(mat4).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 From the Selection list, choose Metal Boundaries.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
12 |
MICROWAVE OVEN
5 Locate the Analysis Methodology section. From the Methodology options list, choose
Fast.
For the electromagnetic part of the problem, begin by defining the input port.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 In the Settings window for Port, locate the Boundary Selection section.
3 From the Selection list, choose Port Boundary.
4 Locate the Port Properties section. From the Wave excitation at this port list, choose
On.
5 In the Pin text field, type 500.
6 Locate the Port Mode Settings section. Specify the E0 vector as
0
cos(pi*y/dg)[V/m]
13 |
MICROWAVE OVEN
In order to solve for the temperature in the potato only, deselect all the other domains.
H E A T TR A N S F E R I N S O L I D S ( H T )
On the Physics toolbar, click Electromagnetic Waves, Frequency Domain (emw) and
choose Heat Transfer in Solids (ht).
1 In the Model Builder window, under Component 1 (comp1) click Heat Transfer in Solids
(ht).
2 Select Domain 3 only.
Initial Values 1
1 In the Model Builder window, under Component 1 (comp1)>Heat Transfer in Solids (ht)
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
Solution 1
1 On the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 1 node.
14 |
MICROWAVE OVEN
Temperature (ht)
The Graphics window shows the temperature distribution on the surface of the potato
after 5 s. Change the unit to degC to reproduce Figure 4.
1 In the Model Builder window, expand the Results>Temperature (ht) node, then click
Surface 1.
2 In the Settings window for Surface, locate the Expression section.
3 From the Unit list, choose degC.
4 On the 3D plot group toolbar, click Plot.
Modify an existing plot group to plot the resistive heating on the symmetry plane.
text field.
6 Click OK.
Resistive Heating
1 In the Model Builder window, expand the Results>Resistive Heating node.
2 Right-click Isosurface 1 and choose Delete.
3 Right-click Resistive Heating and choose Slice.
4 In the Settings window for Slice, click Replace Expression in the upper-right corner
15 |
MICROWAVE OVEN
The plot should now look like Figure 2. Next, add a nice visualization of the
electromagnetic fields to the temperature plot.
Temperature (ht)
1 In the Model Builder window, under Results right-click Temperature (ht) and choose
Slice.
2 In the Settings window for Slice, click Replace Expression in the upper-right corner
Derived Values
Make a volume integral of the microwave heating to find out how much of the energy
is absorbed in the potato.
1 On the Results toolbar, click More Derived Values and choose Integration>Volume
Integration.
2 In the Settings window for Volume Integration, locate the Data section.
3 From the Data set list, choose Study 1/Solution Store 1.
4 Click Replace Expression in the upper-right corner of the Expression section. From
16 |
MICROWAVE OVEN
TA BL E
Data Sets
1 On the Results toolbar, click Cut Point 3D.
2 In the Settings window for Cut Point 3D, locate the Point Data section.
3 In the x text field, type 0.134.
4 In the y text field, type 0.
5 In the z text field, type 0.047.
1D Plot Group 4
1 On the Results toolbar, click 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Cut Point 3D 1.
4 Click to expand the Title section. From the Title type list, choose Manual.
5 In the Title text area, type Temperature in potato.
6 Locate the Plot Settings section. Select the x-axis label check box.
7 In the associated text field, type Time (s).
8 On the 1D plot group toolbar, click Point Graph.
9 In the Settings window for Point Graph, click Replace Expression in the upper-right
corner of the y-axis data section. From the menu, choose Component 1>Heat Transfer
in Solids>Temperature>T - Temperature.
10 Locate the y-Axis Data section. From the Unit list, choose degC.
11 On the 1D plot group toolbar, click Plot.
17 |
MICROWAVE OVEN
18 |
MICROWAVE OVEN
M i c r o w a v e Fi lt er on PCB
This model analyzes the transfer function of a low-pass filter on a printed circuit board.
Introduction
Microstrip filters can be fabricated directly on a printed circuit board (PCB) with a
microstrip line going from the input to the output. Along the microstrip line there are
a number of stubs of certain lengths and widths. The design of the filter involves
choosing the impedance of the microstrip line, the impedance of the stub microstrips,
and the stub lengths. This particular filter is based on a textbook example from Ref. 2.
It is also used as example in Ref. 1, which contains results from other simulation tools
and methods and is freely available online. The filter has a seven-pole low-pass
Chebyshev response with a cutoff frequency of 1 GHz. The entire layout including the
dielectric layer is imported from an ODB++(X) file using the ECAD Import feature.
Model Definition
The model uses the Electromagnetic Waves interface that solves the vector Helmholtz
wave equation. The PCB layout of the Chebyshev filter is imported using the ECAD
Import feature, which creates the entire 3D geometry from the layout and stackup
information of an ODB++ file.
The cutoff frequency of the filter is 1 GHz by design, and the dielectric layer of the
PCB has a relative permittivity of 10.8. The metal layers are modeled as perfect electric
conductors with zero thickness, thereby avoiding a dense meshing of thin conductive
layers. The ECAD import has an option that causes the import engine to ignore the
thickness of metal layers and insert them as faces on the dielectric layer. The width of
the microstrip line is 0.1 mm and the width of the stubs is 5 mm.
The characteristics of the filter are sensitive to the placement and length of the stubs;
therefore this model also analyzes the change in filter characteristics as a function of
mechanical deformation. This is done by adding Solid Mechanics and Moving Mesh
interfaces. The Moving Mesh interface is necessary to enable the Electromagnetic
Waves interface to account for the deformation of the PCB. The deformation comes
from a uniform load across the board with fixed input and output faces.
Because the filter cutoff should be close to 1 GHz, the frequency is swept from
750 MHz up to 1.5 GHz. The first solution step performs this sweep for the
1 |
Figure 1: The graph shows the total displacement of the PCB, due to the load.
Although the PCB deformation is fairly large the S-parameter curve does not change
that much. The cutoff frequency is shifted less than 10 MHz when the force is applied.
Figure 2 displays the difference between the S-parameter curves with and without an
2 |
applied force.
Figure 2: A comparison between the S11 and S21 parameters before and after a force of
40 N has been applied on the PCB. The red and turquoise lines correspond to the
S-parameter curves for the filter with an applied force.
References
1. D.V. Tosic and M. Potrebic, Software Tools for Research and Education,
Microwave Review, vol. 12, no. 2, p. 45, 2006.
2. J.-S.G. Hong and M.J. Lancaster, Microstrip Filters for RF/Microwave
Applications, John Wiley & Sons, 2001.
Modeling Instructions
From the File menu, choose New.
3 |
NEW
The following steps define the parameters for the frequency sweep.
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
fstart
750[MHz]
7.500E8 Hz
Start frequency
fstop
1.5[GHz]
1.500E9 Hz
Stop frequency
fstep
50[MHz]
5.000E7 Hz
Frequency step
GEOMETRY 1
The major part of the geometry is imported using the ECAD Import tool, which
automatically reads a PCB layout and extrudes the layers to a 3D geometry.
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
4 |
Import 1 (imp1)
1 On the Geometry toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 From the Geometry import list, choose ECAD file (ODB++).
4 Click Browse.
5 Browse to the models Model Library folder and double-click the file
pcb_microwave_filter_with_stress.xml.
It is located in the same folder as specified in the Model Library Path on page 3.
This path originates from the folder models under the COMSOL Multiphysics
installation folder.
6 From the Type of import list, choose Metal shell.
This will import all metal layers as faces. which drastically reduce the problem size
in this model without compromising the accuracy.
7 Click Import.
The PCB is now placed within the drawn block. To see the PCB, select wireframe
rendering.
8 Click the Wireframe Rendering button on the Graphics toolbar.
5 |
The next step is to add boundary faces for the input and output ports.
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
The size and position of the rectangle, as defined below, will perfectly fit it to the
microstrip line.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 1.125.
4 In the Height text field, type 1.27.
5 Locate the Position section. In the xw text field, type -4.77.
6 In the yw text field, type -0.635.
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
6 |
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
ADD MATERIAL
The relative permittivity is modified to agree with the value used in Ref. 1. The FR4
material is selected to provide parameters for the solid mechanics simulation.
3 In the Settings window for Material, locate the Material Contents section.
4 In the table, enter the following settings:
Property
Relative permittivity
Name
Value
Unit
Property group
epsilonr
10.8
Basic
5 On the Model toolbar, click Add Material to close the Add Material window.
7 |
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 10 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Wave excitation at this port list, choose On.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 16 only.
These boundaries represent the microstrip line and the ground plane on the PCB.
MESH 1
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
Size 1
1 Right-click Component 1 (comp1)>Mesh 1>Free Triangular 1 and choose Size.
2 In the Settings window for Size, locate the Element Size section.
3 Click the Custom button.
4 Locate the Element Size Parameters section. Select the Maximum element size check
box.
5 In the associated text field, type 1.
8 |
Size
1 In the Model Builder window, right-click Mesh 1 and choose Free Tetrahedral.
2 In the Settings window for Size, locate the Element Size section.
3 From the Predefined list, choose Coarse.
4 Click the Custom button.
5 Locate the Element Size Parameters section. In the Minimum element size text field,
type 0.1.
Allow the Minimum element size to be as small as 0.1 mm to generate a fine mesh
along the thin microstrip.
6 Click the Build All button.
DEFINITIONS
Probe plotting is a convenient technique to plot while solving, which is very useful for
parameter sweeps. It is possible to discover problems before the solution step has
finished, and then stop the sweep to save time. It is also useful in situations when the
solver does more steps than it stores in the output. The probe plot will contain all steps
that the solver takes.
1 In the Model Builder window, under Component 1 (comp1) right-click Definitions and
field.
3 Locate the Expression section. Click emw.S11dB - S-parameter, dB, 11 component in
the upper-right corner of the section. In the Model Builder window, right-click
Definitions and choose Global Variable Probe.
4 In the Settings window for Global Variable Probe, type S21 in the Variable name text
field.
5 Locate the Expression section. Click emw.S21dB - S-parameter, dB, 21 component in
9 |
The following instructions adds physics from the Solid Mechanics and the Moving
Mesh interfaces for the simulations of the deformed PCB.
ADD PHYSICS
1 On the Model toolbar, click Add Physics to open the Add Physics window.
2 Go to the Add Physics window.
3 In the Add physics tree, select Structural Mechanics>Solid Mechanics (solid).
4 Click Add to Component in the window toolbar.
ADD PHYSICS
10 |
1 In the Model Builder window, under Component 1 (comp1) click Solid Mechanics
(solid).
2 Select Domain 2 only.
DEFINITIONS
Parameters
1 In the Model Builder window, under Global>Definitions click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
fload
40[N]
40.00 N
Load on PCB
D E F I N I T I O N S O F P C B VO L U M E
The following steps describe how to measure the volume of the PCB and then copy
and paste the value in a parameter definition
GEOMETRY 1
Import 1 (imp1)
1 Click the Select Domains button on the Graphics toolbar.
2 In the Model Builder window, under Component 1 (comp1)>Geometry 1 click Import
1 (imp1).
3 On the object imp1, select Domain 1 only.
4 On the Geometry toolbar, click Measure.
Click the Measure button from the toolbar. The volume of the PCB domain is
displayed in the Messages window.
Copy the volume of the PCB from the Messages table.
DEFINITIONS
Parameters
1 In the Model Builder window, under Global>Definitions click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
11 |
3 In the table, enter the following settings (by pasting in the previously copied
volume):
Name
Expression
Value
Description
3357.0[mm^3]
3.357E-6 m
Volume of PCB
On the Physics toolbar, click Solid Mechanics (solid) and choose Moving Mesh (ale).
Prescribed Deformation 1
1 On the Physics toolbar, click Domains and choose Prescribed Deformation.
2 Select Domain 2 only.
3 In the Settings window for Prescribed Deformation, locate the Prescribed Mesh
Displacement section.
4 In the dx text-field array, type u on the first row.
5 In the dy text-field array, type v on the 2nd row.
6 In the dz text-field array, type w on the 3rd row.
Free Deformation 1
1 On the Physics toolbar, click Domains and choose Free Deformation.
2 Select Domain 1 only.
Body Load 1
1 On the Physics toolbar, click Domains and choose Body Load.
2 Select Domain 2 only.
12 |
3 In the Settings window for Body Load, locate the Force section.
4 Specify the FV vector as
0
-fload/V
Fixed Constraint 1
1 On the Physics toolbar, click Boundaries and choose Fixed Constraint.
2 Select Boundaries 6, 10, 12, and 1517 only.
STUDY 1
Solve for
Discretization
Solid Mechanics
physics
Moving Mesh
physics
ADD STUDY
1 On the Model toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Custom Studies>Preset
Studies for Some Physics Interfaces>Stationary.
4 Find the Physics interfaces in study subsection. In the table, enter the following
settings:
Physics
Solve
13 |
STUDY 2
Solve for
Discretization
Solid Mechanics
physics
Moving Mesh
physics
section.
10 Select the Plot check box.
11 From the Plot group list, choose S-Parameter (emw).
12 On the Study toolbar, click Compute.
TABLE
S-Parameter (emw) 1
To compare the S-parameters for the initial and the stressed PCB, add the S-parameter
from the first different solutions.
1 On the 1D plot group toolbar, click Global.
14 |
Stress (solid)
1 In the Model Builder window, expand the Results>Stress (solid) node.
You may expand all subnodes under Results>Stress (solid) to reach Deformation.
2 In the Model Builder window, expand the Results>Stress (solid)>Surface 1 node.
3 Right-click Deformation and choose Disable.
4 In the Model Builder window, under Results>Stress (solid) click Surface 1.
5 In the Settings window for Surface, click Replace Expression in the upper-right corner
15 |
16 |
P la s m o ni c W i r e G rat i n g
Introduction
A plane electromagnetic wave is incident on a wire grating on a dielectric substrate.
The model computes transmission and reflection coefficients for the refraction,
specular reflection, and first order diffraction.
Model Definition
Figure 1 shows the considered grating, with a gold wire on a dielectric material with
refractive index n. The grating constant, or the distance between the wires, is d. A
plane-polarized wave traveling through a medium with refractive index n is incident
on the grating, at an angle in a plane perpendicular to the grating.
n
Au
Figure 1: The modeled grating. The model considers a unit cell of a slice through this
geometry. The grating is assumed to consist of an infinite number of infinitely long wires.
If the wavelengths involved in the model are sufficiently short compared to the grating
constant, one or several diffraction orders can be present. The diagram in Figure 2
1 |
shows two transmissive paths taken by light incident on adjacent cells of the grating,
exactly one grating constant apart.
d sin()
n
n
m
d sin(m)
Figure 2: The geometric path lengths of two transmitted parallel beams. The optical path
length is the geometric path length multiplied by the local refractive index.
The criterion for positive interference is that the difference in optical path length along
the two paths equals an integer number of vacuum wavelengths, or:
m 0 = d ( n sin m n sin )
(1)
(2)
Note that for the special cases of perpendicular and grazing incidence, the right-hand
side of the inequality evaluates to dn and d(n + n) respectively.
2 |
d sin(m)
d sin()
m
n
n
(3)
(4)
The model uses n = 1 for air and n = 1.2 for the dielectric substrate. Allowing for
arbitrary angles of incidence and with a grating constant d = 400 nm, Equation 2 sets
the validity limit to vacuum wavelengths greater than 440 nm. The model uses
0 = 441 nm. For the wire, a complex-valued permittivity of 1.75 5.4i
approximates that of gold at the corresponding frequency.
The performance of the grating depends on the polarization of the incident wave.
Therefore both a transverse electric (TE) and a transverse magnetic (TM) case are
considered. The TE wave has the electric field component in the z direction, out of the
modeling xy-plane. For the TM wave, the electric field vector is pointing in the
xy-plane and perpendicular to the direction of propagation, whereas the magnetic field
3 |
has only a component in the z direction. The angle of incidence is for both cases swept
from 0 to /2, with a pitch of /40.
4 |
5 |
6 |
MODEL
DESCRIPTION
na
nb
alpha
Angle of incidence
Internal
-1
Internal
beta
Refraction angle
7 |
MODEL
DESCRIPTION
Internal
-1
Internal
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
8 |
Expression
Value
Description
na
1.0000
nb
1.2
1.2000
Refractive index,
dielectric
400[nm]
4.0000E-7 m
Grating constant
lam0
441[nm]
4.4100E-7 m
Vacuum wavelength
f0
c_const/lam0
6.7980E14 1/s
Frequency
alpha
0.0000
Angle of incidence
beta
asin(na*sin(alpha)/nb)
0.0000 rad
Refraction angle
Create the geometry entirely in terms of the grating constant, for easy scalability.
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type d.
4 In the Height text field, type 3*d.
5 Click the Build Selected button.
Rectangle 2 (r2)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type d.
4 In the Height text field, type 3*d.
5 Locate the Position section. In the y text field, type -3*d.
6 Click the Build Selected button.
Circle 1 (c1)
1 On the Geometry toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type d/5.
9 |
The geometry now consists of two rectangular domains for the air and the dielectric,
and a circle centered on their intersection. You can remove the line through the circle
if you first create a union of the objects.
Union 1 (uni1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Union.
2 Click the Build Selected button.
10 |
Before setting up the materials, define which constitutive relations you want to use in
the Electromagnetic Waves interface.
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
Name
Value
Unit
Property group
Refractive index
na
Refractive index
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Select Domain 1 only.
3 In the Settings window for Material, locate the Material Contents section.
11 |
Name
Value
Unit
Property group
Refractive index
nb
Refractive index
Material 3 (mat3)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Select Domain 3 only.
3 In the Settings window for Material, locate the Material Contents section.
4 In the table, enter the following settings:
Property
Name
Value
Unit
Property group
epsilonBis
5.4
Dielectric losses
epsilonPrim
-1.75
Dielectric losses
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
In the first version of this model, you will assume a TE-polarized wave. This means that
Ex and Ey will be zero throughout the geometry, and that you consequently only need
to solve for Ez.
1 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Components section.
2 From the Electric field components solved for list, choose Out-of-plane vector.
Now define the excitation port. A periodic port assumes that the structure is periodic
and simplifies the setup of ports for the diffraction orders.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
12 |
Notice that you define the electric field by only setting the amplitude. A phase factor
should not be entered.
6 Locate the Port Mode Settings section. Specify the E0 vector as
0
The order in which you set up the ports will determine how the S-parameters are
labeled. You have just created Port 1 for the excitation. If you set up the next port for
the transmission of the purely refracted beam, the S21-parameter will contain
information on the zero order transmission.
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 2 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Periodic.
5 Locate the Port Mode Settings section. Specify the E0 vector as
0
z
The angle of incidence on the exit side corresponds to the angle of incidence an
incident wave on that side would have to provide the correct propagation angle in
the material. Notice that this also means that the sign is opposite that on the entry
side.
13 |
Port 1
Continue with the ports for the reflected diffraction orders. Since these diffraction
orders are not propagating at normal incidence, you have to add the ports manually.
Diffraction Order 1
1 On the Physics toolbar, click Attributes and choose Diffraction Order.
2 In the Settings window for Diffraction Order, locate the Port Mode Settings section.
3 From the Components list, choose Out-of-plane vector.
4 In the m text field, type -1.
Diffraction Order 2
1 Right-click Component 1 (comp1)>Electromagnetic Waves, Frequency Domain
(emw)>Port 1>Diffraction Order 1 and choose Duplicate.
2 In the Settings window for Diffraction Order, locate the Port Properties section.
3 In the Port name text field, type 4.
4 Locate the Port Mode Settings section. In the m text field, type 1.
The transmitted diffraction orders are propagating at normal incidence. Thus, you can
create them automatically by clicking the Compute Diffraction Orders button.
Port 2
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw) click Port 2.
2 In the Settings window for Port, locate the Port Mode Settings section.
3 Click Compute Diffraction Orders. You now find the Diffraction Order ports as
subfeatures to Port 2.
Periodic Condition 1
1 On the Physics toolbar, click Boundaries and choose Periodic Condition.
2 Select Boundaries 1, 3, 7, and 8 only.
3 In the Settings window for Periodic Condition, locate the Periodicity Settings
section.
4 From the Type of periodicity list, choose Floquet periodicity.
The wave vector in the direction for the periodicity is used by the periodic port.
Thus, you can use that wave vector also for the Floquet periodic condition.
14 |
5 From the k-vector for Floquet periodicity list, choose From periodic port.
MESH 1
The periodic boundary conditions perform better if the mesh is identical on the
periodicity boundaries. This is especially important when dealing with vector degrees
of freedom, as will be the case in the TM version of this model.
1 In the Model Builder window, under Component 1 (comp1) click Mesh 1.
2 In the Settings window for Mesh, locate the Mesh Settings section.
3 From the Sequence type list, choose User-controlled mesh.
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Free
Triangular 1 and choose Delete.
Size
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Size.
2 In the Settings window for Size, locate the Element Size section.
3 From the Predefined list, choose Extra fine.
Edge 1
1 In the Model Builder window, right-click Mesh 1 and choose More Operations>Edge.
2 Select Boundaries 1 and 3 only.
Copy Edge 1
1 Right-click Mesh 1 and choose More Operations>Copy Edge.
2 Select Boundary 3 only.
3 In the Settings window for Copy Edge, locate the Destination Boundaries section.
4 Select the Active toggle button.
5 Select Boundary 8 only.
Copy Edge 2
1 Right-click Mesh 1 and choose More Operations>Copy Edge.
2 Select Boundary 1 only.
3 In the Settings window for Copy Edge, locate the Destination Boundaries section.
4 Select the Active toggle button.
5 Select Boundary 7 only.
15 |
Free Triangular 1
1 Right-click Mesh 1 and choose Free Triangular.
2 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Free
Triangular 1 and choose Build All.
To set up the study to sweep for the angle of incidence, some modifications of the
solver is required.
STUDY 1
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
Parameter unit
alpha
5 Click Range.
6 In the Range dialog box, type 0 in the Start text field.
7 In the Stop text field, type pi/2-pi/40.
8 In the Step text field, type pi/40.
9 Click Replace.
10 On the Study toolbar, click Compute.
RESULTS
16 |
like Figure 4.
Rename the plot group to make it clear that it shows the TE solution.
5 Right-click Results>Electric Field (emw) and choose Rename.
6 In the Rename 2D Plot Group dialog box, type 2D Plot Group TE in the New label
text field.
7 Click OK.
Add a 1D plot to look at the various orders of reflection and transmission versus the
angle of incidence.
1D Plot Group 2
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, click to expand the Title section.
3 From the Title type list, choose Manual.
4 In the Title text area, type Reflection and Transmission of TE Wave.
5 Locate the Plot Settings section. Select the x-axis label check box.
6 In the associated text field, type Angle of incidence (degrees).
7 Select the y-axis label check box.
8 In the associated text field, type Reflection and transmission coefficients.
9 On the 1D plot group toolbar, click Global.
Since the diffraction orders are not propagating for all angles, use a logic expression
to set the S-parameter to zero when the wave is evanescent and to multiply with one
when it is propagating.
10 In the Settings window for Global, locate the y-Axis Data section.
11 In the table, enter the following settings:
Expression
Unit
abs(emw.S11)^2
abs(emw.S21)^2
abs(emw.S31)^2*(imag(emw.beta_3)==0)
abs(emw.S41)^2*(imag(emw.beta_4)==0)
17 |
Description
Expression
Unit
abs(emw.S51)^2*(imag(emw.beta_5)==0)
abs(emw.S61)^2*(imag(emw.beta_6)==0)
Description
12 Locate the x-Axis Data section. From the Parameter list, choose Expression.
13 In the Expression text field, type alpha*180/pi.
14 Click to expand the Coloring and style section. Locate the Coloring and Style section.
Find the Line markers subsection. From the Marker list, choose Cycle.
15 Click to expand the Legends section. From the Legends list, choose Manual.
16 In the table, enter the following settings:
Legends
R<sub>0</sub>
T<sub>0</sub>
R<sub>-1</sub>
R<sub>1</sub>
T<sub>-1</sub>
T<sub>1</sub>
text field.
20 Click OK. The plot should now look like Figure 6.
The remaining instructions show to alter the physics so that you solve for an incident
TM wave.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
1 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Components section.
2 From the Electric field components solved for list, choose In-plane vector. You will
Port 1
The easiest way to specify a TM wave is to define the magnetic field, since only the z
component is used.
18 |
Diffraction Order 1
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw)>Port 1 click Diffraction Order 1.
2 In the Settings window for Diffraction Order, locate the Port Mode Settings section.
3 From the Components list, choose In-plane vector.
Diffraction Order 2
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw)>Port 1 click Diffraction Order 2.
2 In the Settings window for Diffraction Order, locate the Port Mode Settings section.
3 From the Components list, choose In-plane vector.
Port 2
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw) click Port 2.
2 In the Settings window for Port, locate the Port Mode Settings section.
3 From the Input quantity list, choose Magnetic field.
4 Specify the H0 vector as
0
5 Click Compute Diffraction Orders to change components for the diffraction orders
19 |
ADD STUDY
1 On the Model toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Preset Studies>Frequency
Domain.
4 Click Add Study in the window toolbar.
5 On the Model toolbar, click Add Study to close the Add Study window.
STUDY 2
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
Parameter unit
alpha
5 Click Range.
6 In the Range dialog box, type 0 in the Start text field.
7 In the Stop text field, type pi/2-pi/40.
8 In the Step text field, type pi/40.
9 Click Replace.
10 On the Study toolbar, click Compute.
RESULTS
20 |
text field.
7 Click OK. You have now reproduced Figure 5.
For the transmission and the reflection of the TM waves, copy and reuse the 1D plot
for the TE waves.
1D Plot Group TE 1
1 In the Model Builder window, under Results right-click 1D Plot Group TE and choose
Duplicate.
2 In the Settings window for 1D Plot Group, locate the Title section.
3 In the Title text area, type Reflection and Transmission of TM Wave.
4 Locate the Data section. From the Data set list, choose Study 2/Solution 2.
5 On the 1D plot group toolbar, click Plot.
6 Right-click Results>1D Plot Group TE 1 and choose Rename.
7 In the Rename 1D Plot Group dialog box, type 1D Plot Group TM in the New label
text field.
8 Click OK. Compare the resulting plot with that in Figure 7.
21 |
22 |
Quarter-Wave Transformer
Introduction
Transmission lines are used when the frequency of the electromagnetic signals is so
high that the wave nature of the signals must be taken into account. A consequence of
the wave nature is that the signals are reflected if there are abrupt changes of the
characteristic impedance along the transmission line. Similarly, the load impedance,
ZL, at the end of the transmission line must match its characteristic impedance, Z0.
Otherwise there will be reflections from the transmission lines end.
A quarter-wave transformer (see Figure 1) is a component that can be inserted
between the transmission line and the load to match the load impedance to the
transmission lines characteristic impedance. To get this functionality, the transformer
must be a quarter of a wavelength long and the relation between the impedances
involved must be
Z in
Z------- = -----Z
ZL
(1)
If the length and the impedance requirements are fulfilled, the load impedance will not
give rise to any reflections.
Typically, the characteristic impedance of transmission lines, Z0, is 50 . Thus, Zin in
Equation 1 should be set to
Z in = Z 0 = 50
(2)
1 |
Q U A R T E R - W A V E TR A N S F O R M E R
Model Definition
The 1D geometry of the model consists of two line intervals. Each line interval
represents a separate transmission line, with different electrical parameters (distributed
capacitance and inductance) and lengths.
To excite and terminate the transmission lines, use lumped ports. This also makes it
easy to obtain the reflection (S11) and transmission (S21) coefficients for the system.
2 |
Q U A R T E R - W A V E TR A N S F O R M E R
Figure 2: Absolute value of the voltage versus the x-coordinate. The quarter-wave
transformer starts at x-coordinate 0.02 m.
3 |
Q U A R T E R - W A V E TR A N S F O R M E R
Figure 3: Spectral response for the transmission line. Notice that the transmission
coefficient (S21) peaks at the frequency (1 GHz) for which the transformer is a
quarter-wave long. At that frequency the reflection coefficient (S11) is zero (approaches
negative infinity which the dB scale used in the graph).
Modeling Instructions
From the File menu, choose New.
NEW
4 |
Q U A R T E R - W A V E TR A N S F O R M E R
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
DEFINITIONS
First add some parameters that defines the electrical and geometrical properties of the
transmission lines.
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
L1
2.5e-6[H/m]
2.5000E-6 H/m
Distributed inductance,
first transmission line
C1
1e-9[F/m]
1.0000E-9 F/m
Distributed capacitance,
first transmission line
1[GHz]
1.0000E9 Hz
Frequency
wl1
1/(f*sqrt(L1*C1))
0.020000 m
Wavelength, first
transmission line
d1
wl1
0.020000 m
Z1
sqrt(L1/C1)
50.000
Characteristic impedance,
first transmission line
ZL
4*Z1
200.00
Terminating impedance
Z2
sqrt(Z1*ZL)
100.00
Characteristic impedance,
second transmission line
C2
C1
1.0000E-9 F/m
Distributed capacitance,
second transmission line
L2
C2*Z2^2
1.0000E-5 H/m
Distributed inductance,
second transmission line
wl2
1/(f*sqrt(L2*C2))
0.010000 m
Wavelength, second
transmission line
d2
wl2/4
0.0025000 m
Length, second
transmission line
hmax
d2/10
2.5000E-4 m
Maximum discretization
step
5 |
Q U A R T E R - W A V E TR A N S F O R M E R
GEOMETRY 1
Interval 1 (i1)
1 On the Geometry toolbar, click Interval.
2 In the Settings window for Interval, locate the Interval section.
3 From the Number of intervals list, choose Many.
4 In the Points text field, type 0, d1, d1+d2.
5 Click the Build All Objects button.
TR A N S M I S S I O N L I N E ( T L )
Assign the first transmission line the distributed capacitance and inductance C1 and L1,
respectively.
Now define the second transmission line by adding a transmission line equation feature
to the second interval.
Replace the default absorbing boundary condition with lumped ports. With the
lumped ports, it is easy to excite the transmission line and also to plot the S-parameters,
6 |
Q U A R T E R - W A V E TR A N S F O R M E R
that is, the reflection and transmission coefficient, for the transmission line.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 1 only.
Select that this port shall be excited. You can use the default voltage for the port.
3 In the Settings window for Lumped Port, locate the Port Properties section.
4 From the Wave excitation at this port list, choose On.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 3 only.
This lumped port should have a different characteristic impedance than the first
lumped port and the two transmission lines.
3 In the Settings window for Lumped Port, locate the Settings section.
4 In the Zref text field, type ZL.
MESH 1
Let the mesh have a maximum subinterval that is one tenth of the quarter-wave part
of the transmission line.
1 In the Model Builder window, under Component 1 (comp1) click Mesh 1.
2 In the Settings window for Mesh, locate the Mesh Settings section.
3 From the Sequence type list, choose User-controlled mesh.
Size
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Size.
2 In the Settings window for Size, locate the Element Size section.
3 Click the Custom button.
4 Locate the Element Size Parameters section. In the Maximum element size text field,
type hmax.
STUDY 1
7 |
Q U A R T E R - W A V E TR A N S F O R M E R
2 In the Settings window for Frequency Domain, locate the Study Settings section.
3 In the Frequencies text field, type f.
4 On the Model toolbar, click Compute.
RESULTS
Notice that the left part of the curve is flat, with a unit amplitude, indicating that
there are no standing waves, despite the fact that the second lumped port has a load
impedance that normally would not be matched with the transmission line.
STUDY 1
Parameters
1 In the Model Builder window, under Global>Definitions click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
8 |
Expression
Value
Description
df
500[MHz]
5.0000E8 Hz
fstep
10[MHz]
1.0000E7 Hz
Frequency step
Q U A R T E R - W A V E TR A N S F O R M E R
STUDY 1
1D Plot Group 2
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 On the 1D plot group toolbar, click Global.
3 In the Settings window for Global, click Replace Expression in the upper-right corner
of the y-axis data section. From the menu, choose Model>Component 1>Transmission
Line>Ports>S-parameter, dB>tl.S11dB - S-parameter, dB, 11 component.
4 Click Add Expression in the upper-right corner of the y-axis data section. From the
Modify the y-axis limits to show that S21 actually has its maximum value for the
frequency where S11 is at its minimum value.
6 In the Model Builder window, click 1D Plot Group 2.
7 In the Settings window for 1D Plot Group, click to expand the Axis section.
8 Select the Manual axis limits check box.
9 In the y minimum text field, type -30.
10 In the y maximum text field, type 2.
9 |
Q U A R T E R - W A V E TR A N S F O R M E R
11 Click to expand the Legend section. From the Position list, choose Lower right.
12 On the 1D plot group toolbar, click Plot. You should now have a plot of the spectrum
Notice that the curve in the left part of the plot is not flat. The sinusoidal oscillation
in the absolute value of the voltage is a signature of the standing wave that appears
when there is a reflection point along the transmission line. For the selected
frequency the quarter-wave transformer is not a quarter-wave long and, thus, there
are now reflections.
To demonstrate that the quarter-wave transformer not only should have a matched
length, but also a matched characteristic impedance, set the characteristic impedance
of the second lumped port to 50 ohms.
10 |
Q U A R T E R - W A V E TR A N S F O R M E R
TR A N S M I S S I O N L I N E ( T L )
Lumped Port 2
1 In the Model Builder window, under Component 1 (comp1)>Transmission Line (tl) click
Lumped Port 2.
2 In the Settings window for Lumped Port, locate the Settings section.
3 In the Zref text field, type 50.
11 |
Q U A R T E R - W A V E TR A N S F O R M E R
4 On the 1D plot group toolbar, click Plot. Notice that there is now a standing wave
12 |
Q U A R T E R - W A V E TR A N S F O R M E R
1D Plot Group 2
Notice that there is still a resonance at the center frequency. However, as was already
indicated by the spatial plot, there is considerable reflection also at the resonance
frequency.
13 |
Q U A R T E R - W A V E TR A N S F O R M E R
14 |
Q U A R T E R - W A V E TR A N S F O R M E R
1 |
Model Setup
GEOMETRY
Due to symmetry, it is sufficient to model only one quarter of the sphere. Figure 1
shows the geometry and boundary conditions.
PMC symmetry plane
PEC sphere
Figure 1: The computational domain for computing the RCS of a PEC sphere in free
space. Due to symmetry, it is sufficient to model one quarter of the sphere.
The geometry consists of two concentric spherical shells. The innermost shell, adjacent
to the sphere, represents the free space domain, and the second shell represents a
perfectly matched layer (PML) region that is used to provide an approximately
reflection free termination of the, in reality unbounded, free space domain.
EQUATION
The model is set up and solved using a frequency domain formulation for the scattered
electric field. The incident plane wave travels in the positive x direction, with the
electric field polarized along the z-axis. The governing frequency domain equation can
be written in the form
2
( r 1 ( E i + E sc ) ) k 0 rc ( E i + E sc ) = 0
2 |
where the scattered electric field Esc is the dependent variable and the incident electric
field Ei= (0,0,Ez), with
E z = 1 [ V m ]e
jk0 x
The equation is discretized using second order edge elements (also known as vector
elements, Nedelec elements, or curl-conforming elements). It is well known that in
order to resolve the wave field, one should strive for 10 or more discretization points
per wavelength. The combination of using second-order elements and 8 elements per
wavelength fulfills this criterion with some margin. To respect the geometry, a mesh
that is somewhat finer for the longest wavelengths is required on the surface of the
scatterer. A maximum element size of half the radius is used on those boundaries. The
PML region requires special meshing as described under the section Perfectly Matched
Layer below.
BOUNDARY CONDITIONS
The sphere has perfect electric conductor (PEC) boundaries. The PEC boundary
condition
nE = 0
sets the tangential component of the electric field to zero. It is used for the modeling
of lossless metallic surfaces or as a symmetry type boundary condition. It imposes
symmetry for magnetic fields and magnetic currents and antisymmetry for electric
fields and electric currents.
PEC boundary conditions and perfect magnetic conductor (PMC) boundary
conditions apply on the symmetry planes used to subdivide the sphere model.
The PMC boundary condition
nH = 0
sets the tangential component of the magnetic field and thus also the surface current
density to zero. On external boundaries, this can be interpreted as a high surface
impedance boundary condition or used as a symmetry type boundary condition. It
imposes symmetry for electric fields and electric currents and antisymmetry for
magnetic fields and magnetic currents.
PERFECTL Y MATCH ED LAYER
The PML region, the second concentric shell around the sphere, provides an
approximately reflection free termination of the computational domain by applying a
3 |
Figure 2: A free tetrahedral mesh is used in the free-space region around the sphere, and a
swept mesh is used in the PML region.
The free space region around the sphere is defined to be the far-field domain. This
specifies that a near-field to far-field calculation is done on the boundary of this
domain, which takes the computed electric fields around the sphere and uses the
Stratton-Chu equation to find the scattered electric field infinitely far away from the
origin.
In 3D, this is:
jk
E p = ------ r 0 [ n E r 0 ( n H ) ] exp ( jkr r0 ) dS
4
For scattering problems, the far field in COMSOL is identical to what in physics is
known as the scattering amplitude.
4 |
The radiating or scattering object is located in the vicinity of the origin, while the
far-field point p is taken at infinity but with a well-defined angular position ( , ) .
In the above formulas,
E and H are the fields on the aperturethe surface S enclosing the sphere.
r0 is the unit vector pointing from the origin to the field point p. If the field points
lie on a spherical surface S', r0 is the unit normal to S'.
n is the unit normal to the surface S.
is the wave impedance:
=
5 |
Figure 3: Comparison of the analytic solution and the COMSOL Multiphysics model of the
RCS of a PEC sphere in free space.
Mesh Convergence
For the wavelength corresponding to the first maximum in the RCS plot in Figure 3,
a mesh convergence study is performed to validate that the model converges toward a
unique solution when refining the mesh isotropically. The model is solved in a
parametric sweep over the number of mesh elements per wavelength. In the PML, the
mesh density is not changed in the radial outward direction (that is, in the sweep
direction for the swept mesh). The PML is resolved by 5 element layers in this
direction which is sufficient to resolve the exponential damping in the radial direction.
Thus the error contribution from the PML is not expected to decrease by adding more
element layers. The main error contribution from the PML is due to the fact that it is
not perfectly absorbing because of finite thickness and damping rather than mesh
density. Thus, it is expected to give a contribution to the error in the computed RCS
that will not decrease when refining the mesh.
Figure 4 shows the mesh convergence. The displayed error is the difference between
the RCS from the finite element model and the exact solution from equation 11-247
6 |
Figure 4: Mesh convergence for the difference in backscattering (monostatic) RCS between
the COMSOL model and the exact solution.
Reference
1. C.A. Balanis, Advanced Engineering Electromagnetics, John Wiley & Sons, 1989.
Modeling Instructions
From the File menu, choose New.
7 |
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
r_lda
0.5
0.5000
r0
5[cm]
0.05000 m
Sphere radius
lda
r0/r_lda
0.1000 m
Wavelength
k0
2*pi/lda
62.83 1/m
Wavenumber
f0
c_const/lda
2.998E9 1/s
Frequency
t_air
lda/2
0.05000 m
t_pml
lda/2
0.05000 m
Thickness of PML
h_size
8.000
E0
1[V/m]
1.000 V/m
GEOMETRY 1
First, create a sphere with two layer definitions. The outermost layer represents the
PMLs and the core represents the PEC sphere for RCS analysis. The median layer is
the air domain.
Sphere 1 (sph1)
1 On the Geometry toolbar, click Sphere.
8 |
Thickness (m)
Layer 1
t_pml
Layer 2
t_air
choose View.
Camera
Change only the sign of y in the Position and Up Vector sections:
1 In the Model Builder window, expand the View 2 node, then click Camera.
2 In the Settings window for Camera, locate the Position section.
3 In the y text field, type 1.871.
4 Locate the Up Vector section. In the y text field, type -0.412.
5 Click the Update button.
9 |
GEOMETRY 1
Due to the symmetry of the structure, it is sufficient to model only one quarter of the
sphere. Delete the domains which are not part of the modeling domain.
section.
3 From the Geometric entity level list, choose Domain.
10 |
4 On the object sph1, select Domains 13, 57, and 915 only.
View 1
After removing unnecessary domains, change the view to the first view definition
which gives a better angle showing all layers.
11 |
Now set up the physics. You will solve the model for the scattered field, which requires
background electric field (E-field) information. The background plane wave is
traveling in the positive x direction, with the electric field polarized along the z-axis.
The default boundary condition is perfect electric conductor, which applies to all
exterior boundaries including the boundaries perpendicular to the background E-field
polarization.
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic
Waves, Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Settings section.
3 From the Solve for list, choose Scattered field.
4 Specify the Eb vector as
12 |
E0*exp(-j*k0*x)
5 Locate the Analysis Methodology section. From the Methodology options list, choose
Fast.
Choose wireframe rendering in the current view to get a better view of the interior
parts.
6 Click the Wireframe Rendering button on the Graphics toolbar.
DEFINITIONS
The outermost domains from the center of the sphere are the PMLs.
13 |
3 In the Settings window for Perfectly Matched Layer, locate the Geometry section.
4 From the Type list, choose Spherical.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
14 |
Far-Field Domain 1
1 On the Physics toolbar, click Domains and choose Far-Field Domain.
2 Select Domains 2 and 3 only.
15 |
Far-Field Calculation 1
1 In the Model Builder window, expand the Far-Field Domain 1 node, then click Far-Field
Calculation 1.
2 In the Settings window for Far-Field Calculation, locate the Boundary Selection
section.
3 Click Clear Selection.
4 Select Boundaries 6 and 13 only.
5 Locate the Far-Field Calculation section. Select the Symmetry in the y=0 plane check
box.
6 Select the Symmetry in the z=0 plane check box.
7 From the Symmetry type list, choose Symmetry in H (PEC).
MATERIALS
16 |
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
MESH 1
Free Tetrahedral 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
The maximum mesh size is at most 0.2 wavelengths in free space. In this model, use
0.125 wavelengths.
Size
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Size.
2 In the Settings window for Size, locate the Element Size section.
3 Click the Custom button.
4 Locate the Element Size Parameters section. In the Maximum element size text field,
type lda/h_size.
17 |
Swept 1
In the Model Builder window, right-click Mesh 1 and choose Swept.
Distribution 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Swept
1 and choose Distribution.
2 Right-click Distribution 1 and choose Build All. Compare the mesh with that shown
in Figure 2.
STUDY 1
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
r_lda
range(0.1,0.025,0.8)
Parameter unit
Follow the instructions below to reproduce the plot in Figure 3 First, show the
computed RCS values using square markers.
1D Plot Group 1
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
18 |
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Study 1/Parametric Solutions 1.
4 Click to expand the Title section. From the Title type list, choose Custom.
5 Find the Type and data subsection. Clear the Unit check box.
6 Clear the Description check box.
7 Clear the Type check box.
8 Find the User subsection. In the Prefix text field, type RCS Calculation.
9 Locate the Plot Settings section. Select the x-axis label check box.
10 In the associated text field, type Sphere radius in wavelengths (a/
lambda<sub>0</sub>).
13 Locate the Axis section. Select the y-axis log scale check box.
14 On the 1D plot group toolbar, click Point Graph.
15 Select Point 2 only.
16 In the Settings window for Point Graph, locate the y-Axis Data section.
17 In the Expression text field, type 4*pi*(emw.normEfar^2/E0^2)/(pi*r0^2).
18 Locate the x-Axis Data section. From the Axis source data list, choose Outer solutions.
19 Click to expand the Coloring and style section. Locate the Coloring and Style section.
Find the Line style subsection. From the Line list, choose None.
20 Find the Line markers subsection. From the Marker list, choose Square.
21 From the Positioning list, choose In data points.
19 |
Parameters
1 In the Model Builder window, under Global>Definitions click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 Add the two rows at the end and change the first line according to:
20 |
Name
Expression
Value
Description
r_lda
r1
0.1636
Sphere radius in
wavelengths
r0
5[cm]
0.05000 m
Sphere radius
lda
r0/r_lda
0.3056 m
Wavelength
Name
Expression
Value
Description
k0
2*pi/lda
20.56 1/m
Wavenumber
f0
c_const/lda
9.811E8 1/s
Frequency
t_air
lda/2
0.1528 m
t_pml
lda/2
0.1528 m
Thickness of PML
h_size
8.000
E0
1[V/m]
1.000 V/m
r1
0.16363636363636364
0.1636
RCS1
3.6549540474068576
3.655
MESH 1
Add a new mesh with some tweaks to make sure that the curvature of the sphere is
always resolved.
This is to avoid inverted mesh elements.
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Duplicate.
COMPONENT 1 (COMP1)
Size
1 In the Model Builder window, expand the Component 1 (comp1)>Meshes>Mesh 2
Size 1
1 In the Model Builder window, under Component 1 (comp1)>Meshes>Mesh 2
21 |
box.
10 In the associated text field, type r0/2.
ROOT
Add a new frequency domain study for the mesh convergence analysis.
ADD STUDY
1 On the Model toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Preset Studies>Frequency
Domain.
4 Click Add Study in the window toolbar.
5 On the Model toolbar, click Add Study to close the Add Study window.
STUDY 2
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
h_size
3 6 9 12 15 20
Parameter unit
22 |
1D Plot Group 2
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Study 2/Parametric Solutions 2.
4 Locate the Plot Settings section. Select the x-axis label check box.
5 In the associated text field, type Number of elements / lambda<sub>0</sub>.
6 Select the y-axis label check box.
7 In the associated text field, type Relative error.
8 On the 1D plot group toolbar, click Point Graph.
9 In the Settings window for Point Graph, locate the Selection section.
10 Click Paste Selection.
11 In the Paste Selection dialog box, type 2 in the Selection text field.
12 Click OK.
13 In the Settings window for Point Graph, locate the y-Axis Data section.
14 In the Expression text field, type (4*pi*(emw.normEfar^2/E0^2)/
(pi*r0^2)-RCS1)/RCS1.
15 Locate the x-Axis Data section. From the Axis source data list, choose h_size.
16 From the Parameter list, choose Expression.
17 In the Expression text field, type h_size.
18 Click to expand the Coloring and style section. Locate the Coloring and Style section.
Find the Line style subsection. In the Width text field, type 2.
19 Find the Line markers subsection. From the Marker list, choose Square.
20 From the Positioning list, choose In data points.
21 Click the x-Axis Log Scale button on the Graphics toolbar.
23 |
22 Click the y-Axis Log Scale button on the Graphics toolbar. Compare the convergence
24 |
RF Heating
Introduction
This is a model of an RF waveguide bend with a dielectric block inside. There are
electromagnetic losses in the block as well as on the waveguide walls which cause the
assembly to heat up over time. The material properties of the block are functions of
temperature. The transient thermal behavior, as well as the steady-state solution, are
computed.
Dielectric block
Copper coating
H-bend waveguide
Figure 1: A waveguide bend with a dielectric block inside. Top boundaries of the
waveguide are removed only for visualization.
Model Definition
The waveguide bend shown in Figure 1 is connected to a 100 W power source,
operating at 10 GHz, via a rectangular waveguide operating in the TE10 mode. The
other end of the bend is also connected to a rectangular waveguide operating in the
TE10 mode. The objective of such a bend is primarily to change the direction of
1 |
RF HEATING
2 |
RF HEATING
transfer equation is, on the other hand, is solved transiently. The electromagnetic fields
are only recomputed when the material properties have changed significantly, as
determined by a criterion involving the relative tolerance of the time-dependent solver.
The objective of the analysis is to determine the change in temperature from given
initial conditions and how long these changes take.
A frequency-stationary simulation solves Maxwells equations in the frequency
domain, but it solves the stationary heat transfer equation under the assumption that
all initial transient variations have died out. Although no transient information is
obtainable, this computation is significantly faster than a frequency-transient analysis
and gives the steady-state temperature distribution.
Figure 2: The maximum temperature, evaluated over the volume of the block, is plotted as
a function of temperature.
3 |
RF HEATING
Figure 3: The electric fields (red arrows) magnetic fields (green arrows) and power flow
(blue arrows) are shown inside of the waveguide. The steady-state temperature is plotted
on the block and waveguide walls.
Figure 4: The loss tangent within the dielectric block for the steady-state solution shows that
the variation in temperature affects the material properties.
4 |
RF HEATING
Modeling Instructions
From the File menu, choose New.
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
f0
10[GHz]
1.0000E10 Hz
Current frequency
lda0
c_const/f0
0.029979 m
Wavelength, air
h_max
0.2*lda0
0.0059958 m
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
5 |
RF HEATING
GEOMETRY 1
First, import the geometry of the waveguide including a dielectric block inside the
waveguide.
Import 1 (imp1)
1 On the Model toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 Click Browse.
4 Browse to the models Model Library folder and double-click the file
rf_heating.mphbin.
5 Click Import.
Use the wireframe rendering to see the inner parts of the waveguide.
6 Click the Wireframe Rendering button on the Graphics toolbar.
DEFINITIONS
Create a set of selections for use before setting up the physics. First, create a selection
for the dielectric block.
Explicit 1
1 On the Definitions toolbar, click Explicit.
6 |
RF HEATING
Explicit 2
1 On the Definitions toolbar, click Explicit.
7 |
RF HEATING
Explicit 3
1 On the Definitions toolbar, click Explicit.
8 |
RF HEATING
Explicit 4
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, locate the Input Entities section.
3 From the Geometric entity level list, choose Boundary.
9 |
RF HEATING
4 Select Boundaries 1618, 35, 53, 54, 72, 74, 75, 78, 96, and 97 only.
Explicit 5
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, locate the Input Entities section.
3 From the Geometric entity level list, choose Boundary.
10 |
RF HEATING
To get a better view, suppress some of the boundaries. Furthermore, by assigning the
resulting settings to a View node, you can easily return to the same view later by
clicking the Go to View 2 button on the Graphics toolbar.
View 2
1 In the Model Builder window, right-click Definitions and choose View.
2 Click the Wireframe Rendering button on the Graphics toolbar.
3 On the 3D view toolbar, click Hide Geometric Entities.
4 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
5 From the Geometric entity level list, choose Boundary.
11 |
RF HEATING
On the Physics toolbar, click Electromagnetic Waves, Frequency Domain (emw) and
choose Heat Transfer in Solids (ht).
12 |
RF HEATING
1 In the Model Builder window, under Component 1 (comp1) click Heat Transfer in Solids
(ht).
2 In the Settings window for Heat Transfer in Solids, locate the Domain Selection
section.
3 In the list, select 2.
4 Click Remove from Selection.
5 Select Domains 1 and 3 only.
Heat Flux 1
1 On the Physics toolbar, click Boundaries and choose Heat Flux.
2 In the Settings window for Heat Flux, locate the Boundary Selection section.
3 From the Selection list, choose Waveguide outside surfaces.
4 Locate the Heat Flux section. Click the Convective heat flux button.
5 In the h text field, type 5.
ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN (EMW)
section.
3 From the Selection list, choose Air.
13 |
RF HEATING
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 Select Boundary 15 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
5 From the Wave excitation at this port list, choose On.
6 In the Pin text field, type 100.
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
14 |
RF HEATING
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
MATERIALS
Next, assign material properties on the model. Begin by specifying Aluminum for the
waveguide structure.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Aluminum.
4 Click Add to Component in the window toolbar.
MATERIALS
Aluminum (mat1)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Aluminum
(mat1).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Waveguide.
15 |
RF HEATING
ADD MATERIAL
Air (mat2)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Air (mat2).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Air.
Material 3 (mat3)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 Select Domain 3 only.
3 In the Settings window for Material, locate the Material Contents section.
4 In the table, enter the following settings:
Property
Name
Value
Unit
Property group
epsilo
nPrim
2.1
Dielectric losses
delta
delta
0.001*
(T/
300[K]
)
Dielectric losses
Relative permeability
mur
Basic
Thermal conductivity
0.3[W/
m/
K]*(T/
300[K]
)
W/
(mK)
Basic
Density
rho
2200
kg/m
Basic
Cp
1050
J/
(kgK)
Basic
16 |
RF HEATING
ADD MATERIAL
1 On the Model toolbar, click Add Material to close the Add Material window.
Copper (mat4)
1 In the Model Builder window, under Component 1 (comp1)>Materials click Copper
(mat4).
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Boundary.
4 From the Selection list, choose Waveguide inside surfaces.
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
17 |
RF HEATING
STUDY 1
Step 1: Frequency-Transient
1 In the Model Builder window, under Study 1 click Step 1: Frequency-Transient.
2 In the Settings window for Frequency-Transient, locate the Study Settings section.
3 In the Times text field, type range(0,15,300).
4 Select the Relative tolerance check box.
5 In the associated text field, type 0.001.
6 In the Frequency text field, type f0.
7 In the Model Builder window, click Study 1.
8 In the Settings window for Study, locate the Study Settings section.
9 Clear the Generate default plots check box.
10 On the Model toolbar, click Compute.
RESULTS
Data Sets
Plot the transient response of the peak temperature.
1 On the Results toolbar, click More Data Sets and choose Evaluation>Maximum.
1D Plot Group 1
1 On the Results toolbar, click 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Maximum 1.
4 On the 1D plot group toolbar, click Point Graph.
5 In the Settings window for Point Graph, click Replace Expression in the upper-right
corner of the y-axis data section. From the menu, choose Model>Component 1>Heat
Transfer in Solids>Temperature>T - Temperature.
6 Locate the y-Axis Data section. From the Unit list, choose degC.
7 On the 1D plot group toolbar, click Plot.
18 |
RF HEATING
ADD STUDY
1 On the Model toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
3 Find the Studies subsection. In the Select study tree, select Custom Studies>Empty
Study.
4 Click Add Study in the window toolbar.
5 On the Model toolbar, click Add Study to close the Add Study window.
STUDY 2
Step 1: Frequency-Stationary
1 On the Study toolbar, click Study Steps and choose Stationary>Frequency-Stationary.
2 In the Settings window for Frequency-Stationary, locate the Study Settings section.
3 In the Frequency text field, type f0.
4 On the Study toolbar, click Compute.
RESULTS
Temperature (ht)
The default plots show the distribution of the electric field norm and the temperature.
For the temperature plot, first change the unit to the degree Celsius and then, add
arrow plots of the electric fields, magnetic fields, and power flow.
1 In the Model Builder window, expand the Temperature (ht) node, then click Surface 1.
2 In the Settings window for Surface, locate the Expression section.
3 From the Unit list, choose degC.
4 In the Model Builder window, right-click Temperature (ht) and choose Arrow Volume.
5 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
19 |
RF HEATING
choose Rename.
2 In the Rename 3D Plot Group dialog box, type Loss Tangent (emw) in the New label
text field.
3 Click OK.
in.
20 |
RF HEATING
Substrate
Lumped port
Microstrip line
Metallic post
Air-filled metallic cavity
Figure 1: A tunable evanescent mode cavity filter is composed of a rectangular cavity with
a metallic post, a piezo actuator disk, and slot-coupled microstrip lines. There is a small
gap between the top of the post and the bottom side of the piezo actuator. The front part of
the cavity wall is removed for visualization purposes.
Note: In addition to the RF Module, this model requires one of the Acoustics
Module, the MEMS Module, or the Structural Mechanics Module.
1 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Model Definition
This example starts from a basic rectangular cavity filter, whose resonant frequencies
are given by
c
m 2
n 2
l 2
f nml = ---------------------- -------- + ------- + -----
a
b
d
2 r r
(1)
where a and b are the waveguide aperture dimensions and d is the length of the
waveguide cavity. In this model, the cavity width, height, and length are a = 100 mm,
b = 50 mm, and d = 50 mm, respectively. The resulting resonant frequency of the
dominant mode, TE101, is 3.354 GHz.
By adding a metallic post and creating reactance inside the cavity, the resonance
frequency can be lowered. The cavity is air filled and the height of the post is slightly
smaller than b, which creates a small gap between the top of the post and the cavity
where the electric fields are confined. Two shorted 50 microstrip lines on a dielectric
substrate, fed by a lumped port, are coupled into the cavity. The dimensions and
locations of the slots can be adjusted to improve input matching properties and power
transfer between input and output ports. The air box around the microstrip lines are
enclosed by a scattering boundary condition representing the infinite air space. A
circular aperture at the top of the cavity is closed with a piezo actuator and the bottom
surface of the disk is finished with a layer of a highly conductive material that is several
skin depths in thickness.
Model all metal partsthe cavity walls, post, substrate ground planes, microstrip lines,
and the bottom surface of the piezo deviceas perfect electric conductors (PECs).
The material for the piezo actuator is Lead Zirconate Titanate (PZT-5H). It is
z-polarized and generates mainly z-directional deflection of the device.
Mesh the model using a tetrahedral mesh with approximately five elements per
wavelength in each material at the highest simulation frequency. When the piezo
device deforms due to the input bias, the Moving Mesh interface is used to deform the
mesh for the Electromagnetic Waves physics.
2 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
between the top of the post and the bottom of the piezo device, strong electric fields
are observed.
Figure 2: This plot shows the total piezo displacement when 300 V is applied on the
actuator. The visualization is exaggerated to emphasize the deflection.
The S-parameters plotted in Figure 4 show the effect of the piezo device deflection on
the filters resonant frequency. The tunable frequency range of this model is ~40MHz.
This range can be adjusted by different choices of the piezo disk size and the input bias
voltage.
3 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Figure 3: The dominant mode inside the cavity is observed from the electric field
distribution plot.
Figure 4: The deflection for the piezo device controlled by the input bias can shift the
resonance frequency of the filter.
4 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Modeling Instructions
From the File menu, choose New.
NEW
5 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
following steps.
1 On the Geometry toolbar, click Insert Sequence.
2 Browse to the models Model Library folder and double-click the file
tunable_cavity_filter_geom_sequence.mph.
3 In the Model Builder window, under Component 1 (comp1) right-click Geometry 1 and
DEFINITIONS
Parameters
1 In the Model Builder window, expand the Global>Definitions node, then click
Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
6 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Expression
Value
Description
V0
300[V]
300.00 V
f_min
3[GHz]
3.0000E9 Hz
f_max
3.06[GHz]
3.0600E9 Hz
ldao
c_const/f_max
0.097971 m
h_max
0.2*ldao
0.019594 m
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
Create a set of selections for use when setting up the physics. First, create a selection
for the microstrip feed line.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Feed line in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundaries 17, 22, 47, and 50 only.
7 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Ground in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundaries 6, 16, 28, 39, and 53 only.
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Substrate in the Label text field.
8 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Explicit 4
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Piezo actuator in the Label text field.
9 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Explicit 5
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Piezo fixed edges in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose Edge.
10 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Explicit 6
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Open boundaries in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
11 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
4 Select Boundaries 4, 5, 7, 8, 10, 12, 13, 23, 25, 37, 38, 40, 41, 4345, 56, and 57
only.
1 In the Model Builder window, under Component 1 (comp1) click Solid Mechanics
(solid).
2 In the Settings window for Solid Mechanics, locate the Domain Selection section.
3 From the Selection list, choose Piezo actuator.
ELECTROSTATICS (ES)
On the Physics toolbar, click Solid Mechanics (solid) and choose Electrostatics (es).
1 In the Model Builder window, under Component 1 (comp1) click Electrostatics (es).
2 In the Settings window for Electrostatics, locate the Domain Selection section.
3 From the Selection list, choose Piezo actuator.
Set up the Moving Mesh interface. Because the substrate and air domains are
deflection free regions, do not include them in this physics.
MOVING MESH (ALE)
On the Physics toolbar, click Electrostatics (es) and choose Moving Mesh (ale).
12 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
1 In the Model Builder window, under Component 1 (comp1) click Moving Mesh (ale).
2 In the Settings window for Moving Mesh, locate the Domain Selection section.
3 Click Clear Selection.
4 Select Domain 1 only.
Set up the Electromagnetic Waves, Frequency Domain interface. Suppress the piezo
actuator disk domain.
ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN (EMW)
On the Physics toolbar, click Moving Mesh (ale) and choose Electromagnetic Waves,
Frequency Domain (emw).
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic
Waves, Frequency Domain (emw).
13 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Assign material properties. Use three materials for this model: PZT-5H, air, and a
user-defined substrate.
MATERIALS
On the Model toolbar, click More Windows and choose Add Material.
ADD MATERIAL
14 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
ADD MATERIAL
Air (mat2)
1 On the Model toolbar, click Add Material to close the Add Material window.
2 In the Model Builder window, under Component 1 (comp1)>Materials click Air (mat2).
3 Select Domains 1, 3, and 7 only.
Material 3 (mat3)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Substrate.
4 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
3.38
Basic
15 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Property
Name
Value
Unit
Property group
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Now, set up the physics constraints. Start by assuming the bottom rim of piezo
actuator part is attached on the same size circular aperture of the cavity top and no
deflection is expected.
Fixed Constraint 1
1 On the Physics toolbar, click Edges and choose Fixed Constraint.
2 In the Settings window for Fixed Constraint, locate the Edge Selection section.
3 From the Selection list, choose Piezo fixed edges.
ELECTROSTATICS (ES)
Electric Potential 1
1 On the Physics toolbar, click Boundaries and choose Electric Potential.
2 Select Boundary 29 only.
3 In the Settings window for Electric Potential, locate the Electric Potential section.
16 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Ground 1
1 On the Physics toolbar, click Boundaries and choose Ground.
2 In the Settings window for Ground, locate the Boundary Selection section.
3 From the Selection list, choose Ground.
MOVING MESH (ALE)
Free Deformation 1
1 On the Physics toolbar, click Domains and choose Free Deformation.
2 Select Domain 1 only.
17 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
section.
3 From the Selection list, choose Feed line.
section.
3 From the Selection list, choose Ground.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 14 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
18 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 54 only.
19 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Before generating the mesh, set the physics to control the mesh size.
1 In the Model Builder window, click Electromagnetic Waves, Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
4 In the Maximum element size text field, type h_max.
5 In the Model Builder window, click Electromagnetic Waves, Frequency Domain (emw).
6 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Analysis Methodology section.
7 From the Methodology options list, choose Fast.
MESH 1
Size 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Size 1.
2 In the Settings window for Size, locate the Element Size Parameters section.
3 In the Maximum element size text field, type 10.
4 In the Minimum element size text field, type 2.
Free Tetrahedral 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Free
Tetrahedral 1.
2 In the Settings window for Free Tetrahedral, locate the Domain Selection section.
3 From the Geometric entity level list, choose Domain.
4 From the Selection list, choose All domains.
5 Select Domains 14 and 68 only.
Swept 1
1 In the Model Builder window, right-click Mesh 1 and choose Swept.
2 In the Settings window for Swept, locate the Domain Selection section.
3 From the Geometric entity level list, choose Domain.
4 Select Domain 5 only.
20 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Distribution 1
1 Right-click Component 1 (comp1)>Mesh 1>Swept 1 and choose Distribution.
2 In the Settings window for Distribution, locate the Distribution section.
3 In the Number of elements text field, type 2.
4 Click the Build All button.
STUDY 1
Step 1: Stationary
1 On the Study toolbar, click Study Steps and choose Stationary>Stationary.
2 In the Settings window for Stationary, locate the Physics and Variables Selection
section.
3 In the table, enter the following settings:
Physics interface
Solve for
Discretization
physics
21 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
settings:
Physics interface
Solve for
Discretization
Solid Mechanics
physics
Electrostatics
physics
Moving Mesh
physics
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
V0
-300,300
Parameter unit
Stress (solid)
Replace the default stress plot by displacement plot.
1 In the Settings window for 3D Plot Group, type Displacement in the Label text
field.
Displacement
1 In the Model Builder window, expand the Results>Displacement node, then click
Surface 1.
2 In the Settings window for Surface, click Replace Expression in the upper-right corner
22 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
This plot shows the deflected piezo actuator disk; compare with Figure 2.
23 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
The resulting plot shows strong electric fields resulting from the dominant
resonance at the center of the cavity as well as in the gap between the metallic post
and the ceiling of the cavity.Compare the plot with that shown in Figure 3.
S-Parameter (emw)
1 In the Settings window for 1D Plot Group, click to expand the Legend section.
2 From the Position list, choose Lower left.
3 In the Model Builder window, expand the S-Parameter (emw) node, then click Global
1.
4 In the Settings window for Global, locate the x-Axis Data section.
5 From the Unit list, choose GHz.
6 On the 1D plot group toolbar, click Plot.
The plotted S-parameters show the frequency shift as a function of the input bias on
the piezo actuator; compare with Figure 4.
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 Click Load from File.
4 Browse to the models Model Library folder and double-click the file
tunable_cavity_filter_parameters.txt.
GEOMETRY 1
24 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Cavity in the Label text field.
3 Locate the Size section. In the Width text field, type 100.
4 In the Depth text field, type 50.
5 In the Height text field, type 50.
6 Locate the Position section. From the Base list, choose Center.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Substrate in the Label text field.
3 Locate the Size section. In the Width text field, type 25.
4 In the Depth text field, type 50.
5 In the Height text field, type thickness.
6 Locate the Position section. From the Base list, choose Center.
7 In the x text field, type -37.5.
8 In the z text field, type 25+thickness/2.
Block 3 (blk3)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Air block in the Label text field.
3 Locate the Size section. In the Width text field, type 25.
4 In the Depth text field, type 50.
5 In the Height text field, type 10.
6 Locate the Position section. From the Base list, choose Center.
7 In the x text field, type -37.5.
8 In the z text field, type 30.
9 Click the Build Selected button.
25 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Block 4 (blk4)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Feed_line in the Label text field.
3 Locate the Size section. In the Width text field, type l_feed+w_slot.
4 In the Depth text field, type 3.2.
5 In the Height text field, type thickness.
6 Locate the Position section. From the Base list, choose Center.
7 In the x text field, type -x_slot-l_feed/2.
8 In the z text field, type 25+thickness/2.
26 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Plane Geometry
Add a rectangle for the slot.
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type w_slot.
4 In the Height text field, type l_slot.
5 Locate the Position section. From the Base list, choose Center.
6 In the xw text field, type -x_slot.
7 Click the Build Selected button.
8 Click the Zoom Extents button on the Graphics toolbar.
Generate the 2nd slot coupled microstrip line by mirroring some geometries.
Mirror 1 (mir1)
1 On the Geometry toolbar, click Transforms and choose Mirror.
2 Select the objects wp1, blk3, blk2, and blk4 only.
3 In the Settings window for Mirror, locate the Input section.
4 Select the Keep input objects check box.
5 Locate the Normal Vector to Plane of Reflection section. In the x text field, type 1.
27 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
Add a block for the metal post in the middle of the cavity.
Block 5 (blk5)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Post in the Label text field.
3 Locate the Size section. In the Width text field, type 15.
4 In the Depth text field, type 15.
5 In the Height text field, type 50-gap_post.
6 Locate the Position section. In the x text field, type -7.5.
7 In the y text field, type -7.5.
8 In the z text field, type -25.
9 Click the Build Selected button.
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type Piezo actuator in the Label text field.
3 Locate the Size and Shape section. In the Radius text field, type 21.
4 In the Height text field, type 0.5.
5 Locate the Position section. In the z text field, type 25.
6 Click the Build Selected button.
The inside of the metal post is not part of the modeling domain. Therefore, subtract
it from the cavity.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object blk1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the object blk5 only.
6 Click the Build All Objects button.
7 Click the Zoom Extents button on the Graphics toolbar.
28 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
29 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
30 |
TU N A B L E E V A N E S C E N T M O D E C A V I T Y F I L T E R U S I N G A P I E Z O A C T U A T O R
V i v a ldi Ante nn a
Introduction
A tapered slot antenna, also known as a Vivaldi antenna, is useful for wide-band
applications. Here, an exponential function is used for the taper profile. The objective of
this model is to compute the far-field pattern and to compute the impedance of the
structure. Good matching is observed over a wide frequency band.
Shorted
Substrate
with PEC ground
Lumped port
on a microstrip line
Figure 1: The Vivaldi antenna is realized on a thin dielectric substrate. The entire domain is
bounded by a perfectly matched layer.
Model Definition
In this Vivaldi antenna model, the tapered slot is patterned with a perfect electric
conductor (PEC) ground plane on the top of the dielectric substrate. A simple exponential
function, e0.044x is used to create the tapered slot curves. One end of the slot is open to
air and the other end is finished with a circular slot. On the bottom of the substrate, the
shorted 50 microstrip feed line is modeled as PEC surfaces. The entire modeling
domain is bounded by a perfectly matched layer (PML) which acts like an anechoic
1 |
VIVALDI ANTENNA
chamber absorbing all radiated energy. To excite the antenna, a lumped port is used. The
model is meshed using a tetrahedral mesh with approximately five elements per wavelength
in each material and simulation frequency.
Figure 2: The frequency response SWR of the Vivaldi antenna shows wide-band impedance
matching, better than 2:1 in most of the simulated frequency range.
2 |
VIVALDI ANTENNA
References
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
2. C.A. Balanis, Antenna Theory, John Wiley & Sons, 1997.
Modeling Instructions
From the File menu, choose New.
NEW
3 |
VIVALDI ANTENNA
MODEL WIZARD
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
thickness
60[mil]
0.001524 m
Substrate thickness
w_slot
0.5[mm]
5.000E-4 m
Slot with
f_min
2.0[GHz]
2.000E9 Hz
Minimum frequency in
sweep
f_max
6.5[GHz]
6.500E9 Hz
Maximum frequency in
sweep
f0
f_max
6.500E9 Hz
Current frequency in
sweep
lda0
c_const/f0
0.04612 m
h_max
0.2*lda0
0.009224 m
Here, 'mil' refers to the unit milliinch and c_const is a predefined COMSOL constant
for the speed of light in vacuum.
4 |
VIVALDI ANTENNA
GEOMETRY 1
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 110.
4 In the Depth text field, type 80.
5 In the Height text field, type thickness.
6 Locate the Position section. From the Base list, choose Center.
7 Right-click Component 1 (comp1)>Geometry 1>Block 1 (blk1) and choose Rename.
8 In the Rename Block dialog box, type Substrate in the New label text field.
9 Click OK.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 3.2.
4 In the Depth text field, type 40+w_slot/2.
5 In the Height text field, type thickness.
6 Locate the Position section. From the Base list, choose Center.
7 In the x text field, type -26.
8 In the y text field, type -20+w_slot/4.
9 Right-click Component 1 (comp1)>Geometry 1>Block 2 (blk2) and choose Rename.
10 In the Rename Block dialog box, type Feed line in the New label text field.
11 Click OK.
Next, create a work plane where you will draw the Vivaldi antenna pattern. Use two
parametric curves for the tapered slot.
5 |
VIVALDI ANTENNA
Plane Geometry
1 Click the Zoom Extents button on the Graphics toolbar.
Mirror 1 (mir1)
1 On the Work plane toolbar, click Transforms and choose Mirror.
2 In the Settings window for Mirror, locate the Normal Vector to Line of Reflection section.
3 In the yw text field, type 1.
4 In the xw text field, type 0.
5 Locate the Input section. Select the Keep input objects check box.
6 Select the object pc1 only.
7 Click the Build Selected button.
6 |
VIVALDI ANTENNA
Add a rectangle describing the thin slot connected to the tapered slot.
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 20.
4 In the Height text field, type w_slot.
5 Locate the Position section. In the xw text field, type -35.
6 In the yw text field, type -w_slot/2.
Circle 1 (c1)
1 On the Work plane toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type 12.
4 Locate the Position section. In the xw text field, type -40.5.
Create a union of the circle and the rectangle to remove unnecessary boundaries.
Union 1 (uni1)
1 On the Work plane toolbar, click Booleans and Partitions and choose Union.
2 Select the objects r1 and c1 only.
3 In the Settings window for Union, locate the Union section.
4 Clear the Keep interior boundaries check box.
7 |
VIVALDI ANTENNA
Add a sphere for the PMLs. Use a layer definition to create a shell-type structure.
Sphere 1 (sph1)
1 On the Geometry toolbar, click Sphere.
2 In the Settings window for Sphere, locate the Size section.
3 In the Radius text field, type 110.
4 Click to expand the Layers section. In the table, enter the following settings:
Layer name
Thickness (mm)
Layer 1
30
8 |
VIVALDI ANTENNA
DEFINITIONS
View 1
Hide some domains to get a better view of the interior parts when setting up the physics
and reviewing the mesh.
1 On the 3D view toolbar, click Hide Geometric Entities.
9 |
VIVALDI ANTENNA
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic Waves,
Frequency Domain (emw).
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
Now set up the physics. Use the selections already defined when assigning boundary
conditions.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
10 |
VIVALDI ANTENNA
5 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
6 From the Wave excitation at this port list, choose On.
MATERIALS
Assign material properties for the model. First, use air for all domains.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
1 Go to the Add Material window.
2 In the tree, select Built-In>Air.
3 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
11 |
VIVALDI ANTENNA
Name
Value
Unit
Property group
Relative permittivity
epsilonr
3.38
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
1 In the Settings window for Mesh, locate the Mesh Settings section.
2 From the Element size list, choose Coarse.
3 Click the Build All button.
12 |
VIVALDI ANTENNA
STUDY 1
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
Parameter name
f0
range(f_min,0.5[GHz],f_max)
Parameter unit
type 0.
5 Find the y-planes subsection. In the Planes text field, type 0.
13 |
VIVALDI ANTENNA
Strong electric fields are observed in the slot and microstrip line.
14 |
VIVALDI ANTENNA
3D Plot Group 3
1 In the Model Builder window, under Results click 3D Plot Group 3.
2 In the Settings window for 3D Plot Group, locate the Data section.
3 From the Parameter value (f0 (1/s)) list, choose 5.500E9.
4 On the 3D plot group toolbar, click Plot.
5 Click the Zoom Extents button on the Graphics toolbar.
1D Plot Group 4
1 On the Model toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Study 1/Parametric Solutions 1.
4 On the 1D plot group toolbar, click Global.
5 In the Settings window for Global, locate the x-Axis Data section.
6 From the Axis source data list, choose Outer solutions.
7 Click to expand the Legends section. Clear the Show legends check box.
15 |
VIVALDI ANTENNA
8 Locate the y-Axis Data section. In the table, enter the following settings:
Expression
Unit
Description
(1+abs(emw.S11))/
(1-abs(emw.S11))
SWR
This SWR plot replicates the wide-band frequency response shown in Figure 2.
16 |
VIVALDI ANTENNA
Waveguide Adapter
Introduction
This is a model of an adapter for microwave propagation in the transition between a
rectangular and an elliptical waveguide. Such waveguide adapters are designed to keep
energy losses due to reflections at a minimum for the operating frequencies. To
investigate the characteristics of the adapter, the simulation includes a wave traveling
from a rectangular waveguide through the adapter and into an elliptical waveguide.
The S-parameters are calculated as functions of the frequency. The involved
frequencies are all in the single-mode range of the waveguide, that is, the frequency
range where only one mode is propagating in the waveguide.
Model Definition
The waveguide adapter consists of a rectangular part smoothly transcending into an
elliptical part as seen in Figure 1.
1 |
WAVE G U I D E A D AP T E R
The rectangular port is excited by a transverse electric (TE) wave, which is a wave that
has no electric field component in the direction of propagation. This is what an
incoming wave would look like after traveling through a straight rectangular
waveguide with the same cross section as the rectangular part of the adapter. The
excitation frequencies are selected so that the TE10 mode is the only propagating mode
through the rectangular waveguide. The cutoff frequencies for the different modes can
be achieved analytically from the relation
2
c m
n
( c ) mn = --- ----- + ---
2 a
b
where m and n are the mode numbers, and c is the speed of light. For the TE10 mode,
m = 1 and n = 0. With the dimensions of the rectangular cross section (a = 2.286 cm
and b = 1.016 cm), the TE10 mode is the only propagating mode for frequencies
between 6.6 GHz and 14.7 GHz.
Although the shape of the TE10 mode is known analytically, this model lets you
compute it using a numerical port. This technique is very general, in that it allows the
port boundary to have any shape. The solved equation is
2
2 2
( n H n ) + ( n k 0 )H n = 0
Here Hn is the component of the magnetic field perpendicular to the boundary, n the
refractive index, the propagation constant in the direction perpendicular to the
boundary, and k0 the free space wave number. The eigenvalues are = j.
The same equation is solved separately at the elliptical end of the waveguide. The
elliptical port is passive, but the eigenmode is still used in the boundary condition of
the 3D propagating wave simulation. The dimensions of the elliptical end of the
waveguide are such that the frequency range for the lowest propagating mode overlaps
that of the rectangular port.
With the stipulated excitation at the rectangular port and the numerically established
mode shapes as boundary conditions, the following equation is solved for the electric
field vector E inside the waveguide adapter:
1
2
j
( r E ) k 0 r --------- E = 0
where r denotes the relative permeability, j the imaginary unit, the conductivity,
the angular frequency, r the relative permittivity, and 0 the permittivity of free space.
The model uses the following material properties for free space: = 0 and r = r = 1.
2 |
WAVE G U I D E AD AP T E R
Results
Figure 2 shows a single-mode wave propagating through the waveguide.
Figure 2: The x component of the propagating wave inside the waveguide adapter at the
frequency 10 GHz.
Naming the rectangular port Port 1 and the elliptical port Port 2, the S-parameters
describing the reflection and transmission of the wave are defined as follows:
( ( E c E 1 ) E 1 ) dA 1
1
---------------------------------------------------------------S 11 = Port
( E 1 E 1 ) dA 1
Port 1
( E c E 2 ) dA 2
Port 2
S 21 = ----------------------------------------------*
( E 2 E 2 ) dA 2
Port 2
Here Ec is the calculated total field. E1 is the analytical field for the port excitation,
and E2 is the eigenmode calculated from the boundary mode analysis and normalized
3 |
WAVE G U I D E A D AP T E R
with respect to the outgoing power flow. Figure 3 shows the S11 and S21 parameters
as functions of the frequency.
Figure 3: The S11 parameter and S21 parameter (in dB) as a function of the frequency.
This parameter describes the reflections when the waveguide adapter is excited at the
rectangular port and a measure of the part of the wave that is transmitted through the
elliptical port when the waveguide adapter is excited at the rectangular port, respectively.
Modeling Instructions
From the File menu, choose New.
NEW
4 |
WAVE G U I D E AD AP T E R
section.
3 In the Mode analysis frequency text field, type 7[GHz]. The exact value of this
frequency is not important. What matters is that it should be above the cutoff
frequency for the fundamental mode, but below that for the next mode. This setting
ensures that the boundary mode analysis finds the fundamental mode.
Add another boundary mode analysis, for the second port.
section.
3 In the Port name text field, type 2.
4 In the Mode analysis frequency text field, type 7[GHz].
Finally, add the 3D equation for the propagating wave in the waveguide.
Import 1 (imp1)
1 On the Model toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
5 |
WAVE G U I D E A D AP T E R
3 Click Browse.
4 Browse to the models Model Library folder and double-click the file
waveguide_adapter.mphbin.
5 Click Import.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
Air (mat1)
By default, the first material you add applies on all domains so you need not alter any
settings.
E LE C TROM A GN E TI C WAVES , FRE QU E NC Y D OM AI N ( EM W )
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 In the Settings window for Port, locate the Port Properties section.
3 From the Type of port list, choose Numeric.
4 From the Wave excitation at this port list, choose On.
5 Select Boundary 13 only. The wave enters the adapter through the port with a
6 |
WAVE G U I D E AD AP T E R
Port 2
1 On the Physics toolbar, click Boundaries and choose Port.
2 In the Settings window for Port, locate the Port Properties section.
3 From the Type of port list, choose Numeric.
4 In the Port name text field, type 2.
5 Select Boundary 6 only. This is the exit port, the one with an elliptical cross-section.
MESH 1
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
Now set up the study to find the boundary modes and use them when computing the
field distribution over a range of frequencies.
section.
3 In the Search for modes around text field, type 50. This value should be in the vicinity
of the value that you expect the fundamental mode to have. If you do not know this
in advance, you can experiment with some different values or estimate one from
analytical formulas valid for cross-sections resembling yours.
4 From the Transform list, choose Out-of-plane wave number.
section.
3 In the Search for modes around text field, type 50.
4 From the Transform list, choose Out-of-plane wave number.
7 |
WAVE G U I D E A D AP T E R
The default plot shows the norm of the electric field on slices through the
waveguide; you can simplify and improve this plot.
First, delete the Multislice plot.
2 In the Model Builder window, expand the Electric Field (emw) node.
3 Right-click Multislice 1 and choose Delete.
4 Right-click Electric Field (emw) and choose Slice.
5 In the Settings window for Slice, click Replace Expression in the upper-right corner
8 |
WAVE G U I D E AD AP T E R
8 On the 3D plot group toolbar, click Plot. The plot now shows the x-component of
the electric field at the highest frequency, 10 GHz (compare with Figure 2). If you
would like to see the field for other frequencies, you can select them by clicking on
the Electric Field (emw) plot group.
Proceed by checking the plot of the S-parameters as functions of the frequency.
S-Parameter (emw)
1 Select the S-Parameter (emw) plot group under Results in Model Builder. The plot
9 |
WAVE G U I D E A D AP T E R
10 |
WAVE G U I D E A D AP T E R
SMA receptacle
Metal screw
1 |
Model Definition
The Wilkinson power divider is a three-port device composed of 50 and 70.7
microstrip lines on a dielectric substrate with a ground plane and a 100 resistor
mounted between two ports. The model also includes a metal enclosure, screws, and
SMA receptacles connected to each port representing a complete package of a power
divider shown in Figure 1. Model all microstrip lines, the SMA receptacles, screws, and
the metal package using perfect electric conductor (PEC) boundaries. The SMA
receptacle and screw domains enclosed by these PEC boundaries are not part of the
model analysis. The relative dielectric constant, r, of the 60 mil substrate is 3.38. The
boundaries facing the dielectric-filled coaxial connector of the SMA receptacles are
specified as coaxial lumped ports. The 100 resistor is realized via a uniform lumped
port with100 characteristic impedance.
Figure 2: The E-field norm plot shows that the input is evenly split between the two output
ports.
2 |
The S-parameters plotted in Figure 3 show the frequency response of the Wilkinson
power divider. Good input impedance matching characteristics are observed and the
coupled power at each output port is about 3 dB around 3 GHz.
Figure 3: The S-parameters show very good input matching at 3 GHz and evenly divided
power at the two output ports.
References
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
2. R.E. Collin, Foundation of Microwave Engineering, McGraw-Hill, 1992.
Modeling Instructions
From the File menu, choose New.
3 |
NEW
Parameters
1 On the Model toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Name
Expression
Value
Description
r_ring
8.5[mm]
0.0085000 m
w_subs
30[mm]
0.030000 m
Width, substrate
l_subs
26[mm]
0.026000 m
Length, susbtrate
r_inner
0.635[mm]
6.3500E-4 m
r_outer
2.05[mm]
0.0020500 m
l_sma
8[mm]
0.0080000 m
Length, SMA
f_min
2[GHz]
2.0000E9 Hz
Minimum frequency in
sweep
f_max
4[GHz]
4.0000E9 Hz
Maximum frequency in
sweep
lda0
c_const/f_max
0.074948 m
Wavelength, air
h_max
0.2*lda0
0.014990 m
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
GEOMETRY 1
4 |
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type w_subs.
4 In the Depth text field, type l_subs.
5 In the Height text field, type 1.524.
6 Locate the Position section. From the Base list, choose Center.
7 In the z text field, type -0.762.
8 Right-click Component 1 (comp1)>Geometry 1>Block 1 (blk1) and choose Rename.
9 In the Rename Block dialog box, type Substrate in the New label text field.
10 Click OK.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type w_subs.
4 In the Depth text field, type l_subs.
5 In the Height text field, type 20.
6 Locate the Position section. From the Base list, choose Center.
7 In the z text field, type 2.
8 Click the Build Selected button.
5 |
Add a work plane for drawing the layout of the power divider.
13 On the Geometry toolbar, click Work Plane.
Plane Geometry
Add two circles to create the ring strip part.
Circle 1 (c1)
1 On the Work plane toolbar, click Primitives and choose Circle.
2 Click the Zoom Extents button on the Graphics toolbar.
3 In the Settings window for Circle, locate the Size and Shape section.
4 In the Radius text field, type r_ring.
5 Right-click Component 1 (comp1)>Geometry 1>Work Plane 1 (wp1)>Plane
Geometry>Circle 1 (c1) and choose Rename.
6 In the Rename Circle dialog box, type Ring outer in the New label text field.
7 Click OK.
Circle 2 (c2)
1 On the Work plane toolbar, click Primitives and choose Circle.
6 |
2 In the Settings window for Circle, locate the Size and Shape section.
3 In the Radius text field, type r_ring-1.87.
4 Right-click Component 1 (comp1)>Geometry 1>Work Plane 1 (wp1)>Plane
Geometry>Circle 2 (c2) and choose Rename.
5 In the Rename Circle dialog box, type Ring inner in the New label text field.
6 Click OK.
Rectangle 1 (r1)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 2.
4 In the Height text field, type 3.
5 Locate the Position section. In the xw text field, type -1.
6 In the yw text field, type -9.
7 Right-click Component 1 (comp1)>Geometry 1>Work Plane 1 (wp1)>Plane
Geometry>Rectangle 1 (r1) and choose Rename.
8 In the Rename Rectangle dialog box, type Ring cut in the New label text field.
9 Click OK.
Difference 1 (dif1)
1 On the Work plane toolbar, click Booleans and Partitions and choose Difference.
2 Select the object c1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the objects r1 and c2 only.
7 |
Rectangle 2 (r2)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 2.
4 Locate the Position section. In the xw text field, type -1.
5 In the yw text field, type -8.
6 Right-click Component 1 (comp1)>Geometry 1>Work Plane 1 (wp1)>Plane
Geometry>Rectangle 2 (r2) and choose Rename.
7 In the Rename Rectangle dialog box, type Lumped element in the New label text
field.
8 Click OK.
Rectangle 3 (r3)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 3.2.
4 In the Height text field, type 5.
8 |
5 Locate the Position section. From the Base list, choose Center.
6 In the yw text field, type 10.5.
Rectangle 4 (r4)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 3.2.
4 In the Height text field, type 2.
5 Locate the Position section. In the xw text field, type -7.
6 From the Base list, choose Center.
7 In the yw text field, type -12.
Rectangle 5 (r5)
1 On the Work plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 3.2.
4 In the Height text field, type 6.
5 Locate the Position section. In the xw text field, type -8.6.
6 In the yw text field, type -11.
7 Locate the Rotation Angle section. In the Rotation text field, type -28.
Mirror 1 (mir1)
1 On the Work plane toolbar, click Transforms and choose Mirror.
2 Select the objects r4 and r5 only.
3 In the Settings window for Mirror, locate the Input section.
4 Select the Keep input objects check box.
Create a union of all objects except the small rectangle for the resistor (r2) to remove
unnecessary boundaries.
Union 1 (uni1)
1 On the Work plane toolbar, click Booleans and Partitions and choose Union.
2 Select the objects mir1(1), dif1, mir1(2), r3, r4, and r5 only.
3 In the Settings window for Union, locate the Union section.
4 Clear the Keep interior boundaries check box.
9 |
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Radius text field, type r_inner.
4 In the Height text field, type l_sma+2.
5 Locate the Position section. In the x text field, type -7.
6 In the y text field, type -l_subs/2-l_sma.
7 Locate the Axis section. From the Axis type list, choose y-axis.
8 Right-click Component 1 (comp1)>Geometry 1>Cylinder 1 (cyl1) and choose Rename.
9 In the Rename Cylinder dialog box, type Coax inner in the New label text field.
10 Click OK.
Cylinder 2 (cyl2)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
10 |
Block 3 (blk3)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size section.
3 In the Width text field, type 12.7.
4 In the Depth text field, type 12.7.
5 In the Height text field, type 1.65.
6 Locate the Position section. From the Base list, choose Center.
7 In the x text field, type -7.
8 In the y text field, type -(l_subs+1.65)/2.
9 Locate the Axis section. From the Axis type list, choose y-axis.
10 Right-click Component 1 (comp1)>Geometry 1>Block 3 (blk3) and choose Rename.
11 In the Rename Block dialog box, type Flange in the New label text field.
12 Click OK.
Create a union of a couple of objects, the SMA connector, and the flange to remove
unnecessary boundaries.
Union 1 (uni1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Union.
2 Select the objects cyl2 and blk3 only.
3 In the Settings window for Union, locate the Union section.
4 Clear the Keep interior boundaries check box.
Cylinder 3 (cyl3)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
11 |
field.
11 Click OK.
Copy 1 (copy1)
1 On the Geometry toolbar, click Transforms and choose Copy.
2 Select the objects uni1, cyl1, and cyl3 only.
3 In the Settings window for Copy, locate the Displacement section.
4 In the x text field, type 7,14.
Rotate 1 (rot1)
1 On the Geometry toolbar, click Transforms and choose Rotate.
12 |
Cylinder 4 (cyl4)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, locate the Size and Shape section.
3 In the Height text field, type 8.
4 Locate the Position section. In the x text field, type -12.
5 In the y text field, type -10.
6 In the z text field, type -8.
7 Right-click Component 1 (comp1)>Geometry 1>Cylinder 4 (cyl4) and choose Rename.
8 In the Rename Cylinder dialog box, type Screw in the New label text field.
9 Click OK.
Array 1 (arr1)
1 On the Geometry toolbar, click Transforms and choose Array.
2 Select the object cyl4 only.
3 In the Settings window for Array, locate the Size section.
4 In the x size text field, type 2.
5 In the y size text field, type 2.
6 Locate the Displacement section. In the x text field, type 24.
7 In the y text field, type 20.
13 |
DEFINITIONS
Create a set of selections to use when setting up the physics. Begin with the microstrip
line boundaries including the substrate ground plane.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, locate the Input Entities section.
3 From the Geometric entity level list, choose Boundary.
14 |
Explicit 2
1 On the Definitions toolbar, click Explicit.
15 |
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 Select Domains 9, 14, and 19 only.
3 Right-click Component 1 (comp1)>Definitions>Explicit 3 and choose Rename.
4 In the Rename Explicit dialog box, type Coax dielectric in the New label text field.
5 Click OK.
Add a selection for the domains consisting of metal. These domains are not part of the
model analysis.
Explicit 4
1 On the Definitions toolbar, click Explicit.
16 |
2 Select Domains 48, 10, 12, 13, 1618, 20, and 2226 only.
Define the model domain, which is the complement of the metal volume selection.
Complement 1
1 On the Definitions toolbar, click Complement.
2 In the Settings window for Complement, locate the Input Entities section.
3 Under Selections to invert, click Add.
4 In the Add dialog box, select Metal volume in the Selections to invert list.
5 Click OK.
6 Right-click Component 1 (comp1)>Definitions>Complement 1 and choose Rename.
7 In the Rename Complement dialog box, type Model domain in the New label text
field.
8 Click OK.
View 1
Suppress some boundaries to get a view of the interior while setting the physics and
mesh.
17 |
18 |
section.
3 From the Selection list, choose Microstrip line.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 95 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Type of lumped port list, choose Coaxial.
5 From the Wave excitation at this port list, choose On.
19 |
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 51 only.
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Type of lumped port list, choose Coaxial.
Lumped Port 3
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
20 |
3 In the Settings window for Lumped Port, locate the Lumped Port Properties section.
4 From the Type of lumped port list, choose Coaxial.
Lumped Element 1
1 On the Physics toolbar, click Boundaries and choose Lumped Element.
21 |
3 In the Settings window for Lumped Element, locate the Settings section.
4 In the Zelement text field, type 100[ohm].
MATERIALS
Next, assign material properties. First, specify air for all domains.
ADD MATERIAL
1 On the Model toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-In>Air.
4 Click Add to Component in the window toolbar.
MATERIALS
1 On the Model toolbar, click Add Material to close the Add Material window.
Override the material for the substrate domains with a dielectric material of r = 3.38.
Material 2 (mat2)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Substrate.
22 |
4 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
3.38
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
Material 3 (mat3)
1 In the Model Builder window, right-click Materials and choose Blank Material.
2 In the Settings window for Material, locate the Geometric Entity Selection section.
3 From the Selection list, choose Coax dielectric.
4 Locate the Material Contents section. In the table, enter the following settings:
Property
Name
Value
Unit
Property group
Relative permittivity
epsilonr
2.1
Basic
Relative permeability
mur
Basic
Electrical conductivity
sigma
S/m
Basic
23 |
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
STUDY 1
24 |
3 In the Model Builder window, expand the Electric Field (emw) node, then click
Multislice 1.
4 In the Settings window for Multislice, locate the Multiplane Data section.
5 Find the x-planes subsection. In the Planes text field, type 0.
6 Find the y-planes subsection. In the Planes text field, type 0.
7 Find the z-planes subsection. From the Entry method list, choose Coordinates.
8 In the Coordinates text field, type 0.
9 Click to expand the Range section. Select the Manual color range check box.
10 In the Maximum text field, type 1000.
The resulting plot shows the E-field equally split between Port 2 and Port 3.
Compare with Figure 2.
S-Parameter (emw)
The reproduced plot shows the calculated S-parameters. Compare with Figure 3.
25 |
26 |