RF Application Library Manual
RF Application Library Manual
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:
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.
Z0/√2
Isolated Output
Z0 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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
Reference
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
4 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
GEOMETRY 1
1 In the Model Builder window, under Component 1 (comp1) click Geometry 1.
2 In the Settings window for Geometry, locate the Units section.
3 From the Length unit list, choose mm.
Rectangle 1 (r1)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape 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 Right-click Rectangle 1 (r1) and choose Build Selected.
5 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
Rectangle 2 (r2)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape 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 and Shape 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 Right-click Rectangle 3 (r3) and choose Build Selected.
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Block, locate the Size and Shape section.
3 In the Width text field, type l_s.
4 In the Depth text field, type l_s.
5 In the Height text field, type thickness.
6 Locate the Position section. From the Base list, choose Center.
7 In the z text field, type thickness/2.
8 Click the Build All Objects button.
Union 1 (uni1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Union.
2 In the Settings window for Union, type Substrate in the Label text field.
3 Locate the Union section. 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 Resulting objects
selection check box.
6 From the Show in physics list, choose All levels.
7 Click the Build All Objects button.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Package in the Label text field.
3 Locate the Size and Shape section. 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 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
DEFINITIONS
Create a selection for the microstrip lines.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Microstrip Line in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundary 13 only.
View 1
Hide three boundaries to get a better view of the interior parts when reviewing the
mesh.
9 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
View 1
1 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
2 From the Geometric entity level list, choose Boundary.
3 Select Boundaries 1, 2, and 4 only.
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.
10 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
Lumped Port 3
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 15 only.
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 Home 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.
11 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
MATERIALS
1 On the Home toolbar, click Add Material to close the Add Material window.
Create a dielectric material of εr = 3.38 overriding air in the substrate.
Material 2 (mat2)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Substrate in the Label text field.
3 Locate the Geometric Entity Selection section. From the Selection list, choose
Substrate.
4 Locate the Material Contents section. In the table, enter the following settings:
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Build All.
12 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
STUDY 1
RESULTS
1 In the Settings window for 3D Plot Group, locate the Data section.
2 From the Parameter value (freq (Hz)) list, choose 3E9.
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 thickness/2.
13 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
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.
Compare the resulting plot with that shown in Figure 3.
14 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
1D Plot Group 4
1 On the Home 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.
Global 1
On the 1D Plot Group 4 toolbar, click Global.
1D Plot Group 4
1 In the Settings window for Global, locate the y-Axis Data section.
2 In the table, enter the following settings:
15 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
Global Evaluation 1
On the Results toolbar, click Global Evaluation.
Derived Values
1 In the Settings window for Global Evaluation, locate the Data section.
2 From the Parameter selection (freq) list, choose From list.
3 In the Parameter values (freq (Hz)) list, select 3E9.
4 In the Parameter values (freq (Hz)) list, select 3E9.
5 Locate the Expression section. In the Expression text field, type
arg(emw.S21)-arg(emw.S31).
16 | BRANCH-LINE COUPLER
Solved with COMSOL Multiphysics 5.2
Model Definition
This example considers three geometries:
The cavity walls are assumed to be a good conductor, such as copper, with an electric
conductivity of 5.7·107 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 results 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 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.
Another way to calculate Q-factor at the dominant mode is via equations in Ref. 1. For
a rectangular cavity, the dominant mode is TE101, at which the cavity provides the
lowest resonant frequency. The Q-factor and resonant frequency at this mode is
π 2 π 2
Q TE101 = -----------------------------, f TE101 = ------------------ --- + ---
1.1107η 1
2π με a c
R s 1 + -------
a
2b
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
2.40492 2
Q TM 010 = --------------------------, f TM 010 = ------------------ ---------------------
1.2025η 1
2π με a
R s 1 + ---
a
h
1.0041η 2.744
Q TM 011 = ----------------------, f TM 011 = ---------------------
Rs 2πa με
ωr μ
Rs = ----------
2σ
These two analytical approaches are compared with the Q-factor obtained from
Eigenfrequency analysis.
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
Figure 1: Arrow plots of electric and magnetic fields. Slice plot of electric field.
Reference
1. C.A. Balanis, Advanced Engineering Electromagnetics, John Wiley & Sons, 1989.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Eigenfrequency.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 Click Load from File.
4 Browse to the application’s Application Library folder and double-click the file
cavity_resonators_parameters.txt.
Since air and lossy wall materials will be used on multiple components, add them on
the global material node. They will be linked to each individual component later on.
5 On the Home toolbar, click Windows and choose Add Material.
ADD MATERIAL
1 Go to the Add Material window.
2 In the tree, select Built-In>Air.
3 Click Add to Global Materials.
GLOBAL DEFINITIONS
Air (mat1)
On the Home toolbar, click Add Material to close the Add Material window.
Material 2 (mat2)
1 In the Model Builder window, under Global Definitions right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Lossy Wall in the Label text field.
3 Click to expand the Material properties section. Locate the Material Properties
section. In the Material properties tree, select Basic Properties>Relative Permittivity.
4 Click Add to Material.
5 In the Material properties tree, select Basic Properties>Relative Permeability.
GEOMETRY 1
Create a block for the rectangular cavity.
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size and Shape 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.
6 Click the Build All Objects button.
MATERIALS
Assign material properties on the model by linking the global material already created.
First, apply air to all domains.
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Material Link.
DEFINITIONS
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.
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 application’s Application 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
choose Free Tetrahedral.
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 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.
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 Select the Desired number of eigenfrequencies check box.
4 In the associated text field, type 1.
5 Select the Search for eigenfrequencies around check box.
6 In the associated 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:
RESULTS
1 In the Model Builder window, under Results right-click Electric Field (emw) and
choose Arrow Volume.
2 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Electric>emw.Ex,emw.Ey,emw.Ez -
Electric field.
3 On the Electric Field (emw) toolbar, click Plot.
4 In the Model Builder window, right-click Electric Field (emw) and choose Arrow
Volume.
5 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Magnetic>emw.Hx,emw.Hy,emw.Hz -
Magnetic field.
6 Locate the Arrow Positioning section. Find the z grid points subsection. In the Points
text field, type 1.
7 Locate the Coloring and Style section. From the Color list, choose White.
8 On the Electric Field (emw) toolbar, click Plot.
9 Click the Zoom Extents button on the Graphics toolbar.
Compare the resulting plot 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 Home toolbar, click Add Plot Group and choose 3D Plot Group.
2 In the Model Builder window, right-click 3D Plot Group 2 and choose Surface.
3 In the Settings window for Surface, click Replace Expression in the upper-right corner
of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Heating and losses>emw.Qsrh - Surface
resistive losses.
4 Locate the Coloring and Style section. From the Color table list, choose
ThermalEquidistant.
5 In the Model Builder window, click 3D Plot Group 2.
6 In the Settings window for 3D Plot Group, type Surface Losses (emw) in the
Label text field.
ROOT
Next, set up a model for the cylindrical cavity.
ADD PHYSICS
1 On the Home 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 Home toolbar, click Add Physics to close the Add Physics window.
ADD STUDY
1 On the Home 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 Empty Study.
You will copy the settings from the existing study later on.
4 Click Add Study in the window toolbar.
5 On the Home 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.
MATERIALS
Assign material properties on the second model. Apply air to all domains.
1 In the Model Builder window, under Component 2 (comp2) right-click Materials and
choose Material Link.
DEFINITIONS
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.
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 2
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 application’s Application 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
choose Free Tetrahedral.
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 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 1
1 Select both Study 1> Step 1: Eigenfrequency 1 and Study 1> Parametric Sweep 1 using
shift-key. Copy them and paste on Study 2.
STUDY 2
1 In the Model Builder window, right-click Study 2 and choose Paste Multiple Items.
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:
RESULTS
4 In the Model Builder window, right-click Electric Field (emw2) and choose Arrow
Volume.
5 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
2>Electromagnetic Waves, Frequency Domain
2>Magnetic>emw2.Hx,emw2.Hy,emw2.Hz - Magnetic field.
6 Locate the Arrow Positioning section. Find the z grid points subsection. In the Points
text field, type 1.
7 Locate the Coloring and Style section. From the Color list, choose White.
8 On the Electric Field (emw2) toolbar, click Plot.
9 Click the Zoom Extents button on the Graphics toolbar.
The plot should now look like that in Figure 1, middle.
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 Home 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) (sol8).
4 Right-click 3D Plot Group 4 and choose Surface.
5 In the Settings window for Surface, click Replace Expression in the upper-right corner
of the Expression section. From the menu, choose Model>Component
2>Electromagnetic Waves, Frequency Domain 2>Heating and losses>emw2.Qsh -
Surface losses.
6 Locate the Coloring and Style section. From the Color table list, choose
ThermalEquidistant.
7 In the Model Builder window, click 3D Plot Group 4.
8 In the Settings window for 3D Plot Group, type Surface Losses (emw2) in the
Label text field.
3 Locate the Coloring and Style section. From the Color list, choose Blue.
4 On the Surface Losses (emw2) toolbar, click Plot.
ROOT
Now add a model for the spherical cavity.
ADD PHYSICS
1 On the Home 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 Home toolbar, click Add Physics to close the Add Physics window.
ADD STUDY
1 On the Home 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 Empty Study.
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.
4 Click the Build All Objects button.
MATERIALS
Assign material properties on the third model. Apply air to all domains.
1 In the Model Builder window, under Component 3 (comp3) right-click Materials and
choose Material Link.
DEFINITIONS
Add variables and two integration coupling operators.
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.
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 3
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 application’s Application 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
choose Free Tetrahedral.
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 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 2
1 Select both Study 2> Step 1: Eigenfrequency 1, and Study 2>Parametric Sweep 1 using
shift-key. Copy them and paste on Study 3.
STUDY 3
1 In the Model Builder window, right-click Study 3 and choose Paste Multiple Items.
Step 1: Eigenfrequency 1
1 In the Model Builder window, under Study 3 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:
RESULTS
Again, add a surface plot of the surface losses and an arrow plot of the surface current
(Figure 2, bottom).
3D Plot Group 6
1 On the Home 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 3/Parametric Solutions 3 (12) (sol14).
4 Right-click 3D Plot Group 6 and choose Surface.
5 In the Settings window for Surface, click Replace Expression in the upper-right corner
of the Expression section. From the menu, choose Model>Component
3>Electromagnetic Waves, Frequency Domain 3>Heating and losses>emw3.Qsh -
Surface losses.
6 Locate the Coloring and Style section. From the Color table list, choose
ThermalEquidistant.
7 In the Model Builder window, click 3D Plot Group 6.
8 In the Settings window for 3D Plot Group, type Surface Losses (emw3) in the
Label text field.
Derived Values
Finish by evaluating the Q-factor and resonant frequency. Compare them with those
values in Table 1, Table 2 and Table 3.
Global Evaluation 1
On the Results toolbar, click Global Evaluation.
Derived Values
1 In the Settings window for Global Evaluation, locate the Data section.
2 From the Data set list, choose Study 1/Parametric Solutions 1 (sol2).
3 From the Eigenfrequency selection list, choose First.
25 Locate the Data section. From the Data set list, choose Study 2/Parametric Solutions
2 (6) (sol8).
26 Click the Evaluate button.
27 From the Data set list, choose Study 3/Parametric Solutions 3 (12) (sol14).
28 Click the Evaluate button.
Connecting a 3D Electromagnetic
Wave Model to an Electrical Circuit
Introduction
An application 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.
1 2 3
Zsource Zload
Model Definition
The geometry in this example 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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 In the Select physics tree, select AC/DC>Electrical Circuit (cir).
5 Click Add.
6 Click Study.
7 In the Select study tree, select Preset Studies for Selected Physics Interfaces>Frequency
Domain.
8 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home 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
Solved with COMSOL Multiphysics 5.2
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
Create the geometry of the coaxial cable using two cylinders.
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type Coax outer in the Label text field.
3 Locate the Size and Shape section. 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.
Cylinder 2 (cyl2)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type Coax inner in the Label text field.
3 Locate the Size and Shape section. 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.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object cyl1 only.
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
Solved with COMSOL Multiphysics 5.2
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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.
Voltage Source 1
1 On the Electrical Circuit 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:
4 Locate the Device Parameters section. From the Source type list, choose AC-source.
Resistor 1
1 On the Electrical Circuit 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
Solved with COMSOL Multiphysics 5.2
4 Locate the Device Parameters section. In the R text field, type Z_coax.
Resistor 2
1 On the Electrical Circuit toolbar, click Resistor.
2 In the Settings window for Resistor, locate the Node Connections section.
3 In the table, enter the following settings:
4 Locate the Device Parameters section. In the R text field, type Z_coax.
External I Vs. U 1
1 On the Electrical Circuit 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:
4 Locate the External Device section. From the V list, choose Lumped port voltage
(emw/lport1).
External I Vs. U 2
1 On the Electrical Circuit 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:
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
Solved with COMSOL Multiphysics 5.2
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 Home 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 Home toolbar, click Add Material to close the Add Material window.
MESH 1
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 Element size list, choose Fine.
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
Solved with COMSOL Multiphysics 5.2
STUDY 1
Solution 1 (sol1)
1 On the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 1 (sol1) 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
Solved with COMSOL Multiphysics 5.2
1 In the Model Builder window, under Results right-click Electric Field (emw) and
choose Arrow Volume.
2 Right-click Electric Field (emw) and choose Arrow Volume.
3 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Magnetic>emw.Hx,emw.Hy,emw.Hz -
Magnetic field.
4 Locate the Coloring and Style section. From the Color list, choose Green.
5 Right-click Electric Field (emw) and choose Arrow Volume.
6 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Energy and
power>emw.Poavx,...,emw.Poavz - Power flow, time average.
7 Locate the Coloring and Style section. From the Color list, choose Blue.
Compare the resulting plot with that shown in Figure 2.
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
Solved with COMSOL Multiphysics 5.2
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 mode—with the electric and magnetic fields normal
to the direction of propagation along the cable—modeling a 2D cross section suffices
to compute the fields and the impedance. For this example, assume perfect conductors
and a lossless dielectric with relative permittivity εr = 2.4. The inner and outer radii are
0.5 mm and 3.43 mm, respectively.
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
The voltage and current in the direction out of the plane are positive for integration
paths oriented as in Figure 2.
2
- . C
1 2
1
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
1 μ0 ro
- log ----- ≈ 74.5 Ω
Z 0,analytic = ------ ---------
2π ε r ε ri
0
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 approaches the analytic value of 74.5 Ω.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Mode Analysis.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
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:
MATERIALS
Define a dielectric material for the region between the conductors.
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Insulator in the Label text field.
3 Locate the Material Contents section. In the table, enter the following settings:
DEFINITIONS
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.
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.
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:
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 = Vi − Vo (see Equation 1) and
to have a positive current value correspond to a current in the positive z direction.
The emw. prefix gives the correct physics-interface scope for the electric and
magnetic field vector components.
MESH 1
Use the default mesh.
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Build All.
STUDY 1
RESULTS
1 In the Model Builder window, under Results right-click Electric Field (emw) and
choose Arrow Surface.
2 In the Settings window for Arrow Surface, locate the Arrow Positioning section.
3 Find the x grid points subsection. In the Points text field, type 21.
4 Find the y grid points subsection. In the Points text field, type 21.
5 Locate the Coloring and Style section. Select the Scale factor check box.
6 Use the slider to adjust the arrow lengths.
7 From the Color list, choose White.
8 On the Electric Field (emw) toolbar, click Plot.
9 Click the Zoom Extents button on the Graphics toolbar.
Compare the resulting plot with that shown in Figure 2.
To find out the integration contour orientations, plot the tangent vector, t1, along
the boundaries as follows:
10 In the Model Builder window, right-click Electric Field (emw) and choose Arrow Line.
11 In the Settings window for Arrow Line, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Geometry>tx,ty - Tangent.
12 Locate the Coloring and Style section. In the Number of arrows text field, type 50.
A comparison with Equation 1 reveals that the line integral for the voltage computes
the potential difference Vi − Vo. 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
the result plot.
Global Evaluation 1
On the Results toolbar, click Global Evaluation.
Derived Values
1 In the Settings window for Global Evaluation, click Replace Expression in the
upper-right corner of the Expression section. From the menu, choose
Model>Component 1>Definitions>Variables>Z0_model - Characteristic impedance.
2 Click the Evaluate button.
TA BL E
1 Go to the Table window.
The result, roughly 74.65 Ω, is within 0.2% of the analytic value, 74.53 Ω.
This example considers a pulse propagating down a coaxial transmission line for three
different termination types: short, open, and matched. The signal propagation time is
deduced from the reflected waves detected at the input port.
Model Definition
The model setup, schematically shown 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
At one end of the coaxial cable, or coax for short, a lumped port boundary condition
excites the structure; specify a transient excitation pulse, V0(t), by using a Gaussian
pulse-windowed sine function. Apply the excitation as a current of magnitude
I(t) = V 0(t) / Zref flowing tangentially to the excitation boundary. Here Zref refers to
the specified characteristic impedance between the voltage generator and the model.
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
Solved with COMSOL Multiphysics 5.2
At the other end of the coax, consider, in turn, three different boundary conditions:
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
Solved with COMSOL Multiphysics 5.2
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
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D Axisymmetric.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Transient
(temw).
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
Solved with COMSOL Multiphysics 5.2
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Time Dependent.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Next, define the excitation, V0(t), in terms of a Gaussian pulse and a sine function.
Analytic 1 (an1)
1 On the Home toolbar, click Functions and choose Global>Analytic.
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
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Analytic, type V0 in the Function name text field.
3 Locate the Definition section. In the Expression text field, type
gauss_pulse(t)*sin(2*pi*f*t).
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 and Shape section.
3 In the Width text field, type R_coax-r_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
Solved with COMSOL Multiphysics 5.2
DEFINITIONS
Set up a point probe for plotting the electric field component Er while solving. You
will also use this plot to reproduce Figure 2.
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.
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
Solved with COMSOL Multiphysics 5.2
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 3 only.
3 In the Settings window for Lumped Port, locate the Settings section.
4 In the Zref text field, type (Z0_const/2/pi)*log(R_coax/r_coax).
ADD MATERIAL
1 On the Home 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 Home 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
choose Free Triangular.
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 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
Solved with COMSOL Multiphysics 5.2
STUDY 1
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
Solved with COMSOL Multiphysics 5.2
STUDY 1
On the Home toolbar, click Compute.
RESULTS
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.
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
Solved with COMSOL Multiphysics 5.2
DEFINITIONS
1 In the Model Builder window, under Component 1 (comp1)>Definitions>Domain Point
Probe 1 click Point Probe Expression 1 (ppb1).
2 In the Settings window for Point Probe Expression, click to expand the Table and
window settings section.
3 Locate the Table and Window Settings section. From the Output table list, choose New
table.
With these settings you get a plot for the short and open termination cases in the
same plot window.
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 )
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
Solved with COMSOL Multiphysics 5.2
STUDY 1
1 In the Model Builder window, click Study 1.
2 In the Settings window for Study, locate the Study Settings section.
3 Clear the Generate default plots check box.
4 On the Home toolbar, click Compute.
RESULTS
DEFINITIONS
1 In the Model Builder window, under Component 1 (comp1)>Definitions>Domain Point
Probe 1 click Point Probe Expression 1 (ppb1).
2 In the Settings window for Point Probe Expression, click to expand the Table and
window settings section.
3 Locate the Table and Window Settings section. From the Output table list, choose New
table.
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
Solved with COMSOL Multiphysics 5.2
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
On the Home toolbar, click Compute.
RESULTS
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
DOMAIN EQUATIONS
An electromagnetic wave propagating in a coaxial cable is characterized by transverse
electromagnetic (TEM) fields. Assuming time-harmonic fields with complex
amplitudes containing the phase information, you have:
C
E = e r ---- e j ( ωt – kz )
r
C
H = e ϕ ------- e j ( ωt – kz )
rZ
2π
k = ------
λ
In the air, the electric field also has a finite axial component whereas the magnetic field
is purely azimuthal. Thus it is possible to model the antenna using an axisymmetric
transverse magnetic (TM) formulation, and the wave equation becomes scalar in H ϕ :
BOUNDARY CONDITIONS
The boundary conditions for the metallic surfaces are:
n×E = 0
At the feed point, a matched coaxial port boundary condition is used to make the
boundary transparent to the wave. The antenna is radiating into free space, but you
can only discretize a finite region. Therefore, truncate the geometry some distance
from the antenna using a scattering boundary condition allowing for outgoing
spherical waves to pass with very little reflections. A symmetry boundary condition is
automatically applied on boundaries at r = 0.
2 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
obtain maximum transmission into free space. This is quite well fulfilled in the high
frequency range.
Figure 3 shows the antenna radiation pattern in the near-field for three different
frequencies. The effect of the finite diameter of the ground plane is to lift the main lobe
from the horizontal plane. For an infinite ground plane or in the high-frequency limit,
the radiation pattern is symmetric around zero elevation. This is easy to understand, as
an infinite ground plane can be replaced by a mirror image of the monocone below the
plane. Such a biconical antenna is symmetric around zero elevation and has its main
lobe in the horizontal direction. The decreased lobe lifting at higher frequencies is just
about visible in Figure 3.
Figure 4 shows the antenna radiation pattern in the far field for the same frequencies
as the radiation pattern at the boundary in Figure 3.
3 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D Axisymmetric.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GEOMETRY 1
Import 1 (imp1)
1 On the Home toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 Click Browse.
4 Browse to the application’s Application Library folder and double-click the file
conical_antenna.mphbin.
5 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
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.
GLOBAL DEFINITIONS
Prepare for the impedance computation by making a few definitions.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
6 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
DEFINITIONS
Variables 1
1 On the Home 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:
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 Settings window for Explicit, type Air in the Label text field.
3 Select Domain 1 only.
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Dielectric in the Label text field.
The dielectric domain is inside the coaxial cable just below the cone. It is easier to
select it if you zoom in a little.
3 Select Domain 2 only.
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Outer Air Boundaries in the Label text
field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundaries 14 and 15 only.
With all selections and expressions now defined, it is time to set up the materials and
the physics of the model.
7 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
ADD MATERIAL
1 On the Home 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 Home 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, under Component 1 (comp1) 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:
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.
4 From the Type of port list, choose Coaxial.
5 From the Wave excitation at this port list, choose On.
8 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
MESH 1
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Free Triangular.
Because this is a relatively small 2D model, you could easily use one of the finer mesh
defaults and get a good accuracy. The following settings however are designed to
get you a fine mesh mostly where it is needed.
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
Solved with COMSOL Multiphysics 5.2
STUDY 1
11 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
RESULTS
1 In the Model Builder window, expand the Electric Field (emw) node, then click Surface
1.
2 In the Settings window for Surface, click Replace Expression in the upper-right corner
of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Magnetic>Magnetic field>emw.Hphi -
Magnetic field, phi component.
3 Click to expand the Range section. Select the Manual color range check box.
4 In the Minimum text field, type -0.5.
5 In the Maximum text field, type 0.5.
6 Locate the Coloring and Style section. From the Color table list, choose Wave.
7 On the Electric Field (emw) toolbar, click Plot.
8 Click the Zoom Extents button on the Graphics toolbar.
12 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
1D Plot Group 6
On the Home toolbar, click Add Plot Group and choose 1D Plot Group.
Global 1
On the 1D Plot Group 6 toolbar, click Global.
1D Plot Group 6
1 In the Settings window for Global, locate the y-Axis Data section.
2 In the table, enter the following settings:
3 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.
13 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
Line Graph 1
On the Polar Plot Group 7 toolbar, click Line Graph.
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.
5 Select the Description check box.
6 In the associated text field, type Near-field radiation pattern.
7 Locate the θ Angle Data section. From the Parameter list, choose Expression.
8 In the Expression text field, type atan2(z,r).
9 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.
10 Click to collapse the Coloring and style section. Click to expand the Legends section.
Select the Show legends check box.
11 From the Legends list, choose Manual.
14 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
Legends
200 MHz
800 MHz
1.5 GHz
13 Click to collapse the Legends section. On the Polar Plot Group 7 toolbar, click Plot.
14 Click the Zoom Extents button on the Graphics toolbar.
Your near-field radiation plot should look like that in Figure 3.
15 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
16 | CONICAL ANTENNA
Solved with COMSOL Multiphysics 5.2
Polytetrafluoroethylene
(PTFE) ring
Cylindrical cavity
Lumped port fed on a
Dielectric
microstrip line
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 lπ 2
f nml = ---------------------- ------------ + -----
c
2π ε r μ r a d
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 example 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.
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.
Reference
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
GEOMETRY 1
1 In the Model Builder window, under Component 1 (comp1) click Geometry 1.
2 In the Settings window for Geometry, locate the Units section.
3 From the Length unit list, choose mm.
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type Cavity in the Label text field.
3 Locate the Size and Shape section. In the Radius text field, type 25.
4 In the Height text field, type 100.
5 Right-click Cavity and choose Build Selected.
Rectangle 1 (r1)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type l_slot.
Create a substrate.
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Bottom_plate in the Label text field.
3 Locate the Size and Shape section. 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 Right-click Bottom_plate and choose Build Selected.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Bottom_feed in the Label text field.
3 Locate the Size and Shape section. 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 Right-click Bottom_feed and choose Build Selected.
Block 3 (blk3)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Housing in the Label text field.
3 Locate the Size and Shape section. In the Width text field, type 60.
4 In the Depth text field, type 60.
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 Right-click Rotate 1 (rot1) and choose Build Selected.
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.
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 Home 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 Home toolbar, click Add Material to close the Add Material window.
Create a substrate material.
Material 2 (mat2)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Substrate in the Label text field.
3 Select Domains 2, 3, 7, and 8 only.
4 Locate the Material Contents section. In the table, enter the following settings:
Material 3 (mat3)
1 Right-click Materials and choose Blank Material.
2 In the Settings window for Material, type PTFE in the Label text field.
3 Select Domain 6 only.
4 Locate the Material Contents section. In the table, enter the following settings:
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Build All.
STUDY 1
RESULTS
1 In the Settings window for 3D Plot Group, locate the Data section.
2 From the Parameter value (freq (Hz)) list, choose 3.53E9.
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 z-planes subsection. In the Planes text field, type 0.
7 On the Electric Field (emw) toolbar, click Plot.
8 Click the Zoom Extents button on the Graphics toolbar.
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.
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
Perfectly matched layer (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 antenna’s 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
approaches the analytic solution.
1 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
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.
The polar plot in Figure 3 of the far-field pattern in the xy-plane shows the expected
isotropic radiation pattern. The 3D visualization of the far-field intensity in Figure 4
shows the expected torus-shaped pattern.
The impedance as seen by the port is evaluated to be 121 + 28i Ω, which agrees
reasonably with expectations. In the limit as the antenna radius and gap height go to
zero and in the limit of mesh refinement, the model approaches the analytic solution
for a dipole antenna.
2 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
4 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
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:
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_antenna.
5 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
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 Click the Build All Objects button.
6 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
DEFINITIONS
Create a set of selections to be used when setting up the physics. First, create a selection
for the antenna radiator surface.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Antenna in the Label text field.
3 Locate the Input Entities section. 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.
6 Click OK.
7 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
View 1
Hide some domains and boundaries. This helps to see the interior parts when setting
up the physics and reviewing the mesh.
8 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
View 1
1 Select Domains 1 and 2 only.
View 1
1 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
2 From the Geometric entity level list, choose Boundary.
9 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
10 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
2 Click the Zoom In button on the Graphics toolbar a couple of times to see the small
gap between antenna radiators clearly.
3 Select Boundaries 16, 17, 29, and 40 only.
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 x
0 y
1 z
Far-Field Domain 1
On the Physics toolbar, click Domains and choose Far-Field Domain.
Far-Field Calculation 1
1 In the Model Builder window, expand the Far-Field Domain 1 node, then click Far-Field
Calculation 1.
11 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
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 Home 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
1 Go to the Add Material window.
2 In the tree, select Built-In>Copper.
3 Click Add to Component in the window toolbar.
MATERIALS
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 Home toolbar, click Add Material to close the Add Material window.
12 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Build All.
STUDY 1
RESULTS
1 In the Model Builder window, expand the Electric Field (emw) node, then click
Multislice 1.
2 In the Settings window for Multislice, locate the Multiplane Data section.
13 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
1 In the Model Builder window, under Results click 2D Far Field (emw).
2 In the Settings window for Polar Plot Group, click to expand the Axis section.
3 Select the Manual axis limits check box.
4 In the r minimum text field, type 0.
5 In the r maximum text field, type 1.
6 On the 2D Far Field (emw) toolbar, click Plot.
The plotted H-plane pattern is omni-directional (isotropic) on the xy-plane as
shown in Figure 3. The E- and H-plane of a linearly polarized antenna are defined
by the antenna main polarization. The E-plane includes the main polarization that
is Ez in this model while the H-plane is perpendicular to the main polarization.
Derived Values
Finish the result analysis by evaluating the port impedance.
Global Evaluation 2
On the Results toolbar, click Global Evaluation.
Derived Values
1 In the Settings window for Global Evaluation, click Replace Expression in the
upper-right corner of the Expression section. From the menu, choose
Model>Component 1>Electromagnetic Waves, Frequency Domain>Ports>emw.Zport_1
- Lumped port impedance.
2 Click the Evaluate button.
14 | DIPOLE ANTENNA
Solved with COMSOL Multiphysics 5.2
This particular model is intended as a tutorial that shows how to model the full
time-dependent wave equation in dispersive media, such as plasmas and
semiconductors (and any linear medium describable by a sum of Drude-Lorentz
resonant terms). The dispersion of the medium in the frequency domain is assumed to
be of the form
2
ωp
ε r ( ω ) = ε∞ – -------------------------------------
2 2
-, (1)
ω – jΓ i ω – ω i
In this model, the wave equation for the magnetic vector potential
–1 ∂A ∂D
∇ × μ r ( ∇ × A ) + μ 0 σ ------- = μ 0 -------- , (2)
∂t ∂t
D = ε0 ε∞ E + P , (3)
is solved together with an ordinary differential equation for the polarization field,
obtained by a Fourier transformation of Equation 1,
∂2 ∂ 2 2
-------2- + Γ i ----- + ω i P = ε 0 fω p E . (4)
∂t ∂t
Notice that this model is not primarily intended to demonstrate the anomalously high
transmission through hole arrays, but rather to demonstrate temporal dispersion
modeling.
Model Definition
The geometry consists of a single dispersive slab of thickness 1 μm with a slit of width
0.5 μm in it. The wavelength used is 1 μm. Periodic boundary conditions are applied
to make the structure physically appear as an array of slits. The source of
electromagnetic radiation is a plane wave pulse with flat front and Gaussian temporal
shape.
Figure 1: The y-component of the electric field at the input boundary. The left part shows
the incident pulse and the right part shows the reflected pulse.
Figure 2 shows the probe plot of the y-component of the polarization at a point in the
entrance of the slit. Notice the propagation delay between the incident field, shown in
Figure 2: The y-component of the polarization at a point at the entrance of the slit.
Figure 3 shows the probe plot of the y-component of the polarization field at a point
at the rear end of the slit.
Figure 3: The y-component of the polarization at a point at the exit of the slit.
Figure 4 shows a field plot of the y-component of the polarization field after the last
time step (100 fs).
Finally, the out-of-plane component of the magnetic field and, as an overlaid contour
plot, the y-component of the polarization field are shown in Figure 5, after 100 fs.
Figure 5: The out-of-plane component of the magnetic field and the y-component of the
polarization field (contours) after 100 fs.
Reference
1. T. W. Ebbesen H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wolff, “Extraordinary
Optical Transmission Through Sub-wavelength Hole Arrays,” Nature, vol. 391, pp.
667-9, 1998.
Model Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Transient
(temw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Time Dependent.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
Add some parameters that will define the geometry and the properties of the
incident field.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
DEFINITIONS
Variables 1
1 In the Model Builder window, under Component 1 (comp1) right-click Definitions and
choose Variables.
Now add some variables that defines the incident field and the material properties.
2 In the Settings window for Variables, locate the Variables section.
GEOMETRY 1
The geometry is simple, consisting of only three centered rectangles.
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type lambda0.
4 In the Height text field, type 6*lambda0.
5 Locate the Position section. From the Base list, choose Center.
Rectangle 2 (r2)
1 In the Model Builder window, right-click Rectangle 1 (r1) and choose Duplicate.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type 20*lambda0.
Rectangle 3 (r3)
1 In the Model Builder window, under Component 1 (comp1)>Geometry 1 right-click
Rectangle 1 (r1) and choose Duplicate.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Height text field, type 0.5*lambda0.
4 Click the Build All Objects button.
DEFINITIONS
Now, add three integration operator that will be used for probing the field and the
polarization in three different points.
Integration 1 (intop1)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, locate the Source Selection section.
3 From the Geometric entity level list, choose Point.
Integration 2 (intop2)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, locate the Source Selection section.
3 From the Geometric entity level list, choose Point.
Integration 3 (intop3)
1 On the Definitions toolbar, click Component Couplings and choose Integration.
2 In the Settings window for Integration, locate the Source Selection section.
3 From the Geometric entity level list, choose Point.
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 )
1 In the Model Builder window, under Component 1 (comp1) click Electromagnetic
Waves, Transient (temw).
2 In the Settings window for Electromagnetic Waves, Transient, locate the Components
section.
3 From the Electric field components solved for list, choose In-plane vector to solve only
for the in-plane components of the field.
4 From the ε∞ list, choose User defined. From the list, choose Diagonal.
5 In the ε∞ table, enter the following settings:
4 0 0
0 4 0
0 0 1
Drude-Lorentz Polarization 1
1 On the Physics toolbar, click Attributes and choose Drude-Lorentz Polarization.
2 In the Settings window for Drude-Lorentz Polarization, locate the Drude-Lorentz
Dispersion Model section.
3 In the fn text field, type 1.
4 In the ωn text field, type omega_1.
5 In the Γn text field, type gamma_1.
3 In the Settings window for Wave Equation, Electric, locate the Electric Displacement
Field section.
4 From the εr list, choose User defined. Locate the Magnetic Field section. From the μr
list, choose User defined. Locate the Conduction Current section. From the σ list,
choose User defined.
Use scattering boundary conditions to excite the wave and to absorb it.
0 x
E_pulse*E_bnd y
0 z
Periodic Condition 1
1 On the Physics toolbar, click Boundaries and choose Periodic Condition.
MESH 1
Since a periodic boundary condition is used, the mesh should also be the same on the
top and bottom edge. Thus, add first an edge mesh and copy the mesh points to the
opposite edge. Then add a triangular mesh.
Edge 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose More Operations>Edge.
Copy Edge 1
1 In the Model Builder window, 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.
Copy Edge 2
1 Right-click Mesh 1 and choose More Operations>Copy Edge.
2 Select Boundary 10 only.
3 In the Settings window for Copy Edge, locate the Destination Boundaries section.
4 Select the Active toggle button.
Copy Edge 3
1 Right-click Mesh 1 and choose More Operations>Copy Edge.
2 Select Boundary 15 only.
3 In the Settings window for Copy Edge, locate the Destination Boundaries section.
4 Select the Active toggle button.
Size
1 Right-click Mesh 1 and choose Free Triangular.
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 lambda0/6.
STUDY 1
Solution 1 (sol1)
1 On the Study toolbar, click Show Default Solver, to be able to make some
modifications of the solver settings.
Force the solver to use a fixed small step size that resolves the temporal field
oscillations.
2 In the Model Builder window, expand the Solution 1 (sol1) node, then click
Time-Dependent Solver 1.
3 In the Settings window for Time-Dependent Solver, click to expand the Time
stepping section.
4 Locate the Time Stepping section. From the Steps taken by solver list, choose Manual.
DEFINITIONS
Before computing the solution, define the three Global Variable Probes that can be
used for monitoring the computation progress.
Study 1
On the Home toolbar, click Compute.
RESULTS
Electric Field
Modify this surface plot to show the y-component of the Drude-Lorentz polarization.
1 In the Model Builder window, expand the Electric Field node, then click Surface 1.1.
2 In the Settings window for Surface, locate the Expression section.
3 In the Expression text field, type temw.Poscy.
4 On the Electric Field toolbar, click Plot.
The first probe plot should look like the one above.
The second probe plot should look like the one above.
Finally, the third probe plot should look like the one above.
Now, add a surface plot of the z-component of the magnetic field and overlay a
contour plot of the y-component of the Drude-Lorentz polarization.
2D Plot Group 5
1 On the Home toolbar, click Add Plot Group and choose 2D Plot Group.
2 In the Model Builder window, right-click 2D Plot Group 5 and choose Surface.
3 In the Settings window for Surface, locate the Expression section.
4 In the Expression text field, type temw.Hz.
5 Locate the Coloring and Style section. From the Color table list, choose Cyclic.
6 Right-click 2D Plot Group 5 and choose Contour.
7 In the Settings window for Contour, locate the Expression section.
8 In the Expression text field, type temw.Poscy.
9 On the 2D Plot Group 5 toolbar, click Plot.
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-polarized—that is, if the electric field vector is in the same plane as
the Poynting vector and the surface normal—then there are 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.
1 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
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
2n 1 cos θ incident
t s = ------------------------------------------------------------------------------------
n 1 cos θ incident + n 2 cos θ transmitted
2n 1 cos θ incident
t p = ------------------------------------------------------------------------------------
n 1 cos θ transmitted + n 2 cos θ incident
2
R = r
n 2 cos θ transmitted 2
T = ------------------------------------------ t
n 1 cos θ incident
n2
θ B = atan ------
n1
2 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
Figure 2: Electric field, Ey (slice) and power flow (arrows) for TE incidence at 70° inside the
unit cell.
3 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
For the TM case, Figure 3 visualizes the y component of the magnetic-field distribution
instead, again in combination with the power flow.
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 case—is close to zero.
4 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
References
1. C.A. Balanis, Advanced Engineering Electromagnetics, Wiley, 1989.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
Define some parameters that are useful when setting up the mesh and the study.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
6 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
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.
DEFINITIONS
Variables 1
1 On the Home toolbar, click Variables and choose Local Variables.
2 In the Settings window for Variables, locate the Variables section.
7 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
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 and Shape 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:
8 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
E LE C T RO M AG NE TI C WAVE S, FRE Q UE N C Y D O MA IN ( E MW )
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.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
9 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
0 x
exp(-i*kax*x)[V/m] y
0 z
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 x
exp(-i*kbx*x)[V/m] y
0 z
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
Solved with COMSOL Multiphysics 5.2
kax x
0 y
0 z
11 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
MATERIALS
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
choose Blank Material.
2 In the Settings window for Material, type Air in the Label text field.
3 Select Domain 2 only.
4 Locate the Material Contents section. In the table, enter the following settings:
12 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
Material 2 (mat2)
1 Right-click Materials and choose Blank Material.
2 In the Settings window for Material, type Glass in the Label text field.
3 Select Domain 1 only.
4 Locate the Material Contents section. In the table, enter the following settings:
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
Parametric Sweep
1 On the Study toolbar, click Parametric Sweep.
13 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Parametric Sweep, locate the Study Settings section.
3 Click Add.
4 In the table, enter the following settings:
RESULTS
1 In the Model Builder window, expand the Electric Field (emw) node, then click Multislice
1.
2 In the Settings window for Multislice, click Replace Expression in the upper-right corner
of the Expression section. From the menu, choose Model>Component 1>Electromagnetic
Waves, Frequency Domain>Electric>Electric field>emw.Ey - Electric field, y component.
3 Locate the Multiplane Data section. Find the x-planes subsection. In the Planes text 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.
7 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Energy and
power>emw.Poavx,...,emw.Poavz - Power flow, time average.
8 Locate the Arrow Positioning section. Find the y grid points subsection. In the Points text
field, type 1.
9 Locate the Coloring and Style section. From the Color list, choose Green.
10 In the Model Builder window, click Electric Field (emw).
11 In the Settings window for 3D Plot Group, locate the Data section.
12 From the Parameter value (alpha (1)) list, choose 1.2217.
13 On the Electric Field (emw) toolbar, click Plot.
14 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
14 Click the Zoom Extents button on the Graphics toolbar. The plot should look like 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 Home 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.
Global 1
On the 1D Plot Group 2 toolbar, click Global.
1D Plot Group 2
1 In the Settings window for Global, locate the y-Axis Data section.
2 In the table, enter the following settings:
3 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.
4 Find the Line markers subsection. From the Marker list, choose Cycle.
5 In the Model Builder window, click 1D Plot Group 2.
6 In the Settings window for 1D Plot Group, type Reflection and Transmission in
the Label text field.
Global 2
On the Reflection and Transmission toolbar, click Global.
15 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
3 On the Reflection and Transmission toolbar, click Plot. Compare the resulting plots with
Figure 4.
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 have no effect. Thus, the H-field is perpendicular
to the plane of incidence and it is convenient to specify the port mode fields as the H-field.
Port 1
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw) click Port 1.
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 x
exp(-i*kax*x)[A/m] y
0 z
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 x
exp(-i*kbx*x)[A/m] y
0 z
16 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
STUDY 1
Solution 1 (sol1)
1 In the Model Builder window, expand the Study 1>Solver Configurations node.
2 Right-click Solution 1 (sol1) and choose Solution>Copy.
RESULTS
STUDY 1
On the Home toolbar, click Compute.
17 | FRESNEL EQUATIONS
Solved with COMSOL Multiphysics 5.2
RESULTS
4 On the Reflection and Transmission 1 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
Solved with COMSOL Multiphysics 5.2
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 example 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
1 | H- B E N D WAVE G UI D E 2 D
Solved with COMSOL Multiphysics 5.2
example 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
Solved with COMSOL Multiphysics 5.2
Figure 3 shows the S-parameters in a dB scale. The result agrees very well with that of
the 3D model.
3 | H- B E N D WAVE G UI D E 2 D
Solved with COMSOL Multiphysics 5.2
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 | H- B E N D WAVE GU I D E 2 D
Solved with COMSOL Multiphysics 5.2
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 and Shape 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
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Rectangle, locate the Size and Shape 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 x text field, type -0.1.
6 In the y text field, type -0.08.
7 Click the Build All Objects button.
8 Click the Zoom Extents button on the Graphics toolbar.
MATERIALS
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Air in the Label text field.
3 Select Domains 1 and 3 only.
4 Click to expand the Material properties section. Locate the Material Properties
section. In the Material properties tree, select Electromagnetic Models>Refractive
Index>Refractive index, real part (n).
5 Click Add to Material.
6 Locate the Material Contents section. In the table, enter the following settings:
Material 2 (mat2)
1 Right-click Materials and choose Blank Material.
2 In the Settings window for Material, type Silica Glass in the Label text field.
3 Select Domain 2 only.
4 Click to expand the Material properties section. Locate the Material Properties
section. In the Material properties tree, select Electromagnetic Models>Refractive
Index>Refractive index, real part (n).
5 Click Add to Material.
6 Locate the Material Contents section. In the table, enter the following settings:
6 | H- B E N D WAVE GU I D E 2 D
Solved with COMSOL Multiphysics 5.2
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 7 only.
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Rectangular.
7 | H- B E N D WAVE G UI D E 2 D
Solved with COMSOL Multiphysics 5.2
MESH 1
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
RESULTS
1 In the Settings window for 2D Plot Group, locate the Data section.
2 From the Parameter value (freq (Hz)) list, choose 4.225E9.
3 On the Electric Field (emw) toolbar, click Plot.
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 S-Parameter (emw) toolbar, click Plot.
8 | H- B E N D WAVE GU I D E 2 D
Solved with COMSOL Multiphysics 5.2
9 | H- B E N D WAVE G UI D E 2 D
Solved with COMSOL Multiphysics 5.2
10 | H-BEND WAVEGUIDE 2D
Solved with COMSOL Multiphysics 5.2
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 example examines a TE wave,
one that has no electric field component in the direction of propagation. More
specifically, for this example 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 component—a 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.
With air as the inside medium of the waveguide, the transmission is nearly perfect
throughout the range of operation. In this example, to make the simulation and the
results more interesting, the bend is 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
Solved with COMSOL Multiphysics 5.2
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 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 Maxwell’s 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 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 of
2 | H- B E N D WAVE GU I D E 3 D
Solved with COMSOL Multiphysics 5.2
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.
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 example 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
Solved with COMSOL Multiphysics 5.2
On the input boundary, the Port boundary condition lets you choose which mode to
send in. Any reflected waves having the same shape are 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.
4 | H- B E N D WAVE GU I D E 3 D
Solved with COMSOL Multiphysics 5.2
The two dips in S21 closely correspond to cavity resonances of the dielectric region in
the bend. At these frequencies, the transmission is almost perfect. (Without the
dielectric, the transmission would be nearly as good throughout the frequency range.)
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
5 | H- B E N D WAVE G UI D E 3 D
Solved with COMSOL Multiphysics 5.2
GEOMETRY 1
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 and Shape section.
6 | H- B E N D WAVE GU I D E 3 D
Solved with COMSOL Multiphysics 5.2
Rectangle 2 (r2)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape 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
MATERIALS
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Air in the Label text field.
3 Select Domains 1 and 3 only.
4 Click to expand the Material properties section. Locate the Material Properties
section. In the Material properties tree, select Electromagnetic Models>Refractive
Index>Refractive index, real part (n).
5 Click Add to Material.
7 | H- B E N D WAVE G UI D E 3 D
Solved with COMSOL Multiphysics 5.2
6 Locate the Material Contents section. In the table, enter the following settings:
Material 2 (mat2)
1 Right-click Materials and choose Blank Material.
2 In the Settings window for Material, type Silica Glass in the Label text field.
3 Select Domain 2 only.
4 Click to expand the Material properties section. Locate the Material Properties
section. In the Material properties tree, select Electromagnetic Models>Refractive
Index>Refractive index, real part (n).
5 Click Add to Material.
6 Locate the Material Contents section. In the table, enter the following settings:
8 | H- B E N D WAVE GU I D E 3 D
Solved with COMSOL Multiphysics 5.2
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
choose Build All.
If you look closely at the mesh, you can see that it is indeed a bit finer in the bend
than elsewhere.
STUDY 1
9 | H- B E N D WAVE G UI D E 3 D
Solved with COMSOL Multiphysics 5.2
RESULTS
An S-parameter plot gives you a quantitative measure of how much of the wave is
transmitted and reflected at different frequencies.
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 S-Parameter (emw) 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
Solved with COMSOL Multiphysics 5.2
11 | H- B E N D WAVE G UI D E 3 D
Solved with COMSOL Multiphysics 5.2
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
12 | H-BEND WAVEGUIDE 3D
Solved with COMSOL Multiphysics 5.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:
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
The Touchstone file should now contain the complete output from the model. The
new solution data set contains two frequency sweeps, one for each port.
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.
13 | H- B E N D WAVE G UI D E 3 D
Solved with COMSOL Multiphysics 5.2
5 Click Add Expression in the upper-right corner of the y-axis data section. From the
menu, choose Model>Component 1>Electromagnetic Waves, Frequency
Domain>Ports>S-parameter, dB>emw.S22dB - S-parameter, dB, 22 component.
6 Locate the Data section. From the Parameter selection (PortName) list, choose Last.
7 On the S-Parameter (emw) toolbar, click Plot.
14 | H-BEND WAVEGUIDE 3D
Solved with COMSOL Multiphysics 5.2
15 | H- B E N D WAVE G UI D E 3 D
Solved with COMSOL Multiphysics 5.2
16 | H-BEND WAVEGUIDE 3D
Solved with COMSOL Multiphysics 5.2
Hexagonal Grating
Introduction
A hexagonal grating is an infinite structure that is periodic with hexagonal (or
rhomboid unit) cells. Figure 1 shows the hexagonal domain used for this model. The
reflecting perfectly conducting surface consists of regularly spaced protruding
semispheres.
Figure 1: The hexagonal domain, used for computing the diffraction from the hexagonal
grating.
As shown in Figure 2, for a hexagonal cell of side length a, the corresponding unit cell
is a rhomboid with side length 3a . In Figure 2, the side vectors for the hexagonal
cell starts from the point P and are denoted a1 and a2. The angle between a1 and a2
is 120 degrees. Similarly, for the rhomboid unit cell, the primitive vectors are denoted
1 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
u1 and u2 and starts from the hexagon center point Q. The angle between the two
primitive vectors is also 120 degrees.
Figure 2: Schematic showing the hexagonal cells with side length a and side vectors a1 and
a2. The primitive cells are defined by the primitive vectors u1 and u2.
k = k || + k ⊥ , (1)
where k|| is the wave vector component parallel to the periodic boundary and k⊥ is
the component orthogonal to the periodic boundary, the in-plane wave vector
component for diffraction order mn is given by
k ||mn = k || + mG 1 + nG 2 , (2)
where the reciprocal lattice vectors G1 and G2 are defined from the primitive vectors
u1 and u2 as
u2 × n
G 1 = 2π --------------------------------- (3)
u1 ⋅ ( u2 × n )
and
n × u1
G 2 = 2π --------------------------------- , (4)
u2 ⋅ ( n × u1 )
2 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
2 2
k ⊥mn = k – k ||mn , (5)
it is clear that for propagating modes, where the out-of-plane wave vector component
above must be real, the length of the in-plane wave vector component must be smaller
than the material wave number k. Figure 3 shows that it is only the modes inside the
circle with radius k that will be propagating. In the example shown in Figure 3, there
are five modes that will be propagating, in this case the modes m = n = 0 (the reflected
wave), m = -1, n = 0, m = 0, n = -1, m = -1, n = -1, and m = -2, n = -1. All other
modes will be evanescent and damped out.
Figure 3: The reciprocal lattice, showing the reciprocal lattice vectors G1 and G2, the
in-plane wave vector component k||, and the circle with radius k (the material wave
number) enclosing the propagating mode points (larger dots. The dotted hexagon indicates
that also the reciprocal lattice is a hexagonal point lattice. The dashed rhomboid indicates
the unit cell spanned by the reciprocal lattice vectors.
Model Definition
In this model, the unit cell is small compared to the wavelength, so there will only be
two modes that are propagating, the modes m = 0, n = -1 and m = -1, n = -1. For
wavelengths longer than approximately 1.01 μm (the critical wavelength), the mode
m = 0, n = -1 will be evanescent.
First a wavelength sweep will be made for an incident field having the polarization
perpendicular to the plane of incidence (spanned by the wave vector for the incident
wave and the normal to the periodic boundary) (so called s-polarization). Thereafter
3 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
another wavelength sweep is made, but now with the polarization in the plane of
incidence (p-polarization).
Figure 4: The electric field norm and the propagation directions for the incident wave (red
arrows), the reflected wave (blue arrows) and the two diffraction orders (green and yellow
arrows). The wavelength is 1.01 μm, which is close to the critical wavelength for the mode
m = 0, n = -1, and the polarization of the incident wave is perpendicular to the plane of
incidence.
Figure 5 shows the reflectance (for mode m = n = 0) and the diffraction efficiencies for
the diffracted waves. Notice that both the reflectance and the diffraction efficiency for
4 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Figure 5: Diffraction efficiencies for the reflected wave and the diffracted waves. The
polarization of the incident wave is perpendicular to the plane of incidence.
5 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Figure 6 shows a similar plot as Figure 4, but here the polarization of the incident wave
is parallel with the plane of incidence.
Figure 6: Similar plot as in Figure 4, but here the polarization of the incident wave is
parallel to the plane of incidence.
6 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Figure 7 shows that for p-polarization both the reflected wave and the two
m = -1, n = -1 modes show resonances close to the critical wavelength for the
m = 0, n = -1 mode.
Figure 7: Similar plot as in Figure 5, but here the polarization of the incident wave is
parallel to the plane of incidence.
Model Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
7 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
GLOBAL DEFINITIONS
First add some parameters that defines the geometry and the incident electric field.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Notice that the azimuth angle phi above measures the angle for the wave vector of
the incident wave from the x-axis. Here, c_const is a predefined COMSOL
constant for the speed of light in vacuum.
GEOMETRY 1
The geometry consists of an extruded hexagon, with a semisphere removed from it at
the bottom.
8 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Polygon 1 (pol1)
1 On the Geometry toolbar, click Primitives and choose Polygon.
2 In the Settings window for Polygon, locate the Coordinates section.
3 In the xw text field, type a a/2 -a/2 -a -a/2 a/2.
4 In the yw text field, type 0 sqrt(3)/2*a sqrt(3)/2*a 0 -sqrt(3)/2*a
-sqrt(3)/2*a.
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)
h0
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 a1.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object ext1 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 sph1 only.
6 Click the Build All Objects button.
7 Click the Go to Default 3D View button on the Graphics toolbar.
9 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
MATERIALS
On the Home toolbar, click Windows and choose Add Material.
ADD MATERIAL
1 Go to the Add Material window.
2 In the tree, select Built-In>Air.
3 Click Add to Component 1.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
10 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
3 In the Settings window for Port, locate the Port Properties section.
4 From the Type of port list, choose Periodic.
5 From the Wave excitation at this port list, choose On.
6 Locate the Port Mode Settings section. Specify the E0 vector as
-E0*sin(phi) x
E0*cos(phi) y
0 z
Before creating the diffraction orders ports, a reference point must be defined on the
periodic port.
11 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Periodic Port Reference Point, locate the Point Selection
section.
3 Click Clear Selection.
4 Select Point 2 only. This point selection makes the angle previously provided for α2
consistent with the intended angle of incidence for the incident wave.
Port 1
1 In the Model Builder window, under Component 1 (comp1)>Electromagnetic Waves,
Frequency Domain (emw) click Port 1.
2 In the Settings window for Port, locate the Port Mode Settings section.
3 Click Compute Diffraction Orders.
Periodic Condition 1
1 On the Physics toolbar, click Boundaries and choose Periodic Condition.
2 In the Settings window for Periodic Condition, locate the Periodicity Settings
section.
3 From the Type of periodicity list, choose Floquet periodicity.
4 From the k-vector for Floquet periodicity list, choose From periodic port.
12 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Periodic Condition 2
1 On the Physics toolbar, click Boundaries and choose Periodic Condition.
2 In the Settings window for Periodic Condition, locate the Periodicity Settings
section.
3 From the Type of periodicity list, choose Floquet periodicity.
4 From the k-vector for Floquet periodicity list, choose From periodic port.
13 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Periodic Condition 3
1 On the Physics toolbar, click Boundaries and choose Periodic Condition.
2 In the Settings window for Periodic Condition, locate the Periodicity Settings
section.
3 From the Type of periodicity list, choose Floquet periodicity.
4 From the k-vector for Floquet periodicity list, choose From periodic port.
14 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
15 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Build All.
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:
5 Click Range.
6 In the Range dialog box, type 0.9[um] in the Start text field.
7 In the Step text field, type 0.01[um].
8 In the Stop text field, type 1.1[um].
9 Click Replace.
16 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
RESULTS
1 In the Model Builder window, under Results right-click Electric Field (emw) and
choose Arrow Surface.
2 In the Settings window for Arrow Surface, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Component
1>Electromagnetic Waves, Frequency Domain>Ports>emw.kIncx_1,...,emw.kIncz_1 -
Incident wave vector.
3 Locate the Expression section. Select the Description check box.
4 In the associated text field, type Incident wave (red).
5 Right-click Results>Electric Field (emw)>Arrow Surface 1 and choose Duplicate.
6 In the Settings window for Arrow Surface, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Component
1>Electromagnetic Waves, Frequency
Domain>Ports>emw.kModex_1,...,emw.kModez_1 - Port mode wave vector.
7 Locate the Coloring and Style section. From the Color list, choose Blue.
8 Locate the Expression section. Select the Description check box.
9 In the associated text field, type Reflected wave (blue).
10 Right-click Results>Electric Field (emw)>Arrow Surface 2 and choose Duplicate.
11 In the Settings window for Arrow Surface, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Component
1>Electromagnetic Waves, Frequency
Domain>Ports>emw.kModex_2,...,emw.kModez_2 - Port mode wave vector.
12 Locate the Coloring and Style section. From the Color list, choose Green.
17 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
18 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
Now, add a plot showing the reflectance and the diffraction efficiencies for the
diffracted waves.
1D Plot Group 2
1 On the Home 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 (sol2).
Global 1
On the 1D Plot Group 2 toolbar, click Global.
1D Plot Group 2
1 In the Settings window for Global, locate the y-Axis Data section.
2 In the table, enter the following settings:
3 Locate the x-Axis Data section. From the Axis source data list, choose Outer solutions.
4 From the Parameter list, choose Expression.
5 In the Expression text field, type wl.
6 Select the Description check box.
7 In the associated text field, type Wavelength.
8 Click to expand the Legends section. From the Legends list, choose Manual.
9 In the table, enter the following settings:
Legends
Reflectance
m = -1, n = -1, in-plane
m = -1, n = -1, out-of-plane
m = 0, n = -1, in-plane
m = 0, n = -1, out-of-plane
19 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
11 In the Settings window for 1D Plot Group, locate the Title section.
12 From the Title type list, choose None.
13 Locate the Plot Settings section. Select the y-axis label check box.
14 In the associated text field, type Diffraction efficiency.
15 Click to expand the Legend section. From the Position list, choose Upper left.
16 On the 1D Plot Group 2 toolbar, click Plot. Your plot should look like Figure 5.
Port 1
Now, repeat the simulation for an incoming wave with p-polarization (the electric field
polarized in the plane of incidence).
20 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
-H0*sin(phi) x
H0*cos(phi) y
0 z
Study 1
On the Home toolbar, click Compute.
RESULTS
21 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
1D Plot Group 2
22 | HEXAGONAL GRATING
Solved with COMSOL Multiphysics 5.2
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 example.
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.
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:
∇ × ( μ r ∇ × E ) – k 0 ε r – --------- E = 0
–1 2 jσ
ωε 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.
μ 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:
2ω 0 ω m ω2 α ωω m ( ω20 – ω2 ( 1 + α2 ) )
2
- + j -------------------------------------------------------------------------------------
χ xy = ------------------------------------------------------------------------------------- 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,
with an effective loss tangent of 2·10-4 and ΔH = 3.18·103 A/m, are taken for
aluminum garnet from Ref. 2. The applied bias field is set to H0 = 7.96·103 A/m. The
electron gyromagnetic ratio taken from Ref. 2 is 1.759·1011 C/kg.
The plot of the S-parameter from the parametric sweep of sc_ferrite indicates a
minimum for a scale factor of 0.518.
The plot of the S-parameter from the parametric sweep of sc_chamfer indicates a
minimum for a scale factor of about 3.0.
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.
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 application’s Application Library folder and double-click the file
lossy_circulator_3d_geom.mph.
GEOMETRY 1
ADD MATERIAL
1 On the Home 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 Home 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.
9 Locate the Conduction Current section. From the σ list, choose User defined.
One inport for excitation and two outports need to be added next.
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.
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.
The final step in the model set up is to solve it for the nominal frequency and inspect
the results for possible modeling errors.
STUDY 1
RESULTS
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
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.
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:
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,
choose Only last.
7 On the Study toolbar, click Compute.
RESULTS
TA BL E
1 Go to the Table window.
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.
GLOBAL 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:
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:
TABLE
Go to the Table window.
RESULTS
TA BL E
1 Go to the Table window.
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
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.
GLOBAL 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:
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
TABLE
Go to the Table window.
RESULTS
1 In the Settings window for 3D Plot Group, locate the Data section.
2 From the Parameter value (freq (Hz)) list, choose 3E9.
3 On the Electric Field (emw) 1 toolbar, click Plot.
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.
3 Click to select row number 2 in the table.
4 Click Delete twice.
5 Click to select row number 3 in the table.
6 Click Delete twice.
7 Click to select row number 4 in the table.
8 Click Delete twice.
The table should now only contain emw.S11dB, emw.S21dB and emw.S31dB.
RESULTS
Reuse the first study for the port sweep. The study is solved for a single frequency to
keep down simulation time though it is possible to solve for a range of frequencies.
STUDY 1
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.
4 Locate the Output While Solving section. Find the Memory settings for jobs
subsection. From the Keep solutions in memory list, choose All.
5 On the Home toolbar, click Compute.
TABLE
Go to the Table window.
RESULTS
Derived Values
1 In the Settings window for Global Matrix Evaluation, locate the Data section.
2 From the Data set list, choose Study 1/Parametric Solutions 1 (sol3).
3 Click Replace Expression in the upper-right corner of the Expression section. From
the menu, choose Component 1>Electromagnetic Waves, Frequency
Domain>Ports>S-parameter, dB>emw.SdB - S-parameter, dB.
4 Click the Evaluate button.
TABLE
Go to the Table window.
RESULTS
DEFINITIONS
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
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
The geometry is set up using a parameterized approach. This allows you to match the
input impedance to that of the connecting waveguide sections by variation of two
geometric design parameters. Before starting to build the geometry the geometric
design parameters need to be entered. These are two dimensionless numbers used to
scale selected geometric building blocks.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
The lossy ferrite material model is set up by referring to global variables. For
convenience the definitions are stored in an external text file that is imported into the
model. The external text file also contains comments.
Variables 1
1 On the Home toolbar, click Variables and choose Global Variables.
2 In the Settings window for Variables, locate the Variables section.
3 Click Load from File.
4 Browse to the application’s Application Library folder and double-click the file
lossy_circulator_3d_parameters.txt.
The geometry is built by first defining a 2D cross section of the 3D geometry in a work
plane. The 2D geometry is then extruded into 3D.
GEOMETRY 1
Plane Geometry
Start by defining one arm of the circulator, then twice copy and rotate it to build all
three arms.
Rectangle 1 (r1)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type 0.2-0.1/(3*sqrt(3)).
4 In the Height text field, type 0.2/3.
5 Locate the Position section. In the xw text field, type -0.2.
6 In the yw text field, type -0.1/3.
Copy 1 (copy1)
1 On the Work Plane toolbar, click Transforms and choose Copy.
2 Select the object r1 only.
Rotate 1 (rot1)
1 On the Work Plane toolbar, click Transforms and choose Rotate.
2 Select the object copy1 only.
3 In the Settings window for Rotate, locate the Rotation Angle section.
4 In the Rotation text field, type 120.
5 Right-click Rotate 1 (rot1) and choose Build Selected.
Copy 2 (copy2)
1 On the Work Plane toolbar, click Transforms and choose Copy.
2 Select the object r1 only.
Rotate 2 (rot2)
1 On the Work Plane toolbar, click Transforms and choose Rotate.
2 Select the object copy2 only.
3 In the Settings window for Rotate, locate the Rotation Angle section.
4 In the Rotation text field, type -120.
5 Right-click Rotate 2 (rot2) and choose Build Selected.
6 Click the Zoom Extents button on the Graphics toolbar.
Plane Geometry
The geometry should now look as in the below figure.
Union 1 (uni1)
1 On the Work Plane toolbar, click Booleans and Partitions and choose Union.
2 Click in the Graphics window and then press Ctrl+A to select all objects.
Plane Geometry
1 Right-click Union 1 (uni1) and choose Build Selected.
Next build the central connecting region and add the ferrite domain. During these
stages, the geometric design parameters will be used. First build/add a triangle
connecting the arms by subtracting a copy of what has already been drawn from a circle
of proper radius.
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.2/(3*sqrt(3)).
4 Right-click Circle 1 (c1) and choose Build Selected.
Copy 3 (copy3)
1 On the Work Plane toolbar, click Transforms and choose Copy.
2 Select the object uni1 only.
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 object copy3 only.
Plane Geometry
1 Right-click Difference 1 (dif1) and choose Build Selected.
The geometry should now look as in the below figure.
Rotate the newly created triangle 180 degrees and use one scaled copy of it to create
linear fillets for impedance matching. Use another scaled copy to define the ferrite.
Rotate 3 (rot3)
1 On the Work Plane toolbar, click Transforms and choose Rotate.
2 Select the object dif1 only.
3 In the Settings window for Rotate, locate the Rotation Angle section.
4 In the Rotation text field, type 180.
5 Right-click Rotate 3 (rot3) and choose Build Selected.
Copy 4 (copy4)
1 On the Work Plane toolbar, click Transforms and choose Copy.
2 Select the object rot3 only.
Plane Geometry
Apply the scaling for the impedance matching.
Scale 1 (sca1)
1 On the Work Plane toolbar, click Transforms and choose Scale.
2 In the Settings window for Scale, locate the Scale Factor section.
3 In the Factor text field, type sc_chamfer.
4 Select the object copy4 only.
5 Right-click Scale 1 (sca1) and choose Build Selected.
Union 2 (uni2)
1 On the Work Plane toolbar, click Booleans and Partitions and choose Union.
2 Select the objects sca1 and uni1 only.
3 In the Settings window for Union, locate the Union section.
4 Clear the Keep interior boundaries check box.
Plane Geometry
1 Right-click Union 2 (uni2) and choose Build Selected.
The geometry should now look as in the below figure.
Scale 2 (sca2)
1 On the Work Plane toolbar, click Transforms and choose Scale.
2 Select the object rot3 only.
3 In the Settings window for Scale, locate the Scale Factor section.
4 In the Factor text field, type sc_ferrite.
Extrude 1 (ext1)
1 On the Geometry toolbar, click Extrude.
2 In the Settings window for Extrude, locate the Distances from Plane section.
Distances (m)
0.1/3
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.
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 directions—into 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.
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
F y ( X g, Y g )
Y g 1 + --- x
y 1 2
2
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.
Figure 4 displays the dielectric distribution, and shows variation in both directions
defined via the mapping described above.
Figure 3: The norm of the electric field shows the Gaussian beam diffracted by the
metamaterial lens.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D.
2 In the Select physics tree, select Mathematics>Deformed Mesh>Deformed Geometry
(dg).
3 Click Add.
4 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
5 Click Add.
6 Click Study.
7 In the Select study tree, select Custom Studies>Preset Studies for Some Physics
Interfaces>Stationary.
8 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
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:
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Height text field, type 2.
4 Locate the Position section. From the Base list, choose Center.
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.
2 In the Settings window for Explicit, type Lens in the Label text field.
Next, add a set of variables for the shape and the dielectric distribution of the
metamaterial lens.
Variables 1
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:
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.
1 On the Physics toolbar, click Boundaries and choose Prescribed Mesh Displacement.
2 Select Boundaries 15–18 only.
3 In the Settings window for Prescribed Mesh Displacement, locate the Prescribed
Mesh Displacement section.
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.
3 In the Settings window for Surface Current Density, locate the Surface Current
Density section.
4 Specify the Js0 vector as
0 x
0 y
exp(-(y/w0)^2) z
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 Home 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 Home 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
choose Free Triangular.
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 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
in a few times to check the quality of the mesh.
STUDY 1
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:
Add a Frequency Domain study step and set as solved only for Electromagnetic Waves,
Frequency Domain.
Frequency Domain
On the Study toolbar, click Study Steps and choose Frequency Domain>Frequency
Domain.
3 Locate the Physics and Variables Selection section. In the table, enter the following
settings:
RESULTS
1 In the Model Builder window, expand the Electric Field (emw) node, then click Surface
1.
2 In the Settings window for Surface, locate the Coloring and Style section.
3 From the Color table list, choose RainbowLight.
4 In the Model Builder window, right-click Electric Field (emw) and choose Contour.
5 In the Settings window for Contour, locate the Levels section.
6 In the Total levels text field, type 14.
7 Locate the Coloring and Style section. From the Coloring list, choose Uniform.
8 From the Color list, choose Black.
9 Clear the Color legend check box.
10 On the Electric Field (emw) toolbar, click Plot. See Figure 3 to compare the
reproduced plot.
Add a filled contour plot describing the dielectric distribution over the lens.
2D Plot Group 2
1 On the Home toolbar, click Add Plot Group and choose 2D Plot Group.
2 In the Model Builder window, right-click 2D Plot Group 2 and choose Contour.
3 In the Settings window for Contour, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Component
1>Electromagnetic Waves, Frequency Domain>Material properties>emw.epsrAv -
Relative permittivity, average.
4 Locate the Levels section. In the Total levels text field, type 12.
5 Locate the Coloring and Style section. From the Contour type list, choose Filled.
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
Solved with COMSOL Multiphysics 5.2
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 User’s
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:
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
Solved with COMSOL Multiphysics 5.2
Figure 3 shows the temperature in the center of the potato as a function of time for
the first 5 seconds. Due to the low thermal conductivity of the potato, the heat
distributes rather slowly, and the temperature profile after 5 seconds has a strong peak
in the center (see Figure 4). When heating the potato further, the temperature in the
center eventually reaches 100 °C and the water contents start boiling, drying out the
center and transporting heat as steam to outer layers. This also affects the
electromagnetic properties of the potato. The simple microwave absorption and heat
conduction model used here does not capture these nonlinear effects. However, the
model can serve as a starting point for a more advanced analysis.
3 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
Figure 3: Temperature in the center of the potato during the first 5 seconds of heating.
4 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Heat Transfer>Electromagnetic Heating>Microwave
Heating.
3 Click Add.
4 Click Study.
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.
Time Dependent
On the Study toolbar, click Study Steps and choose Time Dependent>Time Dependent.
1 In the Settings window for Time Dependent, locate the Physics and Variables
Selection section.
5 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
GLOBAL DEFINITIONS
First, define a set of parameters for creating the geometry.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
6 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
GEOMETRY 1
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size and Shape 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 and Shape 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
Solved with COMSOL Multiphysics 5.2
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.
The object cyl1 is the cylinder.
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.
The object sph1 is the sphere.
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 Click the Build All Objects button.
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.
8 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
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 View > Selection List.
Explicit 1
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Potato in the Label text field.
3 Select Domain 3 only.
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Plate in the Label text field.
3 Select Domain 2 only.
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Air in the Label text field.
3 Select Domains 1 and 4 only.
Explicit 4
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type No Heat Transfer in the Label text field.
3 Select Domains 1, 2, and 4 only.
Explicit 5
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Port Boundary in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundary 23 only.
Explicit 6
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Symmetry Boundaries in the Label text
field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
9 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
Explicit 7
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Metal Boundaries in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundaries 1, 3–5, 17, and 20–22 only.
MATERIALS
Next, define the materials. Air and Copper are already in the Material Library.
ADD MATERIAL
1 On the Home 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, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Potato in the Label text field.
3 Locate the Geometric Entity Selection section. From the Selection list, choose Potato.
4 Locate the Material Contents section. In the table, enter the following settings:
10 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
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 Go to the Add Material window.
2 In the tree, select Built-In>Copper.
3 Click Add to Component in the window toolbar.
MATERIALS
On the Home 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.
11 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
3 Select the Enable check box.
Set the maximum mesh size to 0.2 wavelengths or smaller.
4 In the Maximum element size text field, type c_const/2.45[GHz]/5.
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 x
0 y
cos(pi*y/dg)[V/m] z
12 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
The Heat Transfer physics will automatically use the electromagnetic heat source from
the Electromagnetic Waves physics thanks to the Electromagnetic Heat Source
coupling feature.
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)
click Initial Values 1.
2 In the Settings window for Initial Values, locate the Initial Values section.
3 In the T text field, type T0.
MESH 1
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
13 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
Solution 1 (sol1)
1 On the Study toolbar, click Show Default Solver.
2 In the Model Builder window, expand the Solution 1 (sol1) node.
3 Right-click Solution 1 (sol1) and choose Other>Solution Store.
4 Move Solution Store 1 to the position below Stationary Solver 1.
5 On the Study toolbar, click Compute.
RESULTS
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 Temperature (ht) toolbar, click Plot.
Modify an existing plot group to plot the resistive heating on the symmetry plane.
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
of the Expression section. From the menu, choose Component 1>Electromagnetic
Waves, Frequency Domain>Heating and losses>emw.Qrh - Resistive losses.
5 Locate the Plane Data section. From the Plane list, choose zx-planes.
14 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
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
of the Expression section. From the menu, choose Component 1>Electromagnetic
Waves, Frequency Domain>Electric>Electric field>emw.Ez - Electric field, z component.
3 Locate the Plane Data section. From the Plane list, choose xy-planes.
4 From the Entry method list, choose Coordinates.
5 In the z-coordinates text field, type 0.1.
6 Right-click Results>Temperature (ht)>Slice 1 and choose Deformation.
7 In the Settings window for Deformation, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Component
1>Electromagnetic Waves, Frequency Domain>Electric>emw.Ex,emw.Ey,emw.Ez -
Electric field.
8 On the Temperature (ht) toolbar, click Plot.
Derived Values
Make a volume integral of the microwave heating to find out how much of the energy
is absorbed in the potato.
Volume Integration 1
On the Results toolbar, click More Derived Values and choose Integration>Volume
Integration.
Derived Values
1 In the Settings window for Volume Integration, locate the Data section.
2 From the Data set list, choose Study 1/Solution Store 1 (sol2).
3 Click Replace Expression in the upper-right corner of the Expression section. From
the menu, choose Component 1>Heat Transfer in Solids>Heat sources>ht.Qtot - Total
heat source.
Select one point in time for the output. Since the material parameters of the potato
are independent of the temperature, it does not matter which time you choose.
15 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
4 Locate the Data section. From the Parameter selection (freq) list, choose From list.
5 Locate the Selection section. From the Selection list, choose Potato.
6 Click the Evaluate button.
TABLE
1 Go to the Table window.
The result is 312 W. Finally, to reproduce Figure 3, create a plot of temperature in
the center of the potato as a function of time.
RESULTS
Cut Point 3D 1
On the Results toolbar, click Cut Point 3D.
Data Sets
1 In the Settings window for Cut Point 3D, locate the Point Data section.
2 In the x text field, type 0.134.
3 In the y text field, type 0.
4 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).
Point Graph 1
On the 1D Plot Group 4 toolbar, click Point Graph.
1D Plot Group 4
1 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.
2 Locate the y-Axis Data section. From the Unit list, choose degC.
3 On the 1D Plot Group 4 toolbar, click Plot.
16 | MICROWAVE OVEN
Solved with COMSOL Multiphysics 5.2
M i c r o w a v e Fi lt er on PCB
This example 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 example 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
Electromagnetic Waves interface without any mechanical deformation. Then a
uniform load of 40 N is applied to the PCB, generating a large deformation of the
board. The Solid Mechanics interface calculates the deformation, and the Moving
Mesh interface applies this deformation to the coordinate system that the
Electromagnetic Waves interface uses. After this step, the frequency sweep is
performed again for the Electromagnetic Waves interface using the parametric solver.
This example accounts for the structural deformation in the sense that it solves for the
electromagnetic fields on the deformed geometry, as if the PCB was manufactured in
the deformed shape—free of stress.
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
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, pp. 45–54, 2006.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
The following steps define the parameters for the frequency sweep.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
GEOMETRY 1
Set mm as the default unit for length.
Block 1 (blk1)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, locate the Size and Shape section.
3 In the Width text field, type 100.
4 In the Depth text field, type 40.
5 In the Height text field, type 15.
6 Locate the Position section. In the x text field, type -5.
7 In the y text field, type -15.
8 In the z text field, type -10.
Import 1 (imp1)
1 On the Geometry toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 From the Source list, choose ECAD file (ODB++).
4 Click Browse.
5 Browse to the application’s Application Library folder and double-click the file
pcb_microwave_filter_with_stress.xml.
It is located in the same folder as specified in the Application Library Path on page
3. This path originates from the folder applications 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.
The next step is to add boundary faces for the input and output ports.
Rectangle 1 (r1)
1 On the Geometry 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 and Shape 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 Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape 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 3.645.
6 In the yw text field, type -0.635.
ADD MATERIAL
1 On the Home 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
1 Go to the Add Material window.
2 In the tree, select Built-In>FR4 (Circuit Board).
3 Click Add to Component in the window toolbar.
MATERIALS
5 On the Home toolbar, click Add Material to close the Add Material window.
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.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Select Boundary 16 only.
MESH 1
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose More Operations>Free Triangular.
2 Select Boundaries 10 and 16 only.
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.
The Maximum element size is reduced to 1 mm to increase the accuracy of the
S-parameter calculations.
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
choose Global Variable Probe.
2 In the Settings window for Global Variable Probe, type S11 in the Variable name text
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
the upper-right corner of the section.
STUDY 1
TABLE
Go to the Table window.
RESULTS
The Probe 1D Plot Group 2 displays the S11- and S21-parameters for the frequency
sweep.
ADD PHYSICS
1 On the Home 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
1 Go to the Add Physics window.
2 In the Add physics tree, select Mathematics>Deformed Mesh>Moving Mesh (ale).
3 Click Add to Component in the window toolbar.
4 On the Home toolbar, click Add Physics to close the Add Physics window.
GLOBAL 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:
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.SIGNAL_LAYER_DIEL, 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.
GLOBAL DEFINITIONS
In the Model Builder window, expand the Component 1 (comp1)>Moving Mesh (ale)
node, then click 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.
3 In the Settings window for Body Load, locate the Force section.
4 Specify the FV vector as
0 x
0 y
-fload/V z
Fixed Constraint 1
1 On the Physics toolbar, click Boundaries and choose Fixed Constraint.
2 Select Boundaries 6, 10, 12, and 15–17 only.
STUDY 1
ADD STUDY
1 On the Home 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
Electromagnetic Waves, Frequency Domain (emw)
STUDY 2
Frequency Domain
On the Study toolbar, click Study Steps and choose Frequency Domain>Frequency
Domain.
TABLE
Go to the Table window.
RESULTS
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.
Global 2
On the S-Parameter (emw) 1 toolbar, click Global.
S-Parameter (emw) 1
1 In the Settings window for Global, locate the Data section.
2 From the Data set list, choose Study 1/Solution 1 (sol1).
3 Click Add Expression in the upper-right corner of the y-axis data section. From the
menu, choose emw.S11dB - S-parameter, dB, 11 component.
4 Click Add Expression in the upper-right corner of the y-axis data section. From the
menu, choose emw.S21dB - S-parameter, dB, 21 component.
5 On the S-Parameter (emw) 1 toolbar, click Plot.
You should now see the plot in Figure 2.
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
of the Expression section. From the menu, choose Component 1>Solid
Mechanics>Displacement>solid.disp - Total displacement.
6 In the Model Builder window, click Stress (solid).
7 In the Settings window for 3D Plot Group, locate the Data section.
8 From the Data set list, choose Study 2/Solution Store 1 (sol3).
9 Locate the Plot Settings section. From the Frame list, choose Spatial (x, y, z).
10 On the Stress (solid) toolbar, click Plot.
You should now see the plot in Figure 1.
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 are reflections from the transmission line’s end.
Z in Z-
-------
- = ------ (1)
Z ZL
If the length and the impedance requirements are fulfilled, the load impedance does
not give rise to any reflections.
Z in = Z 0 = 50 Ω
1 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
λ/4
Transmission line
Zin ZL
Model Definition
The 1D geometry of the example 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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 1D.
2 In the Select physics tree, select Radio Frequency>Transmission Line (tl).
4 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
First add some parameters that defines the electrical and geometrical properties of the
transmission lines.
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
5 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
GEOMETRY 1
Set up the geometry as two intervals.
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.
6 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
4 In the Settings window for Transmission Line Equation, locate the Transmission Line
Equation section.
5 In the L text field, type L2.
6 In the C text field, type C2.
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,
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.
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.
7 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
4 Locate the Element Size Parameters section. In the Maximum element size text field,
type hmax.
STUDY 1
1 In the Model Builder window, under Study 1 click Step 1: Frequency Domain.
2 In the Settings window for Frequency Domain, locate the Study Settings section.
3 In the Frequencies text field, type f.
4 On the Home toolbar, click Compute.
RESULTS
1 In the Model Builder window, expand the Electric Potential (tl) node, then click Line
Graph 1.
2 In the Settings window for Line Graph, locate the y-Axis Data section.
3 In the Expression text field, type abs(V).
4 Click Replace Expression in the upper-right corner of the x-axis data section. From
the menu, choose Model>Geometry>Coordinate>x - x-coordinate.
5 On the Electric Potential (tl) toolbar, click Plot.
You should now have a graph as in Figure 2. 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
8 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
GLOBAL 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:
STUDY 1
RESULTS
1D Plot Group 2
On the Home toolbar, click Add Plot Group and choose 1D Plot Group.
Global 1
On the 1D Plot Group 2 toolbar, click Global.
1D Plot Group 2
1 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.
9 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
2 Click Add 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.S21dB - S-parameter, dB, 21 component.
3 On the 1D Plot Group 2 toolbar, click Plot.
Modify the y-axis limits to show that S21 actually has its maximum value for the
frequency where S11 is at its minimum value.
4 In the Model Builder window, click 1D Plot Group 2.
5 In the Settings window for 1D Plot Group, click to expand the Axis section.
6 Select the Manual axis limits check box.
7 In the y minimum text field, type -30.
8 In the y maximum text field, type 2.
Move the legend panel, so it doesn't cover the curves.
9 Click to expand the Legend section. From the Position list, choose Lower right.
10 On the 1D Plot Group 2 toolbar, click Plot.
You should now have a plot of the spectrum for S11 and S21, similar to the one in
Figure 3.
1 In the Model Builder window, under Results click Electric Potential (tl).
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Parameter selection (freq) list, choose Last.
4 On the Electric Potential (tl) toolbar, click Plot.
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
10 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
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.
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.
Compute the spectral plot again.
STUDY 1
On the Home toolbar, click Compute.
11 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
RESULTS
1 In the Settings window for 1D Plot Group, locate the Data section.
2 From the Parameter selection (freq) list, choose From list.
3 In the Parameter values (freq (Hz)) list, select 1E9.
4 In the Parameter values (freq (Hz)) list, select 1E9.
5 On the Electric Potential (tl) toolbar, click Plot.
Notice that there is now a standing wave also at the center frequency.
12 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
14 | Q U A R T E R - W A V E TR A N S F O R M E R
Solved with COMSOL Multiphysics 5.2
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.
PEC sphere
k E
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
where the scattered electric field Esc is the dependent variable and the incident electric
field Ei= (0,0,Ez), with
– jk0 x
E z = 1 [ V ⁄ m ]e
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
n×E = 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.
n×H = 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.”
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.
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.”
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 ( θ, ϕ ) .
• E and H are the fields on the “aperture”—the 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:
η = μ⁄ε
The unit vector r0 can be interpreted as the direction defined by the angular position
( θ, ϕ ) and Ep is the far field in this direction.
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 does 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
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.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
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.
2 In the Settings window for Sphere, locate the Size section.
3 In the Radius text field, type r0+t_air+t_pml.
4 Click to expand the Layers section. In the table, enter the following settings:
DEFINITIONS
Add a view with a different angle of perspective.
1 In the Model Builder window, under Component 1 (comp1) right-click Definitions and
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.
Choose wireframe rendering to get a better view of the interior parts.
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.
4 On the object sph1, select Domains 1–3, 5–7, and 9–15 only.
DEFINITIONS
View 1
After removing unnecessary domains, change the view to the first view definition
which gives a better angle showing all layers.
0 x
0 y
E0*exp(-j*k0*x) z
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.
3 In the Settings window for Perfectly Matched Layer, locate the Geometry section.
4 From the Type list, choose Spherical.
Far-Field Domain 1
1 On the Physics toolbar, click Domains and choose Far-Field Domain.
2 Select Domains 2 and 3 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 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
Next, assign material properties. Use air for all domains.
ADD MATERIAL
1 On the Home 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 Home toolbar, click Add Material to close the Add Material window.
MESH 1
Use a tetrahedral mesh for the air domains.
Free Tetrahedral 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Free Tetrahedral.
2 In the Settings window for Free Tetrahedral, locate the Domain Selection section.
3 From the Geometric entity level list, choose Domain.
4 Select Domains 2 and 3 only.
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.
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:
RESULTS
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 Home 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 (sol2).
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>).
11 Select the y-axis label check box.
12 In the associated text field, type Normalized monostatic RCS
(sigma<sub>3-D</sub>/pi a<sup>2</sup>).
13 Locate the Axis section. Select the y-axis log scale check box.
Point Graph 1
On the 1D Plot Group 1 toolbar, click Point Graph.
1D Plot Group 1
1 Select Point 2 only.
2 In the Settings window for Point Graph, locate the y-Axis Data section.
3 In the Expression text field, type 4*pi*(emw.normEfar^2/E0^2)/(pi*r0^2).
4 Locate the x-Axis Data section. From the Axis source data list, choose Outer solutions.
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 None.
6 Find the Line markers subsection. From the Marker list, choose Square.
7 From the Positioning list, choose In data points.
Next, proceed to perform the mesh convergence study at the first resonance in the Mie
region.
Start by extending the parameter list with the resonant radius and the associated
theoretical RCS value.
GLOBAL 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 Add the two rows at the end and change the first line according to:
MESH 1
Add a new mesh with some tweaks to make sure that the curvature of the sphere is
always resolved.
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Duplicate.
MESH 2
In the Model Builder window, expand the Component 1 (comp1)>Meshes node.
COMPONENT 1 (COMP1)
In the Model Builder window, expand the Component 1 (comp1)>Meshes node.
MESH 2
Size
1 In the Model Builder window, expand the Component 1 (comp1)>Meshes>Mesh 2
node, then click Size.
2 In the Settings window for Size, locate the Element Size Parameters section.
3 In the Minimum element size text field, type r0/2.
Size 1
1 In the Model Builder window, under Component 1 (comp1)>Meshes>Mesh 2
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 Boundary.
4 Click Paste Selection.
5 In the Paste Selection dialog box, type 7,10 in the Selection text field.
6 Click OK.
7 In the Settings window for Size, locate the Element Size section.
8 Click the Custom button.
9 Locate the Element Size Parameters section. Select the Maximum element size check
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 Home 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 Home toolbar, click Add Study to close the Add Study window.
STUDY 2
A parametric sweep is needed to loop over the mesh sizes.
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:
2 In the Settings window for Frequency Domain, locate the Study Settings section.
3 In the Frequencies text field, type f0.
4 In the Model Builder window, click Study 2.
5 In the Settings window for Study, locate the Study Settings section.
6 Clear the Generate default plots check box.
7 On the Study toolbar, click Compute.
RESULTS
Continue to plot the relative error versus elements per wavelength.
1D Plot Group 2
1 On the Home 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 (sol33).
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.
Point Graph 1
On the 1D Plot Group 2 toolbar, click Point Graph.
1D Plot Group 2
1 In the Settings window for Point Graph, locate the Selection section.
2 Click Paste Selection.
3 In the Paste Selection dialog box, type 2 in the Selection text field.
4 Click OK.
5 In the Settings window for Point Graph, locate the y-Axis Data section.
6 In the Expression text field, type (4*pi*(emw.normEfar^2/E0^2)/
(pi*r0^2)-RCS1)/RCS1.
7 Locate the x-Axis Data section. From the Axis source data list, choose h_size.
8 From the Parameter list, choose Expression.
9 In the Expression text field, type h_size.
10 Click to expand the Coloring and style section. Locate the Coloring and Style section.
In the Width text field, type 2.
11 Find the Line markers subsection. From the Marker list, choose Square.
12 From the Positioning list, choose In data points.
13 Click the x-Axis Log Scale button on the Graphics toolbar.
14 Click the y-Axis Log Scale button on the Graphics toolbar.
15 On the 1D Plot Group 2 toolbar, click Plot.
Compare the convergence plot with that shown in Figure 4.
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
Solved with COMSOL Multiphysics 5.2
The waveguide is made of aluminum. To reduce surface losses, the inside walls are
coated with copper, a high-conductivity metal. The dielectric block is modeled as
having electrical conductivity of σ = 0, relative permeability of μr = 1, and a relative
permittivity of εr = 2.1, with a loss tangent that is a function of temperature,
δ = 0.001(1 + T/300 K). The thermal conductivity of this block is also a function of
temperature, k = 0.3(1 + T/300 K) W/m/K. Furthermore, the density is 2200 kg/m3
and the specific heat is 1050 J/kg/K. These are generic properties representative of a
dielectric material.
At the operating frequency, the skin depth of the copper coating is much smaller than
the dimensions of the waveguide, that is, the electromagnetic fields penetrate a
negligible distance into the walls. This means that the electromagnetic losses can be
localized entirely on the surface, and that there is no need to solve Maxwell’s equations
inside of the walls themselves. Thus, Maxwell’s equations only need to be solved in the
air domain inside of the waveguide, as well as inside of the block. The heat transfer
equation is solved in the block as well as the waveguide walls.
The objective of the analysis is to observe how the assembly of the dielectric block and
waveguide heat up over time, as well as to find the steady-state temperature. The
waveguide is initially assumed to be at a constant temperature throughout. After the
power source is turned on, the electromagnetic fields interact with the highly
conductive interior boundaries of the waveguide, as well as the lossy dielectric block.
The losses in the block and on the walls are sources of heat that raises the temperature.
The block is assumed to be in perfect thermal contact with the walls of the waveguide,
that is, any heat generated in the block is conducted away into the walls. The outside
boundaries of the walls are assumed to be facing ambient air, which leads to free
convective cooling off of these faces. This example uses an averaged heat transfer
coefficient to represent this free convection to ambient air.
The model solves two governing equations: Maxwell’s equations, which describe the
electromagnetic fields, and the heat transfer equation, which describes the
temperature. It is assumed that the operating frequency is much higher than any
thermal transients, and thus it is possible to solve the problem either in a
frequency-transient or a frequency-stationary sense.
2 | RF HEATING
Solved with COMSOL Multiphysics 5.2
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.
Figure 3 plots the fields inside of the waveguide, as well as the temperature of the
assembly, for the steady-state temperature solution after all thermal transients have
died out. The dielectric block shows a significant temperature variation, which affects
the thermal conductivity and loss tangent, plotted in Figure 4.
Figure 2: The maximum temperature, evaluated over the volume of the block, is plotted as
a function of temperature.
3 | RF HEATING
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Heat Transfer>Electromagnetic Heating>Microwave
Heating.
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies for Selected Physics
Interfaces>Frequency-Transient.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
5 | RF HEATING
Solved with COMSOL Multiphysics 5.2
GEOMETRY 1
First, import the geometry of the waveguide including a dielectric block inside the
waveguide.
Import 1 (imp1)
1 On the Home toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 Click Browse.
4 Browse to the application’s Application 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.
2 In the Settings window for Explicit, type Dielectric in the Label text field.
6 | RF HEATING
Solved with COMSOL Multiphysics 5.2
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Air in the Label text field.
7 | RF HEATING
Solved with COMSOL Multiphysics 5.2
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Waveguide in the Label text field.
8 | RF HEATING
Solved with COMSOL Multiphysics 5.2
Explicit 4
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Waveguide inside surfaces in the Label
text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
9 | RF HEATING
Solved with COMSOL Multiphysics 5.2
4 Select Boundaries 16–18, 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, type Waveguide outside surfaces in the
Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
10 | RF HEATING
Solved with COMSOL Multiphysics 5.2
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.
View 2
1 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
2 From the Geometric entity level list, choose Boundary.
11 | RF HEATING
Solved with COMSOL Multiphysics 5.2
12 | RF HEATING
Solved with COMSOL Multiphysics 5.2
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 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.
6 In the Model Builder window, click Heat Transfer in Solids (ht).
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.
13 | RF HEATING
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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 Home 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
Solved with COMSOL Multiphysics 5.2
ADD MATERIAL
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
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, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Dielectric in the Label text field.
3 Select Domain 3 only.
4 Locate the Material Contents section. In the table, enter the following settings:
ADD MATERIAL
1 Go to the Add Material window.
2 In the tree, select Built-In>Copper.
3 Click Add to Component in the window toolbar.
16 | RF HEATING
Solved with COMSOL Multiphysics 5.2
MATERIALS
On the Home 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
Choose the maximum mesh size in the air domain smaller than 0.2 wavelengths using
the parameter h_max that you defined earlier. For the dielectric materials, scale the
mesh size by the inverse of the square root of the relative dielectric constant.
Size 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 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 Air.
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 h_max.
Size 2
1 In the Model Builder window, right-click Mesh 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 h_max/sqrt(2.1).
17 | RF HEATING
Solved with COMSOL Multiphysics 5.2
Free Tetrahedral 1
1 Right-click Mesh 1 and choose Free Tetrahedral.
2 In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Free
Tetrahedral 1 and choose Build All.
STUDY 1
Step 1: Frequency-Transient
1 In the Model Builder window, expand the Study 1 node, then 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.
18 | RF HEATING
Solved with COMSOL Multiphysics 5.2
RESULTS
Data Sets
Plot the transient response of the peak temperature.
Maximum 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.
Point Graph 1
On the 1D Plot Group 1 toolbar, click Point Graph.
1D Plot Group 1
1 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.
2 Locate the y-Axis Data section. From the Unit list, choose degC.
3 On the 1D Plot Group 1 toolbar, click Plot.
Compare the resulting plot with that shown in Figure 2.
Next, add a Frequency-Stationary study to evaluate the peak temperature which can
be observed with the Frequency-Transient study after applying a enough long time so
the peak temperature is saturated.
ADD STUDY
1 On the Home 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 Empty Study.
4 Click Add Study in the window toolbar.
5 On the Home toolbar, click Add Study to close the Add Study window.
19 | RF HEATING
Solved with COMSOL Multiphysics 5.2
STUDY 2
Frequency-Stationary
On the Study toolbar, click Study Steps and choose Stationary>Frequency-Stationary.
Step 1: Frequency-Stationary
1 In the Settings window for Frequency-Stationary, locate the Study Settings section.
2 In the Frequency text field, type f0.
3 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
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Electric>emw.Ex,emw.Ey,emw.Ez -
Electric field.
6 Locate the Arrow Positioning section. Find the x grid points subsection. In the Points
text field, type 40.
7 Find the y grid points subsection. In the Points text field, type 40.
8 Find the z grid points subsection. In the Points text field, type 1.
9 On the Temperature (ht) toolbar, click Plot.
10 Right-click Results>Temperature (ht)>Arrow Volume 1 and choose Duplicate.
11 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Magnetic>emw.Hx,emw.Hy,emw.Hz -
Magnetic field.
12 Locate the Coloring and Style section. From the Color list, choose Green.
13 Right-click Results>Temperature (ht)>Arrow Volume 2 and choose Duplicate.
20 | RF HEATING
Solved with COMSOL Multiphysics 5.2
14 In the Settings window for Arrow Volume, click Replace Expression in the upper-right
corner of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Energy and
power>emw.Poavx,...,emw.Poavz - Power flow, time average.
15 Locate the Coloring and Style section. From the Color list, choose Blue.
Compare the resulting plot with that shown in Figure 3.
Finally, reproduce the plot of the loss tangent on the dielectric block shown in
Figure 4.
To create the plot, reuse the plot group named Isothermal Contours (ht).
21 | RF HEATING
Solved with COMSOL Multiphysics 5.2
22 | RF HEATING
Solved with COMSOL Multiphysics 5.2
Shorted Substrate
Lumped port
Microstrip
line
Metallic post
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.
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
Solved with COMSOL Multiphysics 5.2
Note: In addition to the RF Module, this example requires one of the Acoustics
Module, the MEMS Module, or the Structural Mechanics Module.
Model Definition
This example starts from a basic rectangular cavity filter, whose resonant frequencies
are given by
2 2 2
f nml = ---------------------- -------- + ------- + -----
c mπ nπ lπ
2π ε r μ r a b d
where a and b are the waveguide aperture dimensions and d is the length of the
waveguide cavity. In this example, 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 parts—the cavity walls, post, substrate ground planes, microstrip lines,
and the bottom surface of the piezo device—as 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
Solved with COMSOL Multiphysics 5.2
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 filter’s resonant frequency. The tunable frequency range of this example 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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Structural Mechanics>Piezoelectric Devices.
3 Click Add.
4 In the Select physics tree, select Mathematics>Deformed Mesh>Moving Mesh (ale).
5 Click Add.
6 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
7 Click Add.
8 Click Study.
9 In the Select study tree, select Empty Study.
10 Click Done.
GEOMETRY 1
The model geometry is available as a parameterized geometry sequence in a separate
MPH-file. If you want to build it from scratch, follow the instructions in the section
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
Solved with COMSOL Multiphysics 5.2
Appendix: Geometry Modeling Instructions. Otherwise load it from file with the
following steps.
Cavity (blk1)
1 In the Model Builder window, under Component 1 (comp1)>Geometry 1 right-click
Cavity (blk1) and choose Build All Objects.
2 Click the Wireframe Rendering button on the Graphics toolbar.
GLOBAL DEFINITIONS
Parameters
Add parameters that are not related to geometry.
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
Solved with COMSOL Multiphysics 5.2
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
DEFINITIONS
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
4 Select Boundaries 4, 5, 7, 8, 10, 12, 13, 23, 25, 37, 38, 40, 41, 43–45, 56, and 57
only.
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.
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
Solved with COMSOL Multiphysics 5.2
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.
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
Solved with COMSOL Multiphysics 5.2
Assign material properties. Use three materials for this model: PZT-5H, air, and a
user-defined substrate.
MATERIALS
On the Home toolbar, click Windows and choose Add Material.
ADD MATERIAL
1 Go to the Add Material window.
2 In the tree, select Piezoelectric>Lead Zirconate Titanate (PZT-5H).
3 Click Add to Component in the window toolbar.
MATERIALS
ADD MATERIAL
1 Go to the Add Material window.
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
Solved with COMSOL Multiphysics 5.2
MATERIALS
Air (mat2)
1 On the Home 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, under Component 1 (comp1) 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:
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
Solved with COMSOL Multiphysics 5.2
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)
Next, define the Electrostatics constraints.
1 In the Model Builder window, under Component 1 (comp1) click 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.
4 In the V0 text field, type -V0.
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
Solved with COMSOL Multiphysics 5.2
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.
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
Solved with COMSOL Multiphysics 5.2
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.
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.
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
Solved with COMSOL Multiphysics 5.2
Before generating the mesh, set the physics to control the mesh size.
3 In the Model Builder window, click Electromagnetic Waves, Frequency Domain (emw).
4 In the Settings window for Electromagnetic Waves, Frequency Domain, locate the
Physics-Controlled Mesh section.
5 Select the Enable check box.
6 In the Maximum element size text field, type h_max.
7 Locate the Analysis Methodology section. From the Methodology options list, choose
Fast.
MESH 1
Adjust the minimum element size manually to avoid extremely fine mesh on the piezo
actuator. Then, use swept mesh for the piezo actuator to handle the structural
deformation more efficiently.
Size 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click Size 1.
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
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Size, locate the Element Size Parameters section.
3 In the Minimum element size text field, type 1.
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 1–4 and 6–8 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.
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.
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
Solved with COMSOL Multiphysics 5.2
STUDY 1
Stationary
On the Study toolbar, click Study Steps and choose Stationary>Stationary.
Step 1: Stationary
1 In the Settings window for Stationary, locate the Physics and Variables Selection
section.
2 In the table, enter the following settings:
Frequency Domain
On the Study toolbar, click Study Steps and choose Frequency Domain>Frequency
Domain.
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
Solved with COMSOL Multiphysics 5.2
3 Locate the Physics and Variables Selection section. In the table, enter the following
settings:
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:
RESULTS
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
of the Expression section. From the menu, choose Component 1>Solid
Mechanics>Displacement>solid.disp - Total displacement.
3 On the Displacement toolbar, click Plot.
4 Click the Zoom Extents button on the Graphics toolbar.
This plot shows the deflected piezo actuator disk; compare with Figure 2.
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
Solved with COMSOL Multiphysics 5.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
Solved with COMSOL Multiphysics 5.2
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 S-Parameter (emw) 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.
GLOBAL DEFINITIONS
Load geometrical parameters from a file.
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
Solved with COMSOL Multiphysics 5.2
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 Click Load from File.
4 Browse to the application’s Application Library folder and double-click the file
tunable_cavity_filter_parameters.txt.
GEOMETRY 1
1 In the Model Builder window, under Component 1 (comp1) click Geometry 1.
2 In the Settings window for Geometry, locate the Units section.
3 From the Length unit list, choose mm.
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 and Shape 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 and Shape 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.
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
Solved with COMSOL Multiphysics 5.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 and Shape 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 Right-click Air block and choose Build Selected.
10 Click the Wireframe Rendering button on the Graphics toolbar.
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 and Shape 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.
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
Solved with COMSOL Multiphysics 5.2
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.
Plane Geometry
Add a rectangle for the slot.
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape 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 Right-click Rectangle 1 (r1) and choose Build Selected.
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
Solved with COMSOL Multiphysics 5.2
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.
6 In the z text field, type 0.
7 Click the Build All Objects button.
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 and Shape 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.
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
Solved with COMSOL Multiphysics 5.2
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 Right-click Piezo actuator and choose Build Selected.
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.
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
Solved with COMSOL Multiphysics 5.2
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
Solved with COMSOL Multiphysics 5.2
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 example is to compute the far-field pattern and to compute the impedance of
the structure. Good matching is observed over a wide frequency band.
Substrate with
Shorted 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
1 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
an anechoic 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
Solved with COMSOL Multiphysics 5.2
References
1. D.M. Pozar, Microwave Engineering, John Wiley & Sons, 1998.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
3 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GEOMETRY 1
1 In the Model Builder window, under Component 1 (comp1) click Geometry 1.
2 In the Settings window for Geometry, locate the Units section.
3 From the Length unit list, choose mm.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
GEOMETRY 1
Create a block for the antenna substrate.
Block 1 (blk1)
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 and Shape section. In the Width text field, type 110.
4 In the Depth text field, type 80.
4 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
Block 2 (blk2)
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 and Shape section. 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.
Next, create a work plane where you will draw the Vivaldi antenna pattern. Use two
parametric curves for the tapered slot.
Plane Geometry
1 Click the Zoom Extents button on the Graphics toolbar.
Add a parametric curve using the exponential profile.
Mirror 1 (mir1)
1 On the Work Plane toolbar, click Transforms and choose Mirror.
5 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
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 Right-click Mirror 1 (mir1) and choose Build Selected.
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 and Shape 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.
6 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
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.
5 On the Work Plane toolbar, click Build All.
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, type PML in the Label text field.
3 Locate the Size section. In the Radius text field, type 110.
7 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
4 Click to expand the Layers section. In the table, enter the following settings:
DEFINITIONS
Add a perfectly matched layer.
View 1
Hide some domains to get a better view of the interior parts when setting up the
physics and reviewing the mesh.
8 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
View 1
1 Select Domains 2 and 9 only.
View 1
1 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
2 From the Geometric entity level list, choose Boundary.
9 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
Now set up the physics. Use the selections already defined when assigning boundary
conditions.
10 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
Far-Field Domain 1
On the Physics toolbar, click Domains and choose Far-Field Domain.
Lumped Port 1
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
2 Click the Zoom In button on the Graphics toolbar.
3 Click the Zoom In button on the Graphics toolbar.
4 Select Boundary 20 only.
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 Home 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.
11 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
MATERIALS
On the Home toolbar, click Add Material to close the Add Material window.
Air (mat1)
Override the substrate with a dielectric material of εr = 3.38.
Material 2 (mat2)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Substrate in the Label text field.
3 Select Domains 6 and 7 only.
4 Locate the Material Contents section. In the table, enter the following settings:
MESH 1
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 Element size list, choose Coarse.
12 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
STUDY 1
RESULTS
13 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
Strong electric fields are observed in the slot and microstrip line.
14 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
S-Parameter (emw)
15 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
1D Plot Group 6
On the Home toolbar, click Add Plot Group and choose 1D Plot Group.
16 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
Global 1
On the 1D Plot Group 6 toolbar, click Global.
1D Plot Group 6
1 In the Settings window for Global, click Add Expression in the upper-right corner of
the y-axis data section. From the menu, choose Component 1>Electromagnetic
Waves, Frequency Domain>Ports>emw.VSWR_1 - Voltage standing wave ratio.
2 Click to expand the Legends section. Clear the Show legends check box.
3 On the 1D Plot Group 6 toolbar, click Plot.
This VSWR plot replicates the wide-band frequency response shown in Figure 2.
17 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
18 | VIVALDI ANTENNA
Solved with COMSOL Multiphysics 5.2
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.
The walls of manufactured waveguides are typically plated with a good conductor such
as silver. The model approximates the walls by perfect conductors. This is represented
by the boundary condition n × E = 0.
1 | WAVE G U I D E A D AP T E R
Solved with COMSOL Multiphysics 5.2
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 2
( ν c ) mn = --- ----- + ---
c m n
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 example 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 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:
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
Solved with COMSOL Multiphysics 5.2
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
S 11 = Port 1
----------------------------------------------------------------
-
( 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
Solved with COMSOL Multiphysics 5.2
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
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
4 | WAVE G U I D E AD AP T E R
Solved with COMSOL Multiphysics 5.2
STUDY 1
Finally, add the 3D equation for the propagating wave in the waveguide.
Frequency Domain
On the Study toolbar, click Study Steps and choose Frequency Domain>Frequency
Domain.
5 | WAVE G U I D E A D AP T E R
Solved with COMSOL Multiphysics 5.2
GEOMETRY 1
Import 1 (imp1)
1 On the Home toolbar, click Import.
2 In the Settings window for Import, locate the Import section.
3 Click Browse.
4 Browse to the application’s Application Library folder and double-click the file
waveguide_adapter.mphbin.
5 Click Import.
ADD MATERIAL
1 On the Home 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 Home 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.
Port 1
1 On the Physics toolbar, click Boundaries and choose Port.
2 In the Settings window for Port, locate the Port Properties section.
6 | WAVE G U I D E AD AP T E R
Solved with COMSOL Multiphysics 5.2
The wave enters the adapter through the port with a rectangular cross section.
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.
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.
7 | WAVE G U I D E A D AP T E R
Solved with COMSOL Multiphysics 5.2
2 In the Settings window for Boundary Mode Analysis, locate the Study Settings
section.
3 Select the Search for modes around check box.
4 In the associated text field, type 50.
5 From the Transform list, choose Out-of-plane wave number.
RESULTS
The default plot shows the norm of the electric field on slices through the waveguide; you
can simplify and improve this plot.
8 | WAVE G U I D E AD AP T E R
Solved with COMSOL Multiphysics 5.2
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
of the Expression section. From the menu, choose Model>Component
1>Electromagnetic Waves, Frequency Domain>Electric>Electric field>emw.Ex - Electric
field, x component.
6 Locate the Coloring and Style section. From the Color table list, choose WaveLight.
7 Locate the Plane Data section. In the Planes text field, type 1.
8 On the Electric Field (emw) 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.
S-Parameter (emw)
Select the S-Parameter (emw) plot group under Results in Model Builder. The plot
should closely resemble that in Figure 3.
9 | WAVE G U I D E A D AP T E R
Solved with COMSOL Multiphysics 5.2
10 | WAVE G U I D E A D AP T E R
Solved with COMSOL Multiphysics 5.2
SMA receptacle
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. Except for the microstrip lines and ground plane, model all
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 example analysis, so they are set to PEC by default. The
microstrip lines and ground plane made of 1 oz copper layers are modeled using a
transition boundary condition with 35 μm thickness to address lossy conductive
surfaces due to finite copper conductivity. 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.
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.
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 3D.
2 In the Select physics tree, select Radio Frequency>Electromagnetic Waves, Frequency
Domain (emw).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Frequency Domain.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 On the Home toolbar, click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Here, c_const is a predefined COMSOL constant for the speed of light in vacuum.
GEOMETRY 1
1 In the Model Builder window, under Component 1 (comp1) click Geometry 1.
Block 1 (blk1)
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 and Shape section. 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.
Block 2 (blk2)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Package in the Label text field.
3 Locate the Size and Shape section. 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 Right-click Package and choose Build Selected.
Add a work plane for drawing the layout of the power divider.
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 In the Settings window for Circle, type Ring outer in the Label text field.
3 Click the Zoom Extents button on the Graphics toolbar.
4 Locate the Size and Shape section. In the Radius text field, type r_ring.
Circle 2 (c2)
1 On the Work Plane toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, type Ring inner in the Label text field.
3 Locate the Size and Shape section. In the Radius text field, type r_ring-1.87.
Rectangle 1 (r1)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, type Ring cut in the Label text field.
3 Locate the Size and Shape section. 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.
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.
6 Right-click Difference 1 (dif1) and choose Build Selected.
Rectangle 2 (r2)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, type Lumped element in the Label text field.
3 Locate the Size and Shape section. 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.
Rectangle 3 (r3)
1 On the Work Plane toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type 3.2.
4 In the Height text field, type 5.
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 and Shape 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 and Shape 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.
5 On the Work Plane toolbar, click Build All.
Create the coax SMA receptacle composed of the coax inner and outer conductors, the
SMA connector part, and the flange.
Cylinder 1 (cyl1)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type Coax inner in the Label text field.
3 Locate the Size and Shape section. 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.
Cylinder 2 (cyl2)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type SMA in the Label text field.
3 Locate the Size and Shape section. In the Radius text field, type r_outer+0.6.
4 In the Height text field, type l_sma.
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.
Block 3 (blk3)
1 On the Geometry toolbar, click Block.
2 In the Settings window for Block, type Flange in the Label text field.
3 Locate the Size and Shape section. 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.
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, type Coax dielectric in the Label text field.
3 Locate the Size and Shape section. In the Radius text field, type r_outer.
4 In the Height text field, type l_sma.
5 Locate the Position section. In the x text field, type -7.
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.
2 Select the objects copy1(1), copy1(3), and copy1(5) only.
3 In the Settings window for Rotate, locate the Rotation Angle section.
4 In the Rotation text field, type 180.
Cylinder 4 (cyl4)
1 On the Geometry toolbar, click Cylinder.
2 In the Settings window for Cylinder, type Screw in the Label text field.
3 Locate the Size and Shape section. 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.
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.
8 Click the Build All Objects button.
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, type Microstrip line in the Label text field.
3 Locate the Input Entities section. From the Geometric entity level list, choose
Boundary.
4 Select Boundaries 6, 54, 61, 102, and 147 only.
Explicit 2
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Substrate in the Label text field.
Explicit 3
1 On the Definitions toolbar, click Explicit.
2 In the Settings window for Explicit, type Coax dielectric in the Label text field.
3 Select Domains 9, 14, and 19 only.
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.
2 In the Settings window for Explicit, type Metal volume in the Label text field.
3 Select Domains 4–8, 10, 12, 13, 16–18, 20, and 22–26 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, type Model domain in the Label text field.
3 Locate the Input Entities section. Under Selections to invert, click Add.
4 In the Add dialog box, select Metal volume in the Selections to invert list.
5 In the Selections to invert list, select Metal volume.
6 Click OK.
View 1
Suppress some boundaries to get a view of the interior while setting the physics and
mesh.
View 1
1 In the Settings window for Hide Geometric Entities, locate the Geometric Entity
Selection section.
screws. Add a Transition Boundary Condition to the microstrip line and the substrate
ground plane.
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.
Lumped Port 2
1 On the Physics toolbar, click Boundaries and choose Lumped Port.
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.
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.
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 Home 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.
5 Go to the Add Material window.
6 In the tree, select Built-In>Copper.
7 Click Add to Component in the window toolbar.
MATERIALS
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 Microstrip line.
5 On the Home 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 3 (mat3)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Substrate in the Label text field.
3 Locate the Geometric Entity Selection section. From the Selection list, choose
Substrate.
4 Locate the Material Contents section. In the table, enter the following settings:
Material 4 (mat4)
1 Right-click Materials and choose Blank Material.
2 In the Settings window for Material, type PTFE in the Label text field.
3 Locate the Geometric Entity Selection section. From the Selection list, choose Coax
dielectric.
4 Locate the Material Contents section. In the table, enter the following settings:
MESH 1
1 In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and
choose Build All.
STUDY 1
RESULTS
1 In the Settings window for 3D Plot Group, locate the Data section.
2 From the Parameter value (freq (Hz)) list, choose 3E9.
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.