RF Module Users Guide
RF Module Users Guide
User’s Guide
RF Module User’s Guide
© 1998–2020 COMSOL
Protected by patents listed on www.comsol.com/patents, and U.S. Patents 7,519,518; 7,596,474;
7,623,991; 8,457,932; 9,098,106; 9,146,652; 9,323,503; 9,372,673; 9,454,625; 10,019,544;
10,650,177; and 10,776,541. Patents pending.
This Documentation and the Programs described herein are furnished under the COMSOL Software License
Agreement (www.comsol.com/comsol-license-agreement) and may be used or copied only under the terms
of the license agreement.
COMSOL, the COMSOL logo, COMSOL Multiphysics, COMSOL Desktop, COMSOL Compiler,
COMSOL Server, and LiveLink are either registered trademarks or trademarks of COMSOL AB. All other
trademarks are the property of their respective owners, and COMSOL AB and its subsidiaries and products
are not affiliated with, endorsed by, sponsored by, or supported by those trademark owners. For a list of such
trademark owners, see www.comsol.com/trademarks.
Version: COMSOL 5.6
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:
Chapter 1: Introduction
Chapter 2: RF Modeling
Simplifying Geometries 28
2D Models . . . . . . . . . . . . . . . . . . . . . . . . . 28
3D Models . . . . . . . . . . . . . . . . . . . . . . . . . 30
Using Efficient Boundary Conditions . . . . . . . . . . . . . . . 31
Applying Electromagnetic Sources . . . . . . . . . . . . . . . . 31
Meshing and Solving . . . . . . . . . . . . . . . . . . . . . . 32
CONTENTS |3
S-Parameters and Ports 45
S-Parameters in Terms of Voltage and Electric Field . . . . . . . . . 45
S-Parameter Calculations . . . . . . . . . . . . . . . . . . . . 46
S-Parameter Variables . . . . . . . . . . . . . . . . . . . . . 46
Port Sweeps/Manual Terminal Sweeps and Touchstone Export . . . . . 46
4 | CONTENTS
Part Libraries 70
Material Libraries 71
Maxwell’s Equations 74
Introduction to Maxwell’s Equations . . . . . . . . . . . . . . . 74
Constitutive Relations . . . . . . . . . . . . . . . . . . . . . 75
Potentials. . . . . . . . . . . . . . . . . . . . . . . . . . 77
Electromagnetic Energy . . . . . . . . . . . . . . . . . . . . 77
Material Properties . . . . . . . . . . . . . . . . . . . . . . 78
About the RF Materials Database. . . . . . . . . . . . . . . . . 80
Boundary and Interface Conditions . . . . . . . . . . . . . . . . 80
Phasors . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Special Calculations 83
S-Parameter Calculations . . . . . . . . . . . . . . . . . . . . 83
Far-Field Calculations Theory . . . . . . . . . . . . . . . . . . 86
References . . . . . . . . . . . . . . . . . . . . . . . . . 87
Electromagnetic Quantities 88
CONTENTS |5
Perfect Electric Conductor . . . . . . . . . . . . . . . . . . 111
Perfect Magnetic Conductor . . . . . . . . . . . . . . . . . 112
Port . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
Integration Line for Current . . . . . . . . . . . . . . . . . 122
Integration Line for Voltage . . . . . . . . . . . . . . . . . . 122
Ground . . . . . . . . . . . . . . . . . . . . . . . . . 122
Electric Potential . . . . . . . . . . . . . . . . . . . . . . 122
Circular Port Reference Axis . . . . . . . . . . . . . . . . . 123
Diffraction Order . . . . . . . . . . . . . . . . . . . . . 123
Orthogonal Polarization . . . . . . . . . . . . . . . . . . . 125
Periodic Port Reference Point . . . . . . . . . . . . . . . . . 126
Lumped Port . . . . . . . . . . . . . . . . . . . . . . . 127
Lumped Element . . . . . . . . . . . . . . . . . . . . . . 130
Uniform Element . . . . . . . . . . . . . . . . . . . . . . 131
Electric Field . . . . . . . . . . . . . . . . . . . . . . . 131
Magnetic Field . . . . . . . . . . . . . . . . . . . . . . . 132
Scattering Boundary Condition . . . . . . . . . . . . . . . . 132
Reference Point . . . . . . . . . . . . . . . . . . . . . . 136
Impedance Boundary Condition . . . . . . . . . . . . . . . . 137
Surface Current Density . . . . . . . . . . . . . . . . . . . 139
Surface Magnetic Current Density . . . . . . . . . . . . . . . 139
Surface Roughness . . . . . . . . . . . . . . . . . . . . . 140
Transition Boundary Condition . . . . . . . . . . . . . . . . 141
Periodic Condition . . . . . . . . . . . . . . . . . . . . . 142
Magnetic Current . . . . . . . . . . . . . . . . . . . . . 144
Two-Port Network . . . . . . . . . . . . . . . . . . . . . 144
Two-Port Network Port . . . . . . . . . . . . . . . . . . . 145
Mixed Mode S-Parameters . . . . . . . . . . . . . . . . . . 145
Edge Current . . . . . . . . . . . . . . . . . . . . . . . 146
Electric Point Dipole . . . . . . . . . . . . . . . . . . . . 146
Magnetic Point Dipole . . . . . . . . . . . . . . . . . . . . 146
Line Current (Out-of-Plane) . . . . . . . . . . . . . . . . . 147
Archie’s Law . . . . . . . . . . . . . . . . . . . . . . . . 147
Effective Medium . . . . . . . . . . . . . . . . . . . . . . 148
6 | CONTENTS
Wave Equation, Electric . . . . . . . . . . . . . . . . . . . 153
Initial Values . . . . . . . . . . . . . . . . . . . . . . . 156
Drude–Lorentz Polarization. . . . . . . . . . . . . . . . . . 156
CONTENTS |7
Theory for the Electromagnetic Waves Interfaces 180
Introduction to the Physics Interface Equations . . . . . . . . . . 180
Frequency Domain Equation . . . . . . . . . . . . . . . . . 181
Time Domain Equation . . . . . . . . . . . . . . . . . . . 186
Vector Elements . . . . . . . . . . . . . . . . . . . . . . 188
Eigenfrequency Calculations. . . . . . . . . . . . . . . . . . 189
Gaussian Beams as Background Fields and Input Fields. . . . . . . . 189
Effective Material Properties in Effective Media and Mixtures . . . . . 191
Effective Conductivity in Effective Media and Mixtures . . . . . . . . 192
Effective Relative Permittivity in Effective Media and Mixtures . . . . . 193
Effective Relative Permeability in Effective Media and Mixtures . . . . 194
Archie’s Law Theory . . . . . . . . . . . . . . . . . . . . 195
Reference for Archie’s Law . . . . . . . . . . . . . . . . . . 195
8 | CONTENTS
Voltage-Controlled Voltage Source . . . . . . . . . . . . . . . 215
Voltage-Controlled Current Source. . . . . . . . . . . . . . . 216
Current-Controlled Voltage Source. . . . . . . . . . . . . . . 216
Current-Controlled Current Source . . . . . . . . . . . . . . 217
Switch . . . . . . . . . . . . . . . . . . . . . . . . . . 217
Subcircuit Definition . . . . . . . . . . . . . . . . . . . . 218
Subcircuit Instance . . . . . . . . . . . . . . . . . . . . . 219
NPN BJT and PNP BJT. . . . . . . . . . . . . . . . . . . . 219
n-Channel MOSFET and p-Channel MOSFET . . . . . . . . . . . 220
Mutual Inductance . . . . . . . . . . . . . . . . . . . . . 220
Transformer . . . . . . . . . . . . . . . . . . . . . . . 221
Diode . . . . . . . . . . . . . . . . . . . . . . . . . . 221
External I vs. U . . . . . . . . . . . . . . . . . . . . . . 222
External U vs. I . . . . . . . . . . . . . . . . . . . . . . 223
External I-Terminal . . . . . . . . . . . . . . . . . . . . . 224
SPICE Circuit Import . . . . . . . . . . . . . . . . . . . . 225
SPICE Circuit Export . . . . . . . . . . . . . . . . . . . . 225
Chapter 7: Glossary
CONTENTS |9
10 | C O N T E N T S
1
Introduction
This guide describes the RF Module, an optional add-on package for COMSOL
Multiphysics® with customized physics interfaces and functionality optimized for
the analysis of electromagnetic waves.
This chapter introduces you to the capabilities of this module. A summary of the
physics interfaces and where you can find documentation and model examples is
also included. The last section is a brief overview with links to each chapter in this
guide.
11
About the RF Module
In this section:
The module is useful for component design in virtually all areas where you find
electromagnetic waves, such as:
• Antennas
• Filters, couplers, and power dividers
• Planar circuits and passive devices
• RF interconnects and packages
• Waveguides and cavity resonators
• Frequency-selective surfaces
• Metamaterials
12 | CHAPTER 1: INTRODUCTION
The physics interfaces cover the following types of electromagnetics field simulations
and handle time-harmonic, time-dependent, and eigenfrequency/eigenmode
problems:
Material properties include inhomogeneous and fully anisotropic materials, media with
gains or losses, and complex-valued material properties. In addition to the standard
postprocessing features, the module supports direct computation of S-parameters and
far-field radiation patterns. You can add ports with a wave excitation with specified
power level and mode type, and add PMLs (perfectly matched layers) to simulate
electromagnetic waves that propagate into an unbounded domain. For time-harmonic
simulations, you can use the scattered wave or the total wave.
This module also has interfaces for circuit modeling, a SPICE interface, and support
for importing ECAD drawings.
Independently of the structure size, the module accommodates any case of nonlinear,
inhomogeneous, or anisotropic media. It also handles materials with properties that
vary as a function of time as well as frequency-dispersive materials.
14 | CHAPTER 1: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
AC/DC
Electromagnetic Heating
16 | CHAPTER 1: INTRODUCTION
COMPARING THE TIME DEPENDENT AND FREQUENCY DOMAIN STUDIES
When variations in time are present there are two main approaches to represent the
time dependence. The most straightforward is to solve the problem by calculating the
changes in the solution for each time step; that is, solving using the Time Dependent
study (available with the Electromagnetic Waves, Transient interface). However, this
approach can be time consuming if small time steps are necessary for the desired
accuracy. It is necessary when the inputs are transients like turn-on and turn-off
sequences.
However, if the Frequency Domain study available with the Electromagnetic Waves,
Frequency Domain interface is used, this allows you to efficiently simplify and assume
that all variations in time occur as sinusoidal signals. Then the problem is
time-harmonic and in the frequency domain. Thus you can formulate it as a stationary
problem with complex-valued solutions. The complex value represents both the
amplitude and the phase of the field, while the frequency is specified as a scalar model
input, usually provided by the solver. This approach is useful because, combined with
Fourier analysis, it applies to all periodic signals with the exception of nonlinear
problems. Examples of typical frequency domain simulations are wave-propagation
problems like waveguides and antennas.
For nonlinear problems you can apply a Frequency Domain study after a linearization
of the problem, which assumes that the distortion of the sinusoidal signal is small.
Use a Time Dependent study when the nonlinear influence is strong, or if you are
interested in the harmonic distortion of a sine signal. It can also be more efficient to
use a Time Dependent study if you have a periodic input with many harmonics, like a
square-shaped signal.
DEFAULT FREQUENCY
In the RF module, the Frequency Domain, the Frequency-Domain Modal, and the
Eigenfrequency study step have GHz as the default frequency unit and a default
frequency of 1[GHz].
Even after a model is defined, you can edit to input data, equations, boundary
conditions, geometry — the equations and boundary conditions are still available
through associative geometry — and mesh settings. You can restart the solver, for
example, using the existing solution as the initial condition or initial guess. It is also
easy to add another physics interface to account for a phenomenon not previously
described in a model.
18 | CHAPTER 1: INTRODUCTION
context-based) help, and the application libraries are all accessed through the
COMSOL Desktop.
• In the Model Builder or Physics Builder click a node or window and then
press F1.
• On the main toolbar, click the Help ( ) button.
• From the main menu, select Help>Help.
• Press Ctrl+F1.
• From the File menu select Help>Documentation ( ).
• Press Ctrl+F1.
• On the main toolbar, click the Documentation ( ) button.
• From the main menu, select Help>Documentation.
Once the Application Libraries window is opened, you can search by name or browse
under a module folder name. Click to view a summary of the model or application and
its properties, including options to open it or its associated PDF document.
20 | CHAPTER 1: INTRODUCTION
Opening the Application Libraries Window
To open the Application Libraries window ( ):
To include the latest versions of model examples, from the Help menu
select ( ) Update COMSOL Application Library.
22 | CHAPTER 1: INTRODUCTION
Overview of the User’s Guide
The RF Module User’s Guide gets you started with modeling using COMSOL
Multiphysics. The information in this guide is specific to this module. Instructions how
to use COMSOL in general are included with the COMSOL Multiphysics Reference
Manual.
RF THEORY
The Electromagnetics Theory chapter contains a review of the basic theory of
electromagnetics, starting with Maxwell’s Equations, and the theory for some Special
Calculations: S-parameters, lumped port parameters, and far-field analysis. There is
also a list of Electromagnetic Quantities with their SI units and symbols.
RADIO FREQUENCY
Radio Frequency Interfaces chapter describes:
ELECTRICAL CIRCUIT
AC/DC Interfaces chapter describes The Electrical Circuit Interface, which simulates
the current in a conductive and capacitive material under the influence of an electric
field. All three study types (Stationary, Frequency Domain, and Time Dependent) are
available. The underlying theory is also included at the end of the chapter.
HEAT TRANSFER
Heat Transfer Interfaces chapter describes the Microwave Heating interface, which
combines the physics features of an Electromagnetic Waves, Frequency Domain
interface from the RF Module with the Heat Transfer interface. The predefined
interaction adds the electromagnetic losses from the electromagnetic waves as a heat
source and solves frequency domain (time-harmonic) electromagnetic waves in
conjunction with stationary or transient heat transfer. This physics interface is based on
the assumption that the electromagnetic cycle time is short compared to the thermal
time scale (adiabatic assumption). The underlying theory is also included at the end of
the chapter.
24 | CHAPTER 1: INTRODUCTION
2
RF Modeling
The goal of this chapter is to familiarize you with the modeling procedure in the
RF Module. A number of models available through the RF Module Applications
Libraries also illustrate the different aspects of the simulation process.
In this chapter:
25
• Reduced Order Modeling
• Part Libraries and Material Libraries
26 | CHAPTER 2: RF MODELING
P r e pa r i ng f o r R F Mod el i n g
Several modeling topics are described in this section that might not be found in
ordinary textbooks on electromagnetic theory.
Increasing the complexity of a model to make it more accurate usually makes it more
expensive to simulate. A complex model is also more difficult to manage and interpret
than a simple one. Keep in mind that it can be more accurate and efficient to use several
simple models instead of a single, complex one.
In this section:
• 2D Models
• 3D Models
• Using Efficient Boundary Conditions
• Applying Electromagnetic Sources
• Meshing and Solving
2D Models
The text below is a guide to some of the common approximations made for 2D
models. Remember that the modeling in 2D usually represents some 3D geometry
under the assumption that nothing changes in the third dimension or that the field has
a prescribed propagation component in the third dimension.
CARTESIAN COORDINATES
In this case a cross section is viewed in the xy-plane of the actual 3D geometry. The
geometry is mathematically extended to infinity in both directions along the z-axis,
assuming no variation along that axis or that the field has a prescribed wave vector
component along that axis. All the total flows in and out of boundaries are per unit
length along the z-axis. A simplified way of looking at this is to assume that the
geometry is extruded one unit length from the cross section along the z-axis. The total
flow out of each boundary is then from the face created by the extruded boundary (a
boundary in 2D is a line).
28 | CHAPTER 2: RF MODELING
There are usually two approaches that lead to a 2D cross-section view of a problem.
The first approach is when it is known that there is no variation of the solution in one
particular dimension.
This is shown in the model H-Bend Waveguide 2D, where the electric field only has
one component in the z direction and is constant along that axis. The second approach
is when there is a problem where the influence of the finite extension in the third
dimension can be neglected.
Figure 2-1: The cross sections and their real geometry for Cartesian coordinates and
cylindrical coordinates (axial symmetry).
SIMPLIFYING GEOMETRIES | 29
When using the axisymmetric versions, the horizontal axis represents the
radial (r) direction and the vertical axis the z direction, and the geometry
in the right half plane (that is, for positive r only) must be created.
POLARIZATION IN 2D
In addition to selecting 2D or 2D axisymmetry when you start building the model, the
physics interfaces (The Electromagnetic Waves, Frequency Domain Interface or The
Electromagnetic Waves, Transient Interface) in the Model Builder offers a choice in
the Components settings section. The available choices are Out-of-plane vector,
In-plane vector, and Three-component vector. This choice determines what
polarizations can be handled. For example, as you are solving for the electric field, a
2D TM (out-of-plane H field) model requires choosing In-plane vector as then the
electric field components are in the modeling plane.
3D Models
Although COMSOL Multiphysics fully supports arbitrary 3D geometries, it is
important to simplify the problem. This is because 3D models often require more
computer power, memory, and time to solve. The extra time spent on simplifying a
model is probably well spent when solving it. Below are a few issues that need to be
addressed before starting to implement a 3D model in this module.
• Check if it is possible to solve the problem in 2D. Given that the necessary
approximations are small, the solution is more accurate in 2D, because a much
denser mesh can be used.
• Look for symmetries in the geometry and model. Many problems have planes where
the solution is the same on both sides of the plane. A good way to check this is to
flip the geometry around the plane, for example, by turning it up-side down around
the horizontal plane. Then remove the geometry below the plane if no differences
are observed between the two cases regarding geometry, materials, and sources.
Boundaries created by the cross section between the geometry and this plane need
a symmetry boundary condition, which is available in all 3D physics interfaces.
• There are also cases when the dependence along one direction is known, and it can
be replaced by an analytical function. Use this approach either to convert 3D to 2D
or to convert a layer to a boundary condition.
30 | CHAPTER 2: RF MODELING
Using Efficient Boundary Conditions
An important technique to minimize the problem size is to use efficient boundary
conditions. Truncating the geometry without introducing too large errors is one of the
great challenges in modeling. Below are a few suggestions of how to do this. They
apply to both 2D and 3D problems.
• Many models extend to infinity or can have regions where the solution only
undergoes small changes. This problem is addressed in two related steps. First, the
geometry needs to be truncated in a suitable position. Second, a suitable boundary
condition needs to be applied there. For static and quasistatic models, it is often
possible to assume zero fields at the open boundary, provided that this is at a
sufficient distance away from the sources. For radiation problems, special
low-reflecting boundary conditions need to be applied. This boundary should be in
the order of a few wavelengths away from any source.
A more accurate option is to use perfectly matched layers (PMLs). PMLs are layers
that absorbs all radiated waves with small reflections.
• Replace thin layers with boundary conditions where possible. There are several types
of boundary conditions in COMSOL Multiphysics suitable for such replacements.
For example, replace materials with high conductivity by the perfect electric
conductor (PEC) boundary condition.
• Use boundary conditions for known solutions. For example, an antenna aperture
can be modeled as an equivalent surface current density on a 2D face (boundary) in
a 3D model.
SIMPLIFYING GEOMETRIES | 31
Meshing and Solving
The finite element method approximates the solution within each element, using some
elementary shape function that can be constant, linear, or of higher order. Depending
on the element order in the model, a finer or coarser mesh is required to resolve the
solution. In general, there are three problem-dependent factors that determine the
necessary mesh resolution:
• The first is the variation in the solution due to geometrical factors. The mesh
generator automatically generates a finer mesh where there is a lot of fine
geometrical details. Try to remove such details if they do not influence the solution,
because they produce a lot of unnecessary mesh elements.
• The second is the skin effect or the field variation due to losses. It is easy to estimate
the skin depth from the conductivity, permeability, and frequency. At least two linear
elements per skin depth are required to capture the variation of the fields. If the skin
depth is not studied or a very accurate measure of the dissipation loss profile is not
needed, replace regions with a small skin depth with a boundary condition, thereby
saving elements. If it is necessary to resolve the skin depth, the boundary layer
meshing technique can be a convenient way to get a dense mesh near a boundary.
• The third and last factor is the wavelength. To resolve a wave properly, it is necessary
to use about 10 linear (or five 2nd order) elements per wavelength. Keep in mind
that the wavelength depends on the local material properties.
SOLVERS
In most cases the solver sequence generated by COMSOL Multiphysics can be used.
The choice of solver is optimized for the typical case for each physics interface and
study type in this module. However, in special cases tuning the solver settings can be
required. This is especially important for 3D problems because they can require a large
amount of memory.
• Meshing
• Studies and Solvers
32 | CHAPTER 2: RF MODELING
Periodic Boundary Conditions
The RF Module has a dedicated Periodic Condition. The periodic condition can
identify simple mappings on plane source and destination boundaries of equal shape.
The destination can also be rotated with respect to the source. There are three types
of periodic conditions available (only the first two for transient analysis):
• Continuity — The tangential components of the solution variables are equal on the
source and destination.
• Antiperiodicity — The tangential components have opposite signs.
• Floquet periodicity — There is a phase shift between the tangential components.
The phase shift is determined by a wave vector and the distance between the source
and destination. Floquet periodicity is typically used for models involving plane
waves interacting with periodic structures.
Periodic boundary conditions must have compatible meshes. This can be done
automatically by enabling the Physics-control mesh in the setting for The
Electromagnetic Waves, Frequency Domain Interface or by manually setting up the
correct mesh sequence
To learn how to use the Copy Mesh feature to ensure that the mesh on
the destination boundary is identical to that on the source boundary, see
Plasmonic Wire Grating: Application Library path RF_Module/Tutorials/
plasmonic_wire_grating.
34 | CHAPTER 2: RF MODELING
Scattered Field Formulation
For problems where a known background field is illuminating an object in free space
it is possible to use the scattered field formulation. Since the equation of the
background field is known it can be entered as a model input and does not need to be
computed. Starting from the frequency-domain governing equation:
∇ × ( μ –1 ∇ × E ) – ω 2 εc E = 0
The known background field becomes a source term and the scattered field
formulation thus solves for the relative electric field. A linearly polarized plane wave
background field, a paraxial-approximate Gaussian beam, or a user-defined
background field can be specified. When solving the scattered field formulation the
total, the background, and the relative electric fields are available. The relative field is
the difference between the background field and the total field. It is the relative field
that contributes to the far-field calculation. For more information about the Far-Field
computation, see Far-Field Calculations Theory. The benefit to this approach is that if
the background field is much larger in magnitude than the scattered field, the accuracy
of the simulation improves if the relative field is solved for. Another advantage is that
is becomes very easy to set up a perfectly matched layer surrounding the homogeneous
medium modeling domain.
The drawback to this approach is that the relative field requires some careful
interpretation. The relative electric field can conceptually be decomposed into:
The Escattered component is the scattered field from object. This is the field that is of
interest in a scattering problem. However, the relative field may also consist of a
component that represents a correction to the background field and a cancellation of
the background field. The Ecorrection component can be nonzero when the
background field does not exactly satisfy Maxwell’s equations, such as when the
paraxial Gaussian beam approximation is used for a tightly focused beam. For more
information about the Gaussian beam theory, see Gaussian Beams as Background
Fields and Input Fields. The Ecancellation component will be nonzero and equal to
−Ebackground wherever the total field should be zero, such as in the interior of any
An alternative of using the scattered-field formulation, is to use ports with the Activate
slit condition on interior port setting enabled. Then the domain can be excited by the
interior port and the outgoing field can be absorbed by perfectly matched layers. For
more information about the Port feature and the Activate slit condition on interior port
setting, see Port Properties.
36 | CHAPTER 2: RF MODELING
Modeling with Far-Field Calculations
The far electromagnetic field from, for example, antennas can be calculated from the
near-field solution on a boundary using far-field analysis. The antenna is located in the
vicinity of the origin, while the far-field is taken at infinity but with a well-defined
angular direction ( θ, ϕ ) . The far-field radiation pattern is given by evaluating the
squared norm of the far-field on a sphere centered at the origin. Each coordinate on
the surface of the sphere represents an angular direction.
In this section:
The Far-Field Domain and the Far-Field Calculation nodes get their selections
automatically, if the Perfectly Matched Layer (PML) feature has been defined before
adding the Far-Field Domain feature.
For each variable name entered, the software generates functions and variables, which
represent the vector components of the far electric field. The names of these variables
For example, the name Efar is entered and the geometry is Cartesian with the
independent variables x, y, and z, the generated variables get the names Efarx, Efary,
and Efarz.
If, on the other hand, the geometry is axisymmetric with the independent variables r,
phi, and z, the generated variables get the names Efarr, Efarphi, and Efarz.
In 2D, the software only generates the variables for the nonzero field components. The
physics interface name also appears in front of the variable names so they can vary, but
typically look something like emw.Efarz and so forth.
To each of the generated variables, there is a corresponding function with the same
name. This function takes the vector components of the evaluated far-field direction as
arguments.
The vector components also can be interpreted as a position. For example, assume that
the variables dx, dy, and dz represent the direction in which the far electric field is
evaluated.
The expression
Efarx(dx,dy,dz)
gives the value of the far electric field in this direction. To give the direction as an angle,
use the expression
Efarx(sin(theta)*cos(phi),sin(theta)*sin(phi),cos(theta))
where the variables theta and phi are defined to represent the angular direction
( θ, ϕ ) in radians. The magnitude of the far field and its value in dB are also generated
as the variables normEfar and normdBEfar, respectively.
38 | CHAPTER 2: RF MODELING
The variables are plotted for a selected number of angles on a unit circle (in 2D) or a
unit sphere (in 3D). The angle interval and the number of angles can be manually
specified. For 2D Radiation Pattern plots also the reference direction from which the
angle is measured and the normal to the plane the far field is computed for can be
specified. For 3D Radiation Pattern plots you also specify an expression for the surface
color.
The main advantage with the Radiation Pattern plot, as compared to making a Line
Graph, is that the unit circle/sphere that you use for defining the plot directions, is not
part of your geometry for the solution. Thus, the number of plotting directions is
decoupled from the discretization of the solution domain.
4
See Far-Field Analysis Using Functions and operators.
40 | CHAPTER 2: RF MODELING
FAR-FIELD ANALYSIS USING FUNCTIONS AND OPERATORS
The postprocessing far-field functions are available under
Component>Definitions>Functions. Below you find example models using these
functions and some links to more information.
• Far-field analysis using circular port excitation with a positive azimuthal mode
number
• Scattered field analysis excited by the predefined circularly polarized plane wave type
The function can be used in a 3D Radiation Pattern plot, where the input argument of
the function must be same as the Azimuth angle variable in the Evaluation section in the
settings window.
The suffix of a function name varies based on the circular port mode type, port mode
number and azimuthal mode number in the physics interface. For example, when
using azimuthal mode number 1 in the physics interface and transverse electric (TE)
When the function is used in a radiation pattern plot under a 1D or a polar plot group,
the value of input argument defines the plotting plane regardless of the normal and
reference direction in the evaluation section in the settings window. For example,
norm3DEfar_TE12(0)evaluates the norm of the electric far field for the TE12 mode
for 0-degree azimuthal angle. This is equivalent to plotting this variable on the
xz-plane. Similarly, norm3DEfar_TE12(pi/2) is the evaluation at 90-degree azimuthal
angle, which is equivalent to plotting the variable on the yz-plane.
The 3D far-field norm, the linear superposition of the positive and negative azimuthal
modes scaled by 0.5, is
2 2 2
E r cos mφ + E φ sin m φ + E z cos mφ ,
nx ny
sin ----- ( 2πd x sin θ cos φ + α x ) sin ----- ( 2πsd y sin θ cos φ + α y )
2 2
----------------------------------------------------------------------------- ´ --------------------------------------------------------------------------------
2πd sin θ cos φ + α 2πd y sin θ cos φ + α y
sin --------------------------------------------------- sin ---------------------------------------------------
x x
2 2
nz
sin ----- ( 2πd z sin θ cos φ + α z )
2
´ ----------------------------------------------------------------------------- ,
2πd z sin θ cos φ + α z
sin ---------------------------------------------------
2
The uniform two dimensional array factor operator af2 is simpler than the three
dimensional version, as the third, the z-component factor, is unity.
42 | CHAPTER 2: RF MODELING
The number of input arguments for the array factor operators depends on the
dimension of model component, 2D, 2D Axisymmetric, or 3D.
TABLE 2-2: INPUT ARGUMENTS OF ARRAY FACTOR OPERATOR
Gain
The antenna realized gain is defined as
2
normEfar
G = ----------------------------------
60P delivered
2
where the delivered power, Pdelivered is P in ( 1 – S 11 ) . The gain is available only
when the S-parameter calculation is valid, that is, for the single port excitation case.
44 | CHAPTER 2: RF MODELING
S-Parameters and Ports
In this section:
S 11 S 12 . . S 1n
S 21 S 22 . . .
S = . . . . .
. . . . .
S n1 . . . S nn
where S11 is the voltage reflection coefficient at port 1, S21 is the voltage transmission
coefficient from port 1 to port 2, and so on. The time average power reflection/
transmission coefficients are obtained as |Sij |2.
• For more details about lumped ports, see Lumped Ports with Voltage
Input.
• See Port and Lumped Port for instructions to set up a model.
S-Parameter Variables
This module automatically generates variables for the S-parameters. The port names
(use numbers for sweeps to work correctly) determine the variable names. If, for
example, there are two ports with the numbers 1 and 2 and Port 1 is the inport, the
software generates the variables S11 and S21. S11 is the S-parameter for the reflected
wave and S21 is the S-parameter for the transmitted wave. For convenience, two
variables for the S-parameters on a dB scale, S11dB and S21dB, are also defined using
the following relation:
S 11dB = 20 log 10 ( S 11 )
The model and physics interface names also appear in front of the variable names so
they can vary. The S-parameter variables are added to the predefined quantities in
appropriate plot lists.
46 | CHAPTER 2: RF MODELING
The Frequency Domain Source Sweep study is another way of making efficient port
sweeps. It is available as a preset study for the Electromagnetic Waves, Frequency
Domain interface.
Exporting a Touchstone file can also be performed by right-clicking the Export node
under Results and selecting Touchstone or by selecting Touchstone under Data in the
Results ribbon toolbar (Windows users) or the Results context menu (Mac and Linux
users).
A lumped port specified as an input port calculates the impedance, Zport, and S11
S-parameter for that port. The parameters are directly given by the relations
V port
Z port = -------------
I port
V port – V in
S 11 = ----------------------------
V in
where Vport is the extracted voltage for the port given by the electric field line integral
between the terminals averaged over the entire port. The current Iport is the averaged
total current over all cross sections parallel to the terminals. Ports not specified as input
ports only return the extracted voltage and current.
48 | CHAPTER 2: RF MODELING
Lumped Port Parameters
In transmission line theory voltages and currents are dealt with rather than electric and
magnetic fields, so the lumped port provides an interface between them. The
requirement on a lumped port is that the feed point must be similar to a transmission
line feed, so its gap must be much less than the wavelength. It is then possible to define
the electric field from the voltage as
V = E ⋅ dl = ( E ⋅ ah ) dl
h h
where h is a line between the terminals at the beginning of the transmission line, and
the integration is going from positive (phase) V to ground. The current is positive
going into the terminal at positive V.
I
+V
Js E h
Ground
n
Lumped port boundary
The transmission line current can be represented with a surface current at the lumped
port boundary directed opposite to the electric field.
V
Z = ----
I
and in analogy to this an equivalent surface impedance is defined at the lumped port
boundary
E ⋅ ah
η = -------------------------
Js ⋅ ( –ah )
To calculate the surface current density from the current, integrate along the width, w,
of the transmission line
I = ( n × Js ) ⋅ dl = – ( Js ⋅ ah ) dl
w w
( E ⋅ a h ) dl ( E ⋅ a h ) dl
V h
Z = ---- = ----------------------------------- = η ------------------------------ ≈ η ----
h h
I w
– ( J s ⋅ a h ) dl ( E ⋅ a h ) dl
w w
w
η = Z ----
h
where the last approximation assumed that the electric field is constant over the
integrations. A similar relationship can be derived for coaxial cables
2π
η = Z ----------
b
ln ---
a
The transfer equations above are used in an impedance type boundary condition,
relating surface current density to tangential electric field via the surface impedance.
1 1
n × ( H 1 – H 2 ) + --- n × ( E × n ) = 2 --- n × ( E 0 × n )
η η
where E is the total field and E0 the incident field, corresponding to the total voltage,
V, and incident voltage, V0, at the port.
When using the lumped port as a circuit port, the port voltage is fed as
input to the circuit and the current computed by the circuit is applied as
a uniform current density, that is as a surface current condition. Thus, an
open (unconnected) circuit port is just a continuity condition.
50 | CHAPTER 2: RF MODELING
LUMPED PORT VARIABLES
Each lumped port generates variables that are accessible to the user. Apart from the
S-parameter, a lumped port condition also generates the following variables.
NAME DESCRIPTION
For example, a lumped port with port number 1, defined in the first geometry, for the
Electromagnetic Waves interface with the tag emw, defines the port impedance variable
emw.Zport_1.
Assuming now that e1 and e2 are normalized base vectors, the electric field can be
expanded as
E = E1 e1 + E2 e2 ,
E
1.
E
2
In COMSOL, the Jones vector elements are complex numbers and the
Jones vectors are not normalized.
If the Jones vector elements have the same phase or a π phase difference, the Jones
vector represents a linear polarization state. A phase difference of ±π/2 between the
52 | CHAPTER 2: RF MODELING
two Jones vector elements defines a circular polarization state. For other phase
differences, the Jones vector represents elliptic polarization states.
For more details about the different types of periodic ports, see
• Periodic
• Diffraction Order
• Orthogonal Polarization
Polarization Plots
Default Polarization plots are automatically generated for Periodic ports in 3D and in
2D, if the Electric field components for setting in the Components section for the physics
interface is set to Three-component vector. The Polarization plot includes polarization
ellipses for each diffraction order. The polarization ellipses are generated by plotting
the in-plane Jones vector element versus the out-of-plane Jones vector element for a
phase change of 2π.
Thus, the variable emw.JROOP_0_n1 represents the Jones vector element for the
out-of-plane mode in 3D with mode numbers m = 0 and n = -1 for a port located on
the same boundary as the exciting Periodic port.
There are also variables for the norm of the Jones vector, named as tag.normJX_Z,
where X and Z represents the boundary location and mode number, respectively. Thus,
the variable emw.normJR_0_n1 represents the Jones vector norm in 3D with mode
numbers m = 0 and n = -1 for a port located on the same boundary as the exciting
Periodic port.
The base vectors are available for evaluation on the port boundaries, with variable
names like tag.eJXYYY[xyz]_Z, where tag, X, YYY, and Z have the same meanings as
for the variables discussed above and [xyz] is one of the Cartesian components. Thus,
emw.eJTIPx_p1, is the x-component of the in-plane base vector for a 2D port, with
mode number m = +1, that is not located on the same boundary as the exciting port.
54 | CHAPTER 2: RF MODELING
Lossy Eigenvalue Calculations
In mode analysis and eigenfrequency analysis, it is usually the primary goal to find a
propagation constant or an eigenfrequency. These quantities are often real-valued
although it is not necessary. If the analysis involves some lossy part, like a nonzero
conductivity or an open boundary, the eigenvalue is complex. In such situations, the
eigenvalue is interpreted as two parts (1) the propagation constant or eigenfrequency
and (2) the damping in space and time.
In this section:
• Eigenfrequency Analysis
• Mode Analysis and Boundary Mode Analysis
Eigenfrequency Analysis
The eigenfrequency analysis solves for the eigenfrequency of a model. The
time-harmonic representation of the fields is more general and includes a complex
parameter in the phase
˜ jωt ˜ –λ t
E ( r, t ) = Re ( E ( r T )e ) = Re ( E ( r )e )
where the eigenvalue, (−λ) = −δ + jω, has an imaginary part representing the
eigenfrequency, and a real part responsible for the damping. It is often more common
to use the quality factor or Q-factor, which is derived from the eigenfrequency and
damping
ω
Q fact = ---------
2δ
These situations may require special treatment, especially since it can lead to “singular
matrix” or “undefined value” messages if not treated correctly. Under normal
circumstances, the automatically generated solver settings should handle the cases
described in the table above. However, the following discussion provide some
background to the problem of defining the eigenvalue linearization point. The
complication is not only the nonlinearity itself, it is also the way it enters the equations.
For example the impedance boundary conditions with nonzero boundary conductivity
has the term
ε 0 μ 0 μ rbnd
– ( – λ ) ------------------------------------------ ( n × ( n × H ) )
σ bnd
ε rbnd + -----------------
( – λ )ε 0
where (−λ) = −δ + jω. When the solver starts to solve the eigenfrequency problem it
linearizes the entire formulation with respect to the eigenvalue around a certain
linearization point. By default this linearization point is set to the value provided to the
Search for eigenvalues around field, for the three cases listed in the table above.
Normally, this should be a good value for the linearization point. For instance, for the
56 | CHAPTER 2: RF MODELING
impedance boundary condition, this avoids setting the eigenvalue λ to zero in the
denominator in the equation above. For other cases than those listed in the table
above, the default linearization point is zero.
If the default values for the linearization point is not suitable for your particular
problem, you can manually provide a “good” linearization point for the eigenvalue
solver. Do this in the Eigenvalue node (not the Eigenfrequency node) under the Solver
Sequence node in the Study branch of the Model Builder. A solver sequence can be
generated first. In the Linearization Point section, select the Transform point check box
and enter a suitable value in the Point field. For example, if it is known that the
eigenfrequency is close to 1 GHz, enter the eigenvalue 1[GHz] in the field.
In many cases it is enough to specify a good linearization point and then solve the
problem once. If a more accurate eigenvalue is needed, an iterative scheme is necessary:
1 Specify that the eigenvalue solver only searches for one eigenvalue. Do this either
for an existing solver sequence in the Eigenvalue node or, before generating a solver
sequence, in the Eigenfrequency node.
2 Solve the problem with a “good” linearization point. As the eigenvalue shifts, use
the same value with the real part removed from the eigenvalue or, equivalently, use
the real part of the eigenfrequency.
3 Extract the eigenvalue from the solution and update the linearization point and the
shift.
4 Repeat until the eigenvalue does not change more than a desired tolerance.
˜ jωt – αz
E ( r, t ) = Re ( E ( r )e )
α = δ z + jβ = – λ ,
α = jβ = – λ .
In the table above, lambda is the eigenvalue and k0 is the vacuum wave number.
The suffix to the variables above indicates that the variables are defined for the port
named i.
58 | CHAPTER 2: RF MODELING
• For a list of the studies available by physics interface, see The RF
Module Physics Interface Guide
• Studies and Solvers in the COMSOL Multiphysics Reference Manual
In general electrical circuits connect to other physics interfaces via one or more of three
special circuit features:
• External I vs. U
• External U vs. I
• External I-Terminal
60 | CHAPTER 2: RF MODELING
These features either accept a voltage measurement from the connecting noncircuit
physics interface and return a current from a or the Electrical Circuit interface or the
other way around.
• A choice is made in the Settings window for the noncircuit physics interface feature,
which then announces (that is, includes) the coupling to the Electrical Circuit
interface. Its voltage or current is then included to make it visible to the connecting
circuit feature.
• A voltage or current that has been announced (that is, included) is selected in a
feature node’s Settings window.
• Apply the voltage or current from the connecting “External” circuit feature as an
excitation in the noncircuit physics interface.
1 In the Model Builder, right-click the Study node and select Show Default Solver.
2 Expand the Solver>Dependent Variables node and click the state node, in this
example, Current through device R1 (comp1.currents). The variable name is shown in
the Settings window for State.
62 | CHAPTER 2: RF MODELING
Solving
Some modeling errors lead to the error message The DAE is structurally
inconsistent being displayed when solving. This error typically occurs
from having an open current loop, from connecting voltage sources in
parallel, or connecting current sources in series.
In this respect, the predefined coupling features are also treated as (ideal)
voltage or current sources. The remedy is to close current loops and to
connect resistors in series with voltage sources or in parallel with current
sources.
Postprocessing
The Electrical Circuits interface, unlike most of the other physics interfaces, solves for
a relatively large number of global dependent variables (such as voltages and currents),
instead of solving for a few space-varying fields (such as temperature or displacement).
For this reason, the Electrical Circuit interface does not provide default plots when
computing a study.
The physics interface defines a number of variables that can be used in postprocessing.
All variables defined by the Electrical Circuit interface are of a global scope, and can
be evaluated in a Global Evaluation node (under Derived Values). In addition, the time
evolution or dependency on a parameter can be plotted in a Global plot (under a 1D
Plot Group node).
The physics interface defines a Node voltage variable for each electrical node in the
circuit, with name cir.v_name, where cir is the physics interface Label and <name>
is the node Name. For each two-pin component, the physics interface also defines
variables containing the voltage across it and the current flowing through it.
R Resistor
C Capacitor
L Inductor
V Voltage Source
I Current Source
E Voltage-Controlled Voltage Source
F Current-Controlled Current Source
G Voltage-Controlled Current Source
H Current-Controlled Voltage Source
D Diode
Q NPN BJT and PNP BJT
M n-Channel MOSFET and p-Channel MOSFET
X Subcircuit Instance
According to SPICE specification, the first line in the netlist file is assumed to be the
title of the netlist and it is ignored by the parser.
64 | CHAPTER 2: RF MODELING
SPICE Export
The SPICE Export functionality creates a SPICE netlist file containing a description of
the circuit represented by the physics interface. This functionality can be accessed from
the physics interface context menu (right-click the physics interface node and select
Export SPICE Netlist). After specifying a filename, the circuit is exported and messages
from the export process display in the Messages window. During the export process, a
series of operations are performed:
The title of the exported netlist file is the model’s filename, and the time, date, and
version of COMSOL Multiphysics is added as a comment in the netlist file.
Reference
1. http://bwrc.eecs.berkeley.edu/Classes/IcBook/SPICE/
In this section:
66 | CHAPTER 2: RF MODELING
CHOOSING APPROPRIATE ASYMPTOTIC WAVEFORM EVALUATION (AWE)
EXPRESSIONS
The simulation time may vary depending on the user input for the AWE expressions.
Any model variable works as an AWE expression, so long as it has a smooth curve shape
like a Gaussian pulse as a function of frequency. The absolute value of S21,
abs(comp1.emw.S21), often works as the input for the AWE expression in the case of
two-port bandpass filters. For one-port devices like antennas, the absolute value of S11
is a good choice. If the frequency response of the AWE expression contains an infinite
gradient — the case for the S11 value of an antenna with excellent impedance matching
at a single frequency point — the simulation will take longer to complete because it
requires many data points to describe the sharp dip. When the loss in a one-port device
is negligible, an alternative expression such as sqrt(1-abs(comp1.emw.S11)^2) may
work more efficiently than using abs(comp1.emw.S11) directly.
DATA MANAGEMENT
With a very fine frequency step simulation, the solutions contain a lot of data. As a
result, the model file size will increase tremendously when it is saved. By selecting the
Store fields in output check box in the Values of Dependent Variables section of the
Frequency Domain study step settings, it is possible to define for what part of the model
the computed solution should be saved. When only S-parameters are of interest, it is
not necessary to store all of the field solutions. Instead, only store the field on the
selections for the port and lumped port boundaries, as those will be used for the
S-parameter calculations.
In the Values of Dependent Variables section, change the selection in the Store fields in
output combo box from All to For selections and then add the explicit selections that
include the port and lumped port boundaries. The explicit selection can be easily
created from the port and lumped port feature by clicking Create Selection icon in the
Boundary Selection settings once the selection is specified.
DATA MANAGEMENT
The Store fields in output check box in the Values of Dependent Variables section can be
applied to the Frequency Domain, Modal study — if you are interested only in
68 | CHAPTER 2: RF MODELING
S-parameters. By storing solutions only on port or lumped port boundaries, the saved
model file size will decrease a lot.
The Part Library contains complex shapes, frequently required for RF simulations,
including the following parts:
• Connectors
• Surface mount devices
• Waveguides
The parts are built from partially parameterized sequences of geometry instructions.
For example, you can load the connector_sma_flange2 part into a model and then
specify the outer and inner radii of the coaxial structure. RF parts usually include
selections for the conductive boundary that make it easy to apply PEC boundary
conditions while setting up the physics.
The following RF tutorials use the Part Library to create their geometry
sequences:
70 | CHAPTER 2: RF MODELING
Material Libraries
The RF Module features a Material Library with material properties for substrate
materials to assist in modeling RF, microwave, and millimeter-wave circuit boards.
The RF Material Library contains the material property data from following
companies’ products:
ROGERS CORPORATION
RO4000® Laminates, RT/duroid® Laminates, RO3000® Laminates, XT/duroid®
Copper Clad Laminates, TMM® Laminates, TC™ Series Laminates, Kappa™ 438
Laminates, CuClad® & IsoClad® Series Laminate, DiClad® Series Laminates, CLTE
Series™ Materials, AD Series™ Laminates, and Bondply/Prepreg
ISOLA GROUP
185HR, 254, 370 TURBO®, 370HR, A11, Astra® MT77, DE104, ED130UV,
FR402, FR406, FR406N, FR408, FR408HR, G200, I-Speed®, I-Tera® MT40,
I-Tera® MT40 (RF/MW), IS300MD, IS400, IS400HR, IS410, IS415, IS420,
IS680, IS680 AG, IS680 AG -348, P25N, P95/P25, P96/P26, Tachyon® 100G,
TerraGreen®, and TerraGreen® (RF/MW)
The foregoing materials are the property of Isola Group. For product information
concerning such materials, see https://www.isola-group.com/products/. For other
questions concerning such materials, please contact https://www.isola-group.com/
contact-us/.
PREMIX GROUP
PREPERM® Standard Grades, PREPERM® Radome Grades, PREPERM®
High-temperature Grades, PREPERM® Flexible Grades, and PREPERM® H-series
Grades
The foregoing materials are the property of Premix Oy. For product information
concerning such materials, see https://www.preperm.com/products/raw-materials/
. For other questions concerning such materials, please contact https://
www.preperm.com/contact-us/.
MATERIAL LIBRARIES | 71
72 | CHAPTER 2: RF MODELING
3
Electromagnetics Theory
In this chapter:
• Maxwell’s Equations
• Special Calculations
• Electromagnetic Quantities
See also:
73
Maxwell’s Equations
In this section:
The equations can be formulated in differential or integral form. The differential form
are presented here, because it leads to differential equations that the finite element
method can handle. For general time-varying fields, Maxwell’s equations can be
written as
The first two equations are also referred to as Maxwell–Ampère’s law and Faraday’s
law, respectively. Equation three and four are two forms of Gauss’ law, the electric and
magnetic form, respectively.
∂ρ
∇⋅J = – -----
-
∂t
Out of the five equations mentioned, only three are independent. The first two
combined with either the electric form of Gauss’ law or the equation of continuity
form such an independent system.
Constitutive Relations
To obtain a closed system, the constitutive relations describing the macroscopic
properties of the medium, are included. They are given as
D = ε0 E + P
B = μ0 ( H + M )
J = σE
Here ε0 is the permittivity of vacuum, μ0 is the permeability of vacuum, and σ the
electrical conductivity. In the SI system, the permeability of a vacuum is chosen to be
4π·10−7 H/m. The velocity of an electromagnetic wave in a vacuum is given as c0 and
the permittivity of a vacuum is derived from the relation
1 – 12 1 –9
ε 0 = ----------
2
= 8.854 ⋅ 10 F/m ≈ --------- ⋅ 10 F/m
c0 μ0 36π
The electric polarization vector P describes how the material is polarized when an
electric field E is present. It can be interpreted as the volume density of electric dipole
moments. P is generally a function of E. Some materials can have a nonzero P also
when there is no electric field present.
MAXWELL’S EQUATIONS | 75
The magnetization vector M similarly describes how the material is magnetized when
a magnetic field H is present. It can be interpreted as the volume density of magnetic
dipole moments. M is generally a function of H. Permanent magnets, however, have a
nonzero M also when there is no magnetic field present.
For linear materials, the polarization is directly proportional to the electric field,
P = ε0χeE, where χe is the electric susceptibility. Similarly in linear materials, the
magnetization is directly proportional to the magnetic field, M = χmH, where χm is the
magnetic susceptibility. For such materials, the constitutive relations can be written
D = ε 0 ( 1 + χ e )E = ε 0 ε r E = εE
B = μ 0 ( 1 + χ m )H = μ 0 μ r H = μH
The parameter εr is the relative permittivity and μr is the relative permeability of the
material. These are usually scalar properties but they can, for a general anisotropic
material, be 3-by-3 tensors. The properties ε and μ (without subscripts) are the
permittivity and permeability of the material.
D = ε0 εr E + Dr
The field Dr is the remanent displacement, which is the displacement when no electric
field is present.
Similarly, a generalized form of the constitutive relation for the magnetic field is
B = μ0 μr H + Br
where Br is the remanent magnetic flux density, which is the magnetic flux density
when no magnetic field is present.
e
J = σE + J
B = ∇×A
∂A
E = – ∇V – -------
∂t
The defining equation for the magnetic vector potential is a direct consequence of the
magnetic Gauss’ law. The electric potential results from Faraday’s law.
Electromagnetic Energy
The electric and magnetic energies are defined as
D T
∂D
We = V 0 E ⋅ dD dV = V 0 E ⋅ -------
∂t
dt dV
B T
∂B
Wm = H ⋅ dB dV = - dt dV
0 H ⋅ ------
V 0
V ∂t
The time derivatives of these expressions are the electric and magnetic power
∂D
Pe = V E ⋅ -------
∂t
dV
∂B
Pm = V H ⋅ ------
∂t
- dV
These quantities are related to the resistive and radiative energy, or energy loss,
through Poynting’s theorem (Ref. 3)
∂D ∂B
– V E ⋅ -------
∂t
+ H ⋅ ------- dV = J ⋅ E dV + ( E × H ) ⋅ n dS
∂t V °S
where V is the computation domain and S is the closed boundary of V.
The first term on the right-hand side represents the resistive losses,
Ph = V J ⋅ E dV
MAXWELL’S EQUATIONS | 77
which result in heat dissipation in the material. (The current density J in this
expression is the one appearing in Maxwell–Ampère’s law.)
The second term on the right-hand side of Poynting’s theorem represents the radiative
losses,
Pr =
°S ( E × H ) ⋅ n dS
The quantity S = E × H is called the Poynting vector.
Under the assumption the material is linear and isotropic, it holds that
∂D ∂E
E ⋅ ------- = εE ⋅ ------- = ∂ --- εE ⋅ E
1
∂t ∂t ∂t 2
∂B ∂B
H ⋅ ------- = --- B ⋅ ------- = ∂ ------- B ⋅ B
1 1
∂t μ ∂t ∂ t 2μ
By interchanging the order of differentiation and integration (justified by the fact that
the volume is constant and the assumption that the fields are continuous in time), this
equation results:
∂
V --2- εE ⋅ E + ------
- B ⋅ B dV =
1 1
–
∂t 2μ V J ⋅ E dV + °S ( E × H ) ⋅ n dS
The integrand of the left-hand side is the total electromagnetic energy density
1 1
w = w e + w m = --- εE ⋅ E + ------- B ⋅ B
2 2μ
Material Properties
Until now, there has only been a formal introduction of the constitutive relations.
These seemingly simple relations can be quite complicated at times. There are four
main groups of materials where they require some considerations. A given material can
belong to one or more of these groups. The groups are:
• Inhomogeneous materials
• Anisotropic materials
• Nonlinear materials
• Dispersive materials
For anisotropic materials, the field relations at any point are different for different
directions of propagation. This means that a 3-by-3 tensor is required to properly
define the constitutive relations. If this tensor is symmetric, the material is often
referred to as reciprocal. In these cases, the coordinate system can be rotated in such
a way that a diagonal matrix is obtained. If two of the diagonal entries are equal, the
material is uniaxially anisotropic. If none of the elements have the same value, the
material is biaxially anisotropic (Ref. 2). An example where anisotropic parameters
are used is for the permittivity in crystals (Ref. 2).
Finally, dispersion describes changes in the velocity of the wave with wavelength. In
the frequency domain, dispersion is expressed by a frequency dependence in the
constitutive laws.
The physics-specific domain material properties are by default taken from the material
specification. The material properties are inputs to material laws or constitutive
relations that are defined on the feature level below the physics interface node in the
model tree. There is one editable default domain feature (wave equation) that initially
represents a linear isotropic material. Domains with different material laws are specified
by adding additional features. Some of the domain parameters can either be a scalar or
a matrix (tensor) depending on whether the material is isotropic or anisotropic.
MAXWELL’S EQUATIONS | 79
In a similar way, boundary, edge, and point settings are specified by adding the
corresponding features. A certain feature might require one or several fields to be
specified, while others generate the conditions without user-specified fields.
n2 × ( E1 – E2 ) = 0
n2 ⋅ ( D1 – D2 ) = ρs
n2 × ( H1 – H2 ) = Js
n2 ⋅ ( B1 – B2 ) = 0
where ρs and Js denote surface charge density and surface current density,
respectively, and n2 is the outward normal from medium 2. Of these four conditions,
only two are independent. One of the first and the fourth equations, together with one
of the second and third equations, form a set of two independent conditions.
A consequence of the above is the interface condition for the current density,
∂ρ s
n 2 ⋅ ( J 1 – J 2 ) = – --------
∂t
–n2 × E2 = 0
–n2 × H2 = Js
–n2 ⋅ D2 = ρs
–n2 ⋅ B2 = 0
Phasors
Whenever a problem is time-harmonic the fields can be written in the form
ˆ
E ( r, t ) = E ( r ) cos ( ωt + φ )
Instead of using a cosine function for the time dependence, it is more convenient to
use an exponential function, by writing the field as
ˆ ˆ jφ jωt ˜ jωt
E ( r, t ) = E ( r ) cos ( ωt + φ ) = Re ( E ( r )e e ) = Re ( E ( r )e )
˜
The field E ( r ) is a phasor (phase vector), which contains amplitude and phase
information of the field but is independent of t. One thing that makes the use of
phasors suitable is that a time derivative corresponds to a multiplication by jω,
∂------
E- ˜ jωt
= Re ( jωE ( r )e )
∂t
This means that an equation for the phasor can be derived from a time-dependent
equation by replacing the time derivatives by a factor jω. All time-harmonic equations
in this module are expressed as equations for the phasors. (The tilde is dropped from
the variable denoting the phasor.).
MAXWELL’S EQUATIONS | 81
For example, all plot functions visualize
˜
Re ( E ( r ) )
• S-Parameter Calculations
• Far-Field Calculations Theory
• References
S-Parameter Calculations
For high-frequency problems, voltage is not a well-defined entity, and it is necessary
to define the scattering parameters (S-parameter) in terms of the electric field. To
convert an electric field pattern on a port to a scalar complex number corresponding
to the voltage in transmission line theory an eigenmode expansion of the
electromagnetic fields on the ports needs to be performed. Assume that an eigenmode
analysis has been performed on the ports 1, 2, 3, … and that the electric field patterns
E1, E2, E3, … of the fundamental modes on these ports are known. Further, assume
that the fields are normalized with respect to the integral of the power flow across each
port cross section, respectively. This normalization is frequency dependent unless
TEM modes are being dealt with. The port excitation is applied using the fundamental
eigenmode, the mode with subscript 1. The computed electric field Ec on the port
consists of the excitation plus the reflected field. That is, on the port boundary where
there is an incident wave, the computed field can be expanded in terms of the mode
fields as
Ec = E1 + Si1 Ei ,
i=1
Ec = Si1 Ei
i=1
SPECIAL CALCULATIONS | 83
The S-parameter for the mode with index k is then given by multiplying with the
conjugate of the mode field for mode k and integrating over the port boundary. Since
the mode fields for the different modes are orthogonal, the following relations are
obtained for the S-parameters
*
( ( 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
*
( E c ⋅ E 3 ) dA 3
port 3
S 31 = ----------------------------------------------
-
*
( E 3 ⋅ E 3 ) dA 3
port 3
and so on. To get S22 and S12, excite port number 2 in the same way.
1 *
S av = --- Re ( E × H )
2
The amount of power flowing out of a port is given by the normal component of the
Poynting vector,
1 *
n ⋅ S av = n ⋅ --- Re ( E × H )
2
Below the cutoff frequency the power flow is zero, which implies that it is not possible
to normalize the field with respect to the power flow below the cutoff frequency. But
in this region the S-parameters are trivial and do not need to be calculated.
In the following subsections the power flow is expressed directly in terms of the electric
field for TE, TM, and TEM waves.
E = – Z TE ( n × H )
ωμ
Z TE = -------
β
ω is the angular frequency of the wave, μ the permeability, and β the propagation
constant. The power flow then becomes
1 * 1 * 1 2
n ⋅ S av = --- n ⋅ Re ( E × H ) = – --- Re ( E ⋅ ( n × H ) ) = -------------- E
2 2 2Z TE
TM Waves
For TM waves it holds that
1
H = ----------- ( n × E )
Z TM
β
Z TM = -------
ωε
1 * 1 *
n ⋅ S av = --- n ⋅ Re ( E × H ) = --------------- ( n ⋅ Re ( E × ( n × E ) ) )
2 2Z TM
1 2
= --------------- n × E
2Z TM
TEM Waves
For TEM waves it holds that
1
H = --------------- ( n × E )
Z TEM
μ
Z TEM = ---
ε
SPECIAL CALCULATIONS | 85
The power flow then becomes
1 * 1 2 1 2
n ⋅ S av = --- n ⋅ Re ( E × H ) = ------------------ n × E = ------------------ E
2 2Z TEM 2Z TEM
where the last equality holds because the electric field is tangential to the port.
jk
E p = ------ r 0 × [ n × E – ηr 0 × ( n × H ) ] exp ( jkr ⋅ r 0 ) dS
4π
S
jk
Ep =
4π
λ ------ r 0 × [ n × E – ηr 0 × ( n × H ) ] exp ( jkr ⋅ r 0 ) dS
S
In both cases the integration is performed on a closed boundary. In the scattered field
formulation, where the total electric field is the sum of the background field and the
scattered field, the far field only gets contributions from the scattered field, since the
contributions from the background field cancel out when integrated over all parts of
the closed boundary.
For scattering problems, the far field in COMSOL Multiphysics is identical to what in
physics is known as the “scattering amplitude”.
The antenna 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 antenna.
• 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 impedance:
η = μ⁄ε
Thus the unit vector r0 can be interpreted as the direction defined by the angular
position ( θ, ϕ ) and Ep is the far field in this direction.
Because the far field is calculated in free space, the magnetic field at the far-field point
is given by
r0 × Ep
H p = -------------------
η0
The Poynting vector gives the power flow of the far field:
* 2
r 0 ⋅ S = r 0 ⋅ Re ( E p × H p ) ∼ E p
References
1. D.K. Cheng, Field and Wave Electromagnetics, 2nd ed., Addison-Wesley, 1991.
4. R.K. Wangsness, Electromagnetic Fields, 2nd ed., John Wiley & Sons, 1986.
SPECIAL CALCULATIONS | 87
Electromagnetic Quantities
Table 3-1 shows the symbol and SI unit for most of the physical quantities that are
included with this module.
TABLE 3-1: ELECTROMAGNETIC QUANTITIES
ELECTROMAGNETIC QUANTITIES | 89
90 | CHAPTER 3: ELECTROMAGNETICS THEORY
4
This chapter discusses the physics interfaces found under the Radio Frequency
branch ( ).
In this chapter:
91
The Electromagnetic Waves,
Frequency Domain Interface
The Electromagnetic Waves, Frequency Domain (emw) interface ( ), found under the
Radio Frequency branch ( ) when adding a physics interface, is used to solve for
time-harmonic electromagnetic field distributions.
For this physics interface, the maximum mesh element size should be limited to a
fraction of the wavelength. The domain size that can be simulated thus scales with the
amount of available computer memory and the wavelength. The physics interface
supports the Frequency Domain, Eigenfrequency, Mode Analysis, and Boundary
Mode Analysis study types. The Frequency Domain study type is used for source driven
simulations for a single frequency or a sequence of frequencies. The Eigenfrequency
study type is used to find resonance frequencies and their associated eigenmodes in
resonant cavities.
This physics interface solves the time-harmonic wave equation for the electric field.
When this physics interface is added, these default nodes are also added to the Model
Builder — Wave Equation, Electric, Perfect Electric Conductor, and Initial Values. Then,
from the Physics toolbar, add other nodes that implement, for example, boundary
conditions. You can also right-click Electromagnetic Waves, Frequency Domain to select
physics features from the context menu.
The Mode analysis study type is applicable only for 2D and 2D axisymmetric cross
sections of waveguides and transmission lines where it is used to find allowed
propagating modes. Boundary mode analysis is used for the same purpose in 2D, 2D
axisymmetry, and 3D and applies to boundaries representing waveguide ports.
Physics-Controlled Mesh
The physics-controlled mesh is controlled from the Mesh node’s Settings window (if the
Sequence type is Physics-controlled mesh). In the table in the Physics-Controlled Mesh
section, find the physics interface in the Contributor column and select or clear the
check box in the Use column on the same row for enabling (the default) or disabling
contributions from the physics interface to the physics-controlled mesh.
When the Use check box for the physics interface is selected, this invokes a parameter
for the maximum mesh element size in free space. The physics-controlled mesh
automatically scales the maximum mesh element size as the wavelength changes in
When the Use check box is selected for the physics interface, in the section for the
physics interface below the table, choose one of the four options for the Maximum mesh
element size control parameter — From study (the default), User defined, Frequency, or
Wavelength. When From study is selected, 1/5 of the vacuum wavelength from the
highest frequency defined in study step is used for the maximum mesh element size.
For the option User defined, enter a suitable Maximum element size in free space. For
example, 1/5 of the vacuum wavelength or smaller. When Frequency is selected, enter
the highest frequency intended to be used during the simulation. The maximum mesh
element size in free space is 1/5 of the vacuum wavelength for the entered frequency.
For the Wavelength option, enter the smallest vacuum wavelength intended to be used
during the simulation. The maximum mesh element size in free space is 1/5 of the
entered wavelength.
When Resolve wave in lossy media is selected, the outer boundaries of lossy media
domains are meshed with a maximum mesh element size in free space given by the
minimum value of half a skin depth and 1/5 of the vacuum wavelength.
The maximum mesh element size in dielectric media is that in free space divided by the
square root of the product of the relative permittivity and permeability.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is emw.
For Scattered field select a Background wave type according to the following table:
TABLE 4-1: BACKGROUND WAVE TYPE BASED ON COMPONENT DIMENSION
User Defined
Enter the component expressions for the Background electric field Eb (SI unit: V/m).
The entered expressions must be differentiable.
Gaussian Beam
For Gaussian beam select the Gaussian beam type — Paraxial approximation (the default)
or Plane wave expansion.
If the beam spot radius is smaller than the wavelength, evanescent plane waves need to
be included in the expansion. The evanescent waves decay exponentially in the
propagation direction, why it only makes sense to model such tightly focused beams if
the focal plane coincides with the input boundary. If the focal plane is located inside
For Plane wave expansion select Wave vector distribution type — Automatic (the default)
or User defined. For Automatic also check Allow evanescent waves, to include evanescent
waves in the plane wave expansion. For User defined also enter values for the Wave
vector count Nk (the default value is 13) and Maximum transverse wave number kt,max
(SI unit: rad/m, default value is (2*(sqrt(2*log(10))))/emw.w0). Use an odd
number for the Wave vector count Nk to make sure that a wave vector pointing in the
main propagation direction is included in the plane-wave expansion. The Wave vector
count Nk specifies the number of wave vectors that will be included per transverse
dimension. So for 3D the total number of wave vectors will be Nk·Nk.
A plane wave expansion with a finite number of plane waves included will make the
field periodic in the plane orthogonal to the main propagation direction. If the
separation between the transverse wave vector components, given by 2kt,max/
(Nk − 1), is too small, replicas of the Gaussian beam background field can appear. To
avoid that, increase the value for the Wave vector count Nk.
The number of plane waves included in the expansion can be quite large, especially for
3D. For instance, using the default settings, 2·13·13 = 338 plane waves will be
included (the factor 2 accounts for the two possible polarizations for each wave
vector). Thus, initializing the plane-wave expansion for the Gaussian beam
background field can take some time in 3D.
For more information about the Gaussian beam theory, see Gaussian Beams as
Background Fields and Input Fields.
• Select a Beam orientation: Along the x-axis (the default), Along the y-axis, or for 3D
components, Along the z-axis.
• Enter a Beam radius w0 (SI unit: m). The default is 20π/emw.k0 m (10 vacuum
wavelengths).
• Enter a Focal plane along the axis p0 (SI unit: m). The default is 0 m.
• Enter an Electric field amplitude E0 (SI unit: V/m). The default is 1 V/m.
• Enter a Roll angle (SI unit: rad), which is a right-handed rotation with respect to the
+x direction. The default is 0 rad, corresponding to polarization along the
+z direction.
• Enter a Pitch angle (SI unit: rad), which is a right-handed rotation with respect to
the +y direction. The default is 0 rad, corresponding to the initial direction of
propagation pointing in the +x direction.
Figure 4-1: Schematic of the directions for the wave vector k, the electric field E0, and the
roll, pitch and yaw rotations. The top image represents an initial wave propagating in the
x direction with a polarization along the z direction.
– jmϕ
E b ( r, ϕ, z ) = Ẽ b ( r, z )e ,
where
˜ ( r, z ) = E ( r̂ – jmϕ̂ )e –jkz ,
E b 0
m is the azimuthal mode number (+1 or -1) varying depending on the Circular
polarization type and Direction of propagation settings, and r̂ and ϕ̂ are the unit vectors
in the r and ϕ directions, respectively.
Select the Electric field components solved for — Three-component vector, Out-of-plane
vector, or In-plane vector. Select:
For 2D components, assign a wave vector component to the Out-of-plane wave number
field. For 2D axisymmetric components, assign an integer constant or an integer
parameter expression to the Azimuthal mode number field.
ANALYSIS METHODOLOGY
From the Methodology options list, select one of three solver configurations: Robust,
Intermediate, or Fast (the default).
The settings of each methodology option are found in Solver Configurations and the
subsidiary nodes.
For Use manual port sweep enter a Sweep parameter name to assign a specific name to
the parameter that controls the port number solved for during the sweep. Before
making the port sweep, the parameter must also have been added to the list of
parameters in the Parameters section of the Parameters node under the Global
Definitions node. This process can be automated by clicking the Configure Sweep
Settings button. The Configure Sweep Settings button helps add a necessary port sweep
Select Export Touchstone file and the S-parameters are subject to Touchstone file export.
Click Browse to locate the file, or enter a file name and path. Select an Parameter format
(value pair): Magnitude angle, Magnitude (dB) angle, or Real imaginary.
Enter a Reference impedance for Touchstone file export Zref (SI unit: Ω) that is used only
for the header in the exported Touchstone file. The default is 50 Ω.
DEPENDENT VARIABLES
The dependent variables (field variables) are for the Electric field E and its components
(in the Electric field components fields). The name can be changed but the names of
fields and dependent variables must be unique within a model.
DISCRETIZATION
Select the shape order for the Electric field dependent variable — Linear, Linear type 2,
Quadratic (the default), Quadratic type 2, Cubic, or Cubic type 2. For more information
about the Discretization section, see Settings for the Discretization Sections in the
COMSOL Multiphysics Reference Manual.
DOMAIN
BOUNDARY CONDITIONS
With no surface currents present the boundary conditions
n2 × ( E1 – E2 ) = 0
n2 × ( H1 – H2 ) = 0
need to be fulfilled. Because E is being solved for, the tangential component of the
electric field is always continuous, and thus the first condition is automatically fulfilled.
The second condition is equivalent to the natural boundary condition
and is therefore also fulfilled. These conditions are available (listed in alphabetical
order):
2
∇ × ( μ r– 1 ∇ × E ) – k 0 ε rc E = 0
for the time-harmonic and eigenfrequency problems. The wave number of free space
k0 is defined as
ω
k 0 = ω ε 0 μ 0 = -----
c0
˜
E ( x, y, z ) = E ( x, y ) exp ( – ik z z ) .
˜ 2 ˜
( ∇ – ik z z ) × [ μ r– 1 ( ∇ – ik z z ) × E ] – k 0 ε rc E = 0 ,
Notice that the ansatz above just explains how the wave equation is modified when the
out-of-plane wave vector component kz is not zero. As an example, for a plane wave
with a nonzero out-of-plane wave vector component, the electric field is of course
given by
˜
E ( x, y, z ) = E ( x, y ) exp ( – ik z z ) = A exp ( – i ( k x x + k y y + k z z ) ) ,
where A is a constant amplitude and kx, ky, and kz are the wave vector components.
In 2D axisymmetry, the electric field varies with the azimuthal mode number m as
∇ – im ˜ ˜
----- ϕ × μ r–1 ∇ – i m
----- ϕ × E – k 0 ε rc E = 0 ,
2
r r
ω
Q fact = ---------
2δ
Using the relation εr = n2, where n is the refractive index, the equation can
alternatively be written
2 2
∇ × ( ∇ × E ) – k0 n E = 0
When the equation is written using the refractive index, the assumption is that μr = 1
and σ = 0 and only the constitutive relations for linear materials are available. When
solving for the scattered field the same equations are used but E = Esc + Ei and Esc is
the dependent variable.
The Divergence Constraint subnode is available from the context menu (right-click the
parent node) or from the Physics toolbar, Attributes menu.
• Relative permittivity
• Refractive index
• Loss tangent, loss angle
• Loss tangent, dissipation factor
• Dielectric loss
As an example, the material can be specified to use a Refractive index
material model. Then the real and imaginary parts of the refractive index
can be synchronized to compute the complex relative permittivity, if the
Electric displacement field model is set to Relative permittivity.
When synchronizing to the Loss tangent, loss angle and Loss tangent,
dissipation factor Electric displacement field models, the loss angle δ and the
dissipation factor tanδ, respectively, must be converted to isotropic values.
Relative Permittivity
When Relative permittivity is selected, the default Relative permittivity εr takes values
From material. Select Effective medium to add an Effective Medium subnode, or for User
defined select Isotropic, Diagonal, Symmetric, or Full and enter values or expressions in
the field or matrix.
Refractive Index
When Refractive index is selected, the default Refractive index n and Refractive index,
imaginary part k take the values From material. To specify the real and imaginary parts
of the refractive index and assume a relative permeability of unity and zero
Dielectric Loss
When Dielectric loss is selected, the default Relative permittivity ε′ and Relative
permittivity (imaginary part) ε″ take values From material. For User defined select
M 2
fj ωP
εr ( ω ) = ε∞ + ---------------------------------------
2
ω 0j – ω + iΓ j ω
2
-
j=1
When Drude-Lorentz dispersion model is selected, the default Relative permittivity, high
frequency ε∞ (dimensionless) takes its value From material. For User defined select
Isotropic, Diagonal, Symmetric, or Full and enter a value or expression in the field or
matrix.
In the table, enter values or expressions in the columns for the Oscillator strength,
Resonance frequency (rad/s), and Damping in time (rad/s).
Δε k
ε ( ω ) = ε∞ + ---------------------
1 + iωτ k
-
k
When Debye dispersion model is selected, the default Relative permittivity, high
frequency ε∞ (dimensionless) takes its value From material. For User defined select
Isotropic, Diagonal, Symmetric, or Full and enter a value or expression in the field or
matrix.
MAGNETIC FIELD
Select the Constitutive relation — Relative permeability (the default) or Magnetic losses.
• For Relative permeability the relative permeability μr uses values From material. For
User defined select Isotropic, Diagonal, Symmetric, or Full based on the characteristics
of the magnetic field, and then enter values or expressions in the field or matrix. If
Effective medium is selected, the Effective Medium subnode is available from the
context menu (right-click the parent node) or from the Physics toolbar, Attributes
menu.
• For Magnetic losses the default values for Relative permeability (real part) μ′ and
Relative permeability (imaginary part) μ″ are taken From material. For User defined
enter different values.
CONDUCTION CURRENT
By default, the Electrical conductivity σ (SI unit: S/m) uses values From material.
• For User defined select Isotropic, Diagonal, Symmetric, or Full based on the
characteristics of the current and enter values or expressions in the field or matrix.
• For Linearized resistivity the default values for the Reference temperature Tref (SI
unit: K), Resistivity temperature coefficient α (SI unit: 1/K), and Reference resistivity
ρ0 (SI unit: Ω⋅m) are taken From material. For User defined enter other values or
expressions for any of these variables.
• When Effective medium is selected, the Effective Medium subnode is available from
the context menu (right-click the parent node) or from the Physics toolbar,
Attributes menu.
• When Archie’s law is selected, the Archie’s Law subnode is available from the context
menu (right-click the parent node) or from the Physics toolbar, Attributes menu.
∇⋅J = 0
∇ ⋅ ( σA ) = 0
DIVERGENCE CONSTRAINT
Enter a value or expression for the Divergence condition variable scaling ψ0.
For the Electromagnetic Waves, Transient interface (and the Microwave Plasma
interface available with the Plasma Module) the SI unit is A/m and the default is 1 A/
m.
Initial Values
The Initial Values node adds an initial value for the electric field that can serve as an
initial guess for a nonlinear solver. Add additional Initial Values nodes from the Physics
toolbar.
INITIAL VALUES
Enter values or expressions for the initial values of the components of the Electric field
E (SI unit: V/m). The default values are 0 V/m.
J = σE + J e
Far-Field Domain
To set up a far-field calculation, add a Far-Field Domain node and specify the far-field
domains in its Settings window. Use Far-Field Calculation subnodes (one is added by
default) to specify all other settings needed to define the far-field calculation. If a
Perfectly Matched Layer (PML) node has been added before adding the Far-Field
Domain, all of the domains in the Electromagnetic Waves, Frequency Domain interface
adjacent to the PML are automatically selected by default. If there is no PML, all of
the domains are selected. The selection can be modified. In that case, select only a
homogeneous domain or domain group that is outside of all radiating and scattering
objects and which has the material settings of the far-field medium.
Far-Field Calculation
A Far-Field Calculation subnode is added by default to the Far-Field Domain node and
is used to select boundaries corresponding to a single closed surface surrounding all
radiating and scattering objects. By default, all exterior boundaries of the Far-Field
Domain are selected. If a Perfectly Matched Layer (PML) node has been added before
adding the Far-Field Domain, all exterior boundaries of the Far-Field Domain adjacent
to the PML are selected. Symmetry reduction of the geometry makes it relevant to
select boundaries defining a nonclosed surface. Also use this feature to indicate
symmetry planes and symmetry cuts applied to the geometry, and whether the selected
FAR-FIELD CALCULATION
Enter a Far-field variable name. The default is Efar.
Select as needed the Symmetry in the x=0 plane, Symmetry in the y=0 plane, or
Symmetry in the z=0 plane check boxes to use it your model when calculating the
far-field variable. The symmetry planes have to coincide with one of the Cartesian
coordinate planes.
When a check box is selected, also choose the type of symmetry to use from the
Symmetry type list that appears — Symmetry in E (PMC) or Symmetry in H (PEC). The
selection should match the boundary condition used for the symmetry boundary.
Using these settings, include the parts of the geometry that are not in the model for
symmetry reasons in the far-field analysis.
From the Boundary relative to domain list, select Inside or Outside (the default) to define
if the selected boundaries are defining the inside or outside of the far-field domain (that
is, whether facing away from infinity or toward infinity).
If perfectly matched layers are added to the model after the Far-Field
Domain is configured, then it is necessary to press the Reset Far-Field
Boundaries button to reassign all exterior boundaries.
2
σE
SAR = -------------- .
2ρ
n×E = 0
is a special case of the electric field boundary condition that sets the tangential
component of the electric field to zero. It is used for modeling of a lossless metallic
surface (for example, a ground plane) or as a symmetry type boundary condition. It
imposes symmetry for magnetic fields and “magnetic currents” and antisymmetry for
electric fields and electric currents. It supports induced electric surface currents and
thus any prescribed or induced electric currents (volume, surface, or edge currents)
flowing into a perfect electric conductor boundary is automatically balanced by
induced surface currents.
Js
J
I'
I
Js
Figure 4-2: The perfect electric conductor boundary condition is used on exterior and
interior boundaries representing the surface of a lossless metallic conductor or (on exterior
boundaries) representing a symmetry cut. The shaded (metallic) region is not part of the
model but still carries effective mirror images of the sources. Note also that any current
flowing into the boundary is perfectly balanced by induced surface currents. The
tangential electric field vanishes at the boundary.
CONSTRAINT SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. For information about the
n×H = 0
is a special case of the surface current boundary condition that 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. Electric currents (volume, surface, or edge
currents) are not allowed to flow into a perfect magnetic conductor boundary as that
would violate current conservation. On interior boundaries, the perfect magnetic
conductor boundary condition literally sets the tangential magnetic field to zero,
which in addition to setting the surface current density to zero also makes the
tangential electric field discontinuous.
Js=0
I'
I
J=0
Figure 4-3: The perfect magnetic conductor boundary condition is used on exterior
boundaries representing the surface of a high impedance region or a symmetry cut. The
shaded (high impedance) region is not part of the model but nevertheless carries effective
mirror images of the sources. Note also that any electric current flowing into the boundary
is forbidden as it cannot be balanced by induced electric surface currents. The tangential
magnetic field vanishes at the boundary. On interior boundaries, the perfect magnetic
conductor boundary condition literally sets the tangential magnetic field to zero which in
addition to setting the surface current density to zero also makes the tangential electric
field (and in dynamics the tangential electric field) discontinuous.
Port
Use the Port node where electromagnetic energy enters or exits the model. A port can
launch and absorb specific modes. Use the boundary condition to specify wave type
ports. Ports support S-parameter calculations but can be used just for exciting the
model. This node is not available with the Electromagnetic Waves, Transient interface.
A port assumes that the cross section’s geometry and material is constant in the normal
direction. Furthermore, the port boundary is assumed to be flat, resulting in a constant
normal across the boundary.
In 3D, the following subnodes are available from the context menu, by right-clicking
the Port node, or from the Physics toolbar, Attributes menu:
• Circular Port Reference Axis to determine a reference direction for the modes. This
subnode is selected from the Points submenu when Circular is selected as the type of
port.
• Periodic Port Reference Point to uniquely determine reciprocal lattice vectors. This
subnode is selected from the Points submenu when Periodic is selected as the type of
port.
When Transverse electromagnetic (TEM) is selected as the type of port, the following
subnodes are available from the context menu, by right-clicking the Port node, or from
the Physics toolbar, Attributes menu:
• Ground to define zero potential boundaries for the port mode field calculation in
TEM Boundary Mode Analysis.
• Electric Potential to define positive or negative 1 V potential boundaries for the port
mode field calculation in TEM Boundary Mode Analysis.
PORT PROPERTIES
Enter a unique Port name. Only nonnegative integer numbers can be used as Port name
as it is used to define the elements of the S-parameter matrix and numeric port names
are also required for port sweeps and Touchstone file export.
Periodic ports are available in 3D and 2D. Circular and Coaxial ports are available in 3D
and 2D axisymmetry.
Numeric ports require a Boundary Mode Analysis study type. It should appear before the
frequency domain study node in the study branch of the model tree. If more than one
numeric port is needed, use one Boundary Mode Analysis node per port and assign
each to the appropriate port. Then, it is best to add all the studies; Boundary Mode
Analysis 1, Boundary Mode Analysis 2, …, Frequency Domain 1, manually. Numeric
ports are by default computed for the deformed mesh whereas other types of ports
compute the mode shape using geometry information.
The Boundary Mode Analysis study step stores the frequency fref and
propagation constant βref for which it was run. For a TE, TM, or TEM
mode, the propagation constant β for an arbitrary frequency f is given by
2 2 2 2
β = β ref + k ( 1 – ( f ref ⁄ f ) ) .
In addition, for TE, TM, and TEM modes, the mode field shape is
independent of the frequency. Thus, when making a frequency sweep
including only TE, TM, and TEM modes, the Boundary Mode Analysis
study steps can be done for just one frequency, with the propagation
constants obtained from the expression above for the other frequencies.
For waveguides consisting of multiple dielectric materials, like optical
fibers, where there are no TE, TM, or, TEM modes, the Boundary Mode
Analysis steps must be recomputed for each frequency.
Transverse electromagnetic (TEM) ports require a TEM Boundary Mode Analysis study type.
It should appear before the frequency domain study node in the study branch of the
model tree. Add Ground and Electric Potential subnodes to compute the port mode
field and impedance. These subnodes are available from the context menu (right-click
Then select a Slit type — PEC-backed (the default) or Domain-backed. The PEC-backed
type makes the port on interior boundaries perform as it does on exterior boundaries.
The Domain-backed type can be combined with perfectly matched layers to absorb the
excited mode from a source port and other higher order modes.
Click Toggle Power Flow Direction button to define the power flow for the port. For an
excited port, the power flow should point in to the excited domain and for a listener
When analyzing an arbitrary shape of a port, if the port boundary consists of at least
two separate conducting parts, the port can support a (quasi-) TEM mode. In this case,
it is possible to calculate the port mode characteristic impedance from the results of a
Boundary Mode Analysis and include the impedance mismatch compared to the
user-specified reference impedance in the simulation by selecting the Analyze as a TEM
field check box. The Analyze as a TEM field requires an Integration Line for Voltage
subfeature to calculate the port mode impedance. The Integration Line for Current
subfeature is optional. The impedance is calculated using either the ratio of the voltage
and current when both Integration Line for Voltage and Integration Line for Current
are defined, or the ratio of the voltage and power when only the Integration Line for
Voltage is defined. These subfeatures are available from the context menu (right-click
the Port parent node) or from the Physics toolbar, Attributes menu.
Enter a reference Characteristic impedance Zref (SI unit: Ω). The default is 50 Ω. The
mode field on a port boundary is scaled by the ratio between the computed
characteristic impedance and user-specified Zref.
User Defined
For User defined specify the eigenmode of the incoming wave at the port. Even if Wave
excitation at this port is set to Off, the mode field should represent the incoming wave
that corresponds to the actual outgoing wave. The mode field can be entered with an
arbitrary amplitude and is normalized internally.
• Enter the components of the Electric mode field E0 (SI unit: V/m) or the Magnetic
mode field H0 (SI unit: A/m). The entered expressions must be differentiable.
• Enter the Propagation constant β (SI unit: rad/m). This is frequency dependent for
all but TEM modes and a correct frequency-dependent expression must be used.
Rectangular
For Rectangular specify a unique rectangular mode.
For 2D components, to excite the fundamental mode, select the mode type Transverse
electromagnetic (TEM), since the rectangular port represents a parallel-plate waveguide
port that can support a TEM mode. Only TE modes are possible when solving for the
out-of-plane vector component, and only TM and TEM modes are possible when
solving for the in-plane vector components. There is only a single mode number,
which is selected from a list.
Coaxial
In 2D axisymmetry, Coaxial does not support nonzero azimuthal mode number. The
Azimuthal mode number in the Physics interface should be defined as zero.
For 3D components, enter the Mode number, for example, 11 for a TE11 mode, or 01
for a TM01 mode. When Circular is selected as the type of port in 3D, the Circular Port
Reference Axis subnode is available from the context menu (right-click the parent
node) or from the Physics toolbar, Attributes menu. It defines the orientation of fields
on a port boundary.
Periodic
For Periodic specify parameters for the mode field. When Periodic is selected, the
Diffraction Order port subnode is available from the context menu (right-click the
parent node) or from the Physics toolbar, Attributes menu.
Select a Input quantity — Electric field or Magnetic field and define the mode field
amplitude for the incoming wave at the port. Even if Wave excitation at this port is set
to Off, the mode field amplitude should represent the incoming wave that corresponds
to the actual outgoing wave.
• For 2D components and if the Input quantity is set to Electric field, define the Electric
mode field amplitude. For example, for a TE wave set the x, y, and z components to
0, 0, 1. Similarly, if the Input quantity is set to Magnetic field, define the Magnetic
mode field amplitude. For a TM wave set the x, y, and z components to 0, 0, and 1.
• Define the Angle of incidence, if Wave excitation at this port is On.
For 3D components, if Wave excitation at this port is On, define the Elevation angle of
incidence and Azimuth angle of incidence. The Elevation angle of incidence α1 and Azimuth
angle of incidence α2 are used in the relations
k = k parellel + k perpendicular
The Azimuth angle of incidence is the counterclockwise rotating angle from the
primitive vector a1 around the axis built with Periodic Port Reference Point and n.
For periodic ports with hexagonal port boundaries, the definition of the vector a1 is
slightly different from the default definition. In this case, the unit cell is actually a
rhomboid, with primitive vectors pointing in other directions than the side vectors of
the hexagon. Thus, for a hexagonal periodic port, the vector a1 is defined along one
of the sides of the hexagon, and it is not one of the primitive vectors of the hexagonal
point lattice. The Azimuth angle of incidence α2 is still measured from the vector a1,
even though this vector now refers to a side vector of the hexagonal port boundary and
not a primitive vector.
For 2D components define the Angle of incidence. The Angle of incidence α is defined
by the relation
k × n = k sin αz
where k is the projection of the wave vector in the xy-plane, n is the normalized
normal vector to the boundary, k is the magnitude of the projected wave vector in the
xy-plane, and z is the unit vector in the z direction.
Notice that the mode field defined for the Periodic port assumes
homogeneous isotropic material properties in the domain adjacent to the
selected port boundary.
Default Polarization plots are automatically generated for Periodic ports in 3D and in
2D, if the Electric field components solved for setting in the Components section for the
physics interface is set to Three-component vector. The Polarization plot includes
polarization ellipses for each diffraction order. The polarization ellipse line graphs are
generated by plotting the in-plane Jones vector element versus the out-of-plane Jones
vector element.
For more information about Jones Vectors and Polarization plots, see
Jones Vectors and Polarization Plots.
• Select the Include in automatic diffraction order calculation check box to add
Diffraction Order subnodes to the selected Periodic port, when the Add Diffraction
Orders button is clicked from the exciting Periodic port.
• Define the Refractive index, real part at the boundary.
• Define the Maximum frequency — From study (the default) or User defined. When
From study is selected, the Maximum frequency is taken from the study step associated
with the physics interface. For User defined, enter the maximum frequency fmax (SI
unit: Hz). The default value is 0 Hz. If a single frequency is used, insert the
frequency, or if a frequency sweep is performed, insert the maximum frequency of
the sweep. This parameter is only available when Wave excitation at this port is On.
When all parameters are defined, click the Add Diffraction Orders button from the
exciting Periodic port to automatically create Diffraction Order ports as subnodes to all
Periodic ports having the Include in automatic diffraction order calculation check box
selected.
DE-EMBEDDING PORT
Set the Port Offset doffset (SI unit: m) for the de-embedded S-parameter calculation.
The default is 0 m. The phase of the de-embedded S-parameters is adjusted from the
calculated S-parameters with the propagation constant and the value of doffset.
where m is the source port name, n is the listener port name, β is propagation constant,
and d is offset distance from the port boundary.
CONSTRAINT SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. For information about the
Constraint Settings section, see Constraint Settings in the COMSOL Multiphysics
Reference Manual.
Ground
The Ground subnode is available from the context menu (right-click the Port node)
when Transverse electromagnetic (TEM) is selected as the type of port. Zero potential is
assigned on the boundary selection in the port mode field calculation.
Electric Potential
The Electric Potential subnode is available from the context menu (right-click the Port
node) when Transverse electromagnetic (TEM) is selected as the type of port. Either
positive or negative 1 V is assigned on the boundary selection in the port mode field
calculation.
Diffraction Order
The Diffraction Order port is available in 3D and 2D. When the Type of Port is set to
Periodic under Port Properties, this subnode is available from the context menu
(right-click the Port parent node) or from the Physics toolbar, Attributes menu.
Use the Diffraction Order port to define diffraction orders from a periodic structure.
Normally a Diffraction Order node is added automatically during the Periodic port
setup. Additional Diffraction Order ports subnodes are available from the context menu
(right-click the parent node) or from the Physics toolbar, Attributes menu.
PORT PROPERTIES
Enter a unique Port name. Only nonnegative integer numbers can be used as Port name
as it is used to define the elements of the S-parameter matrix and numeric port names
are also required for port sweeps and Touchstone file export.
Components
Select the Components setting for the port — In-plane vector (the default) or
Out-of-plane vector.
Note that In-plane vector and Out-of-plane vector are based on the plane of diffraction
which is constructed with the diffraction wave vector and the outward normal vector
of the port boundary. The diffraction wave vector is defined by
k diffraction,parallel = k F + MG 1 + NG 2
2 2
k diffraction,perpendicular = k – k diffraction,parallel
In-plane vector lies on the plane of diffraction while Out-of-plane vector is normal to the
plane of diffraction.
For a 2D component, In-plane vector is available when the settings for the physics
interface is set to either In-plane vector or Three-component vector under Components.
Out-of-plane vector is available when the settings for the physics interface is set to either
Out-of-plane vector or Three-component vector under Components.
Enter a value or expression for the Mode phase θin (SI unit: rad). The default is
0 radians. The Mode phase setting is further discussed for the Port feature.
Notice that the mode field defined for the Diffraction Order feature
assumes isotropic material properties in the domain adjacent to the
selected feature boundary.
Use the Orthogonal Polarization port to define port with a mode that is orthogonal to
the mode of the parent Periodic Port. Normally a Orthogonal Polarization node is added
automatically during the Periodic port setup, but the Orthogonal Polarization port
subnode is also available from the context menu (right-click the parent node) or from
the Physics toolbar, Attributes menu. Only one Orthogonal Polarization node can be
added per parent Periodic Port.
PORT PROPERTIES
Enter a unique Port name. Only nonnegative integer numbers can be used as Port name
as it is used to define the elements of the S-parameter matrix and numeric port names
are also required for port sweeps and Touchstone file export.
Diffraction Order
The Orthogonal Polarization port represent a zero-order mode (same as assumed for the
parent Periodic port). Thus, the Diffraction order settings are just for information and
cannot be edited.
Enter a value or expression for the Mode phase θin (SI unit: rad). The default is
0 radians. The Mode phase setting is further discussed for the Port feature.
Notice that the mode field defined for the Orthogonal Polarization feature
assumes isotropic material properties in the domain adjacent to the
selected feature boundary.
The Periodic Port Reference Point is used to uniquely identify two primitive unit cell
vectors, a1 and a2, and two reciprocal lattice vectors, G1 and G2. These reciprocal
vectors are defined in terms of the unit cell vectors, a1 and a2, tangent to the edges
shared between the port and the adjacent periodic boundary conditions. G1 and G2
are defined by the relation
a1 × a2
---------------------
- = n
a1 × a2
a2 × n n × a1
G 1 = 2π --------------------------- and G 2 = 2π ---------------------------
a1 ⋅ a2 × n a1 ⋅ a2 × n
where n is the outward unit normal vector to the port boundary. If there are multiple
points defined in the selection list, only the last point is used.
POINT SELECTION
The primitive unit cell vectors, a1 and a2 are defined from two edges sharing the
Periodic Port Reference Point on a port boundary. The two vectors can have unequal
lengths and are not necessarily orthogonal. They start from the Periodic Port Reference
Point.
For listener (passive, observation, and not excited) ports, if the outward normal vector
on the listener port boundary is opposite to that of the source port, the listener port
reference point needs to be mirrored from the source port reference point based on
the center coordinate of the model domain. For example, if the source port reference
point is at {-1,-1,1} in a cubic domain around the origin, the mirrored listener port
reference point is {1,1,-1}. In this case, if the Azimuth angle of incidence at the source
port boundary is α2, the Azimuth angle of incidence at the listener port boundary π/
2 − α2 and the signs of the diffraction order on the source and listener ports are
opposite. See also Periodic for the angle definition.
If the lattice vectors are collinear with two Cartesian axes, then the lattice vectors can
be defined without the Periodic Port Reference Point. For the port where n points along
a positive Cartesian direction, a1 and a2 are also assigned to point along positive
Cartesian directions. Conversely, for the port where n points along a negative
Lumped Port
Use the Lumped Port node to apply a voltage or current excitation of a model or to
connect to a circuit. A lumped port is a simplification of the port boundary condition.
A Lumped Port condition can only be applied on boundaries that extend between two
metallic boundaries — that is, boundaries where Perfect Electric Conductor, Impedance
Boundary, or Transition Boundary (Electromagnetic Waves, Frequency Domain
interface only) conditions apply — separated by a distance much smaller than the
wavelength.
Select User defined for nonuniform ports, for example, a curved port and enter values
or expressions in the fields — Height of lumped port hport (SI unit: m), Width of lumped
port wport (SI unit: m), and Direction between lumped port terminals ah. In 2D
axisymmetry, Coaxial does not support a nonzero azimuthal mode number. The
Azimuthal mode number should be set to zero in the Out-of-Plane Wave Number section
Notice that the input field for Direction between lumped port terminals ah
is not shown in 2D, when Electric field components solved for is set to
Out-of-plane vector in the Components section for the physics interface. In
this case the Direction between lumped port terminals ah is defined to be
in the out-of-plane direction.
Terminal Type
Select a Terminal type — a Cable port for a voltage driven transmission line and
S-parameter calculation, a Current driven port, or a Circuit port.
For Cable select On or Off from the Wave excitation at this port list to set whether it is
an inport (excitation) or a listener port (observation).
Source Type
If On is selected for the Wave excitation at this port, select Voltage or Power from the
Source type. For the Voltage source type, enter a Voltage V0 (SI unit: V), and Port phase
θin (SI unit: rad). For the Power source type, enter a Power P0 (SI unit: W) that is the
average input power to the lumped port.
Note it is only possible to excite one Cable port at a time if the purpose is to compute
S-parameters. In other cases, for example, when studying microwave heating, more
than one inport might be wanted, but the S-parameter variables cannot be correctly
computed so if several ports are excited, the S-parameter output is turned off.
• For a Cable terminal type enter the Characteristic impedance Zref (SI unit: Ω).
• For a Current terminal type enter a Terminal current I0 (SI unit: A).
CALCULATOR
Enter a Relative permittivity (default is 1) for the dielectric material between the coaxial
inner and outer conductor and then click the Compute Coaxial Line Impedance button
to compute the impedance for the coaxial lumped port. The result is displayed in the
Messages window.
CALCULATE S-PARAMETER
The Calculate S-parameter check box needs to be activated for S-parameter calculation
with the Electromagnetic Waves, Transient and Electromagnetic Waves, Time Explicit
interfaces, while the Cable port in the Electromagnetic Wave, Frequency Domain
calculates S-parameters automatically.
( x – 2 ⁄ f0 )
2
2f 0 – --------------------------
2 ⁄ ( 2f 0 )
2
-
----------
-e ⋅ sin ( 2π ⋅ ( 1 + η f ) ⋅ f 0 )
2π
where the Center frequency f0 defines the location as 2/f0, the standard deviation as 1/
(2f0), and the modulation frequency f0. ηf is the upshift ratio for the sinusoidal
modulation frequency.
The modulating sinusoidal function can be slightly shifted from the center frequency
f0 by a factor of 1+ηf to improve the frequency responses.
CENTER FREQUENCY
When Calculate S-parameter is checked, enter a Center frequency f0 for the input
modulated Gaussian pulse (SI unit: Hz).
Lumped Element
Use a Lumped Element node to mimic the insertion of a capacitor, inductor, or general
impedance between two metallic boundaries. A Lumped Element condition is a passive
lumped port boundary condition which cannot be used as a source. Unlike a Lumped
Port, it does not generate S-parameters. The sign of the current and power of a Lumped
Element is opposite to that of a Lumped Port.
It can only be applied on boundaries that extend between two metallic boundaries —
that is, boundaries where Perfect Electric Conductor, Impedance Boundary, or Transition
Boundary (Electromagnetic Waves, Frequency Domain interface only) conditions
apply — separated by a distance much smaller than the wavelength.
SETTINGS
Select a Lumped element device — User defined (the default), Inductor, Capacitor,
Parallel LC, Series LC, Parallel RLC, or Series RLC. Inductor, Capacitor, Parallel LC, Series LC,
Parallel RLC, and Series RLC are available only in the frequency domain study type.
• For User defined enter a Lumped element impedance Zelement (SI unit: Ω). The
default is 50 Ω.
CALCULATOR
Enter a Relative permittivity (default is 1) for the dielectric material between the coaxial
inner and outer conductor and then click the Compute Coaxial Line Impedance button
to compute the impedance for the coaxial lumped element. The result is displayed in
the Messages window.
Uniform Element
The Uniform Element is available in 3D. When the Type of lumped port is set to
Multielement uniform under Lumped Port Properties, this subnode is available from the
context menu (right-click the Lumped Port parent node) or from the Physics toolbar,
Attributes menu.
PORT PROPERTIES
Enter a unique Uniform element name. Enter values or expressions in the fields —
Direction between uniform element terminals ah. The Direction between uniform element
terminals defines the electric potential polarity, as ah points in the direction of the
electric field.
Electric Field
The Electric Field boundary condition
n × E = n × E0
specifies the tangential component of the electric field. It should in general not be used
to excite a model. Consider using the Port, Lumped Port, or Scattering Boundary
Condition instead. It is provided mainly for completeness and for advanced users who
ELECTRIC FIELD
Enter the value or expression for the components of the Electric field E0 (SI unit: V/
m).
CONSTRAINT SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. For information about the
Constraint Settings section, see Constraint Settings in the COMSOL Multiphysics
Reference Manual.
Magnetic Field
The Magnetic Field node adds a boundary condition for specifying the tangential
component of the magnetic field at the boundary:
n × H = n × H0
MAGNETIC FIELD
Enter the value or expression for the components of the Magnetic field H0 (SI unit: A/
m).
– jk ( n ⋅ r ) – jk ( k ⋅ r )
E = E sc e + E0 e Plane scattered wave
– jk ( n ⋅ r )
e – jk ( k ⋅ r )
E = E sc ------------------------ + E 0 e Cylindrical scattered wave
r
– jk ( n ⋅ r )
e – jk ( k ⋅ r )
E = E sc ------------------------ + E 0 e Spherical scattered wave
rs
– jk ( k ⋅ ( r – r ref ) )
E0 e ,
where rref is a reference point determined as the average point from the point selection
in the Reference Point subnode.
The boundary is only perfectly transparent for scattered (outgoing) waves of the
selected type at normal incidence to the boundary. That is, a plane wave at oblique
incidence is partially reflected and so is a cylindrical wave or spherical wave unless the
wave fronts are parallel to the boundary. For the Electromagnetic Waves, Frequency
Domain interface, the Perfectly Matched Layer feature is available as a general way of
modeling an open boundary.
• For cylindrical waves, specify around which cylinder axis the waves are cylindrical.
Do this by specifying one point at the cylinder axis and the axis direction.
• For spherical waves, specify the center of the sphere around which the wave is
spherical.
The domain material adjacent to the boundary where the Scattering Boundary
Condition is applied can be lossy.
If the problem is solved for the eigenfrequency or the scattered field, the boundary
condition does not include the incident wave.
– jk ( n ⋅ r )
E sc = E sc e Plane scattered wave
– jk ( n ⋅ r )
e
E sc = E sc ------------------------ Cylindrical scattered wave
r
– jk ( n ⋅ r )
e
E sc = E sc ------------------------ Spherical scattered wave
rs
If the Incident field is set to Gaussian beam, select an Input quantity: Electric field
amplitude (the default) or Power. If the Input quantity is Electric field amplitude, enter
the component expressions for the Gaussian beam electric field amplitude Eg0 (SI unit:
V/m). If the Input quantity is set to Power, enter the Input power (SI unit: W in 2D
axisymmetry and 3D and W/m in 2D) and the component expressions for the Gaussian
beam non-normalized electric field amplitude Eg0 (SI unit: V/m). Also edit the Beam
radius w0 (SI unit: m) and the Distance to focal plane p0 (SI unit: m). The default values
are ((10*2)*pi)/emw.k0 and 0 m, respectively. The optical axis for the Gaussian
beam is defined by a line including a reference point on the feature selection with a
direction specified by the Incident wave direction (see below). By default, the reference
point is the average position for the feature selection. However, by adding a Reference
Point subnode any available point (or the average of several selected points) on the
feature selection can be used as the reference point. The focal plane for the Gaussian
beam is located the Distance to focal plane p0 from the reference point in the Incident
wave direction.
If the Incident field is not set to No incident field, edit the Incident wave direction kdir
for the vector coordinates. The default direction is in the opposite direction to the
boundary normal. For 2D axisymmetry, the Incident wave direction kdir should be
parallel or anti-parallel to the symmetry axis.
Select a Scattered wave type for which the boundary is absorbing — Plane wave (the
default), Spherical wave, or Cylindrical wave.
• For the Electromagnetic Waves, Frequency Domain interface, select an Order —First
order (the default) or Second order.
• For Cylindrical wave also enter coordinates for the Source point r0 (SI unit: m) and
Source axis direction raxis (dimensionless). For 2D the Source axis direction is
assumed to be in the z direction, whereas in 2D axisymmetry it is assumed to be
along the axis of rotation.
• For Spherical wave enter coordinates for the Source point r0 (SI unit: m).
MODE ANALYSIS
Expand the Mode Analysis section and check the Subtract propagation constant from
material wave number check box to calculate the wave number for the scattered wave as
where kn is the wave number for the scattered wave propagating in the normal
direction, k is the material wave number, and β is the propagation constant,
determined from the mode analysis. If the check box is unchecked (the default),
kn = k.
For more information about the Gaussian beam theory, see Gaussian
Beams as Background Fields and Input Fields.
Select the Dispersion and absorption model that will be used when calculating the wave
number and attenuation constant for the incident and scattered waves — Low loss
When the Dispersion and absorption model is set to Low loss approximation
the refractive index is calculated from the relative permittivity and the
relative permeability as
n = εr μr .
1 μ0 μr 1
γ = --- σ -----------
- = --- σZ c ,
2 ε0 εr 2
When the Dispersion and absorption model is set to High loss, the real and
the imaginary parts of the complex refractive index is solved for from the
real and the imaginary parts of the relative permittivity, using the relations
2 2
n – κ = ε' r μ r
and
σμ r
2nκ = ε'' r μ r = --------- .
ωε 0
ω
γ = ---- κ .
c
Reference Point
The Reference Point subnode is available only when there is an available incident field
defined in the parent node. Then this subnode is available from the context menu
(right-click the Scattering Boundary Condition parent node) or from the Physics
toolbar, Attributes menu.
In the parent node, the incident field is then defined using the reference position:
– jk ( k ⋅ ( r – r ref ) )
E0 e .
POINT SELECTION
Select the points that should be used when calculating the reference position. The
reference position is calculated as the average position of the selected points. The point
selection is only effective when Definition in the Reference Point section is set to Point
selection.
REFERENCE POINT
Select the Definition for the reference point — Point selection (the default) or User
defined. When User defined is selected, enter the expressions for the components for the
Reference point r0. The Reference point must be a point on the parent feature’s
boundary selection.
μ0 μr
------------ n × H + E – ( n ⋅ E )n = ( n ⋅ E s )n – E s
εc
is used at boundaries where the field is known to penetrate only a short distance
outside the boundary. This penetration is approximated by a boundary condition to
avoid the need to include another domain in the model. Although the equation is
identical to the one in the low-reflecting boundary condition, it has a different
interpretation. The material properties are for the domain outside the boundary and
not inside, as for low-reflecting boundaries. A requirement for this boundary condition
to be a valid approximation is that the magnitude of the complex refractive index
με c
N = -----------
-
μ1 ε1
where μ1 and ε1 are the material properties of the inner domain, is large; that is,
|N| >> 1.
The Surface Roughness subnode is available from the context menu (right-click the
parent node) or from the Physics toolbar, Attributes menu.
Js
J
I'
I
Js
Enter a Source electric field Es (SI unit: V/m). The default is 0 V/m.
–n × H = Js
n2 × ( H1 – H2 ) = Js
specifies a surface current density at both exterior and interior boundaries, respectively.
The current density is specified as a three-dimensional vector, but because it needs to
flow along the boundary surface, COMSOL Multiphysics projects it onto the
boundary surface and neglects its normal component. This makes it easier to specify
the current density and avoids unexpected results when a current density with a
component normal to the surface is given.
For the Surface Current Density subnode, select Side — Upside (the default) or Downside
to define on which side the Surface Current Density is applied. The red arrow visualized
on the selected boundaries always indicates the upside.
n × E = J ms
n 2 × ( E 1 – E 2 ) = – J ms
specifies a surface magnetic current density at both exterior and interior boundaries,
respectively. The magnetic current density is specified as a three-dimensional vector,
but because it needs to flow along the boundary surface, COMSOL Multiphysics
projects it onto the boundary surface and neglects its normal component. This makes
it easier to specify the magnetic current density and avoids unexpected results when a
magnetic current density with a component normal to the surface is given.
CONSTRAINT SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. For information about the
Constraint Settings section, see Constraint Settings in the COMSOL Multiphysics
Reference Manual.
Surface Roughness
The Surface Roughness subnode is available from the context menu (right-click the
Impedance Boundary Condition or Transition Boundary Condition parent node) or
from the Physics toolbar, Attributes menu.
SURFACE ROUGHNESS
Select a Surface roughness model — Sawtooth (the default) or Snowball. For Sawtooth,
enter a value or expression for the Surface roughness ΔRMS (SI unit: m). For Snowball,
enter a Snowball hexagon cell area Ahex (SI unit: m2). Then enter values for the Snowball
radius (SI unit: m) and Number of snowballs in the table, adding as many rows as you
need.
Select the Skin depth type to be Physics-controlled (the default) or User defined. For User
defined, enter the Skin depth δ (SI unit: m).
The Surface Roughness increases the impedance, and consequently it decreases the
surface current density of Impedance Boundary Condition or Transition Boundary
Condition, proportional to the ratio between the impedance of a rough surface and
that of a smooth surface:
SAWTOOTH MODEL
Δ RMS 2
Z rough = Z smooth ⋅ 1 + --- atan ( 1.4 ) ---------------
2
π δ
where ΔRMS is the root mean square of the surface roughness, and δ is the skip depth
of the material (Ref. 1).
2
3 N i ⋅ 4πa i δ δ
2
Z rough = Z smooth ⋅ 1 + ---
2 ------------------------ ⁄ 1 + ----- + -----2-
A hex ai a
i i
where Ahex is the hexagonal area of a unit cell, Ni is the number of snowballs, ai is the
radius of a snowball, and δ is the skip depth of the material (Ref. 2).
For Transition Boundary Condition, select Side — Upside (the default) or Downside to
define on which side the Surface Roughness is applied. The red arrow visualized on the
selected boundaries always indicates the upside.
References
1. E. Hammerstad, O. Jensen, “Accurate Models for Microstrip Computer-Aided
Design”, Microwave symposium Digest, 1980 IEEE MTT-S International, pp. 407–
409, May 1980
( Z S E t1 – Z T E t2 )
J s1 = --------------------------------------------
2 2
-
ZS – ZT
( Z S E t2 – Z T E t1 )
J s2 = --------------------------------------------
2 2
-
ZS – ZT
– jωμ 1
Z S = ------------- ----------------------
k tan ( kd )
– jωμ 1
Z T = ------------- ---------------------
k sin ( kd )
k = ω ( ε + ( σ ⁄ ( jω ) ) )μ
The thickness of the layer should also be less than the radius of curvature
for the boundary.
The Surface Roughness subnode is available from the context menu (right-click the
parent node) or from the Physics toolbar, Attributes menu.
Select the Electrically thick layer check box (unselected by default) to make the two
domains adjacent to the boundary uncoupled. Use this setting, for instance, when the
thickness is greater than three times of the skin depth. When the Electrically thick layer
check box is unselected, enter a Thickness d (SI unit: m). The default is 0.01 m.
Periodic Condition
The Periodic Condition sets up a periodicity between the selected boundaries. The
Destination Selection subnode is available from the context menu (right-click the parent
node) or from the Physics toolbar, Attributes menu.
PERIODICITY SETTINGS
Select a Type of periodicity — Continuity (the default), Antiperiodicity, or Floquet
periodicity. Select:
• Continuity to make the electric field periodic (equal on the source and destination),
• Antiperiodicity to make it antiperiodic, or
• Floquet periodicity (The Electromagnetic Waves, Frequency Domain Interface only)
to use a Floquet periodicity (Bloch-Floquet periodicity).
- For Floquet periodicity also enter the source for the k-vector for Floquet periodicity.
- For User defined specify the components of the k-vector for Floquet periodicity kF
(SI unit: rad/m).
- For From periodic port the k-vector for Floquet periodicity kF is obtained from the
Periodic Port settings.
CONSTRAINT SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. For information about the
Constraint Settings section, see Constraint Settings in the COMSOL Multiphysics
Reference Manual.
ORIENTATION OF SOURCE
For information about the Orientation of Source section, see Orientation of Source and
Destination in the COMSOL Multiphysics Reference Manual.
Magnetic Current
The Magnetic Current node specifies a magnetic line current along one or more edges.
For a single Magnetic Current source, the electric field is orthogonal to both the line
and the distance vector from the line to the field point. For 2D and 2D axisymmetric
models the Magnetic Current node is applied to Points, representing magnetic currents
directed out of the model plane. For 3D models, the Magnetic Current is applied to
Edges.
MAGNETIC CURRENT
Enter a value for the Magnetic current Im (SI unit: V).
Two-Port Network
Use the Two-Port Network node to characterize the response of a two-port network
system such as reflection and transmission using S-parameters.
A Two-Port Network can only be applied on boundaries that extend between two
metallic boundaries — that is, boundaries where Perfect Electric Conductor, Impedance
Boundary, or Transition Boundary conditions apply — separated by a distance much
smaller than the wavelength.
Type of Port
Select a Type of Port — Coaxial or Uniform.
Coaxial does not support nonzero azimuthal mode number. The Azimuthal mode
number in the Physics interface should be defined as zero.
SETTINGS
Enter the Characteristic impedance Zref (SI unit: Ω) for Matrix.
Interpolation
Select a Interpolation — Nearest neighbor, Linear, Piecewise cubic (the default), or Cubic
spline.
The S-parameter values are interpolated within the frequency range specified in the
Touchstone file.
Extrapolation
Select a Extrapolation — Constant (the default), Linear, Nearest function, or Specific
value.
The S-parameter values are extrapolated outside the frequency range specified in the
Touchstone file.
BALANCED PORT
Specify Port name for port a in balanced pair 1, Port name for port b in balanced pair 1,
Port name for port c in balanced pair 2, and Port name for port d in balanced pair 2.
Edge Current
The Edge Current node specifies an electric line current along one or more edges.
EDGE CURRENT
Enter an Edge current I0 (SI unit: A).
DIPOLE SPECIFICATION
Select a Dipole specification — Magnitude and direction or Dipole moment.
DIPOLE PARAMETERS
Based on the Dipole specification selection:
• For Magnitude and direction enter coordinates for the Electric current dipole moment
direction np and Electric current dipole moment, magnitude p (SI unit: A·m).
• For Dipole moment enter coordinates for the Electric current dipole moment p (SI
unit: A·m).
DIPOLE SPECIFICATION
Select a Dipole specification — Magnitude and direction or Dipole moment.
• For Magnitude and direction enter coordinates for the Magnetic dipole moment
direction nm and Magnetic dipole moment, magnitude m (SI unit: m2·A).
• For Dipole moment enter coordinates for the Magnetic dipole moment m (SI unit:
m2·A).
Archie’s Law
This subfeature is available only when Archie’s law is selected as the Electrical
conductivity material parameter in the parent feature (for example, the Wave Equation,
Electric node). Then the subnodes are made available from the context menu
(right-click the parent node) as well as from the Physics toolbar, Attributes menu.
Use the Archie’s Law subnode to provide an electrical conductivity computed using
Archie’s Law. This subnode can be used to model nonconductive effective medium
saturated (or variably saturated) by conductive liquids, using the relation:
n m
σ = sL εp σL
CONDUCTION CURRENTS
By default, the Electrical conductivity σL (SI unit: S/m) for the fluid is defined From
material. This uses the value of the conductivity of the material domain.
• Cementation exponent m
• Saturation exponent n
• Fluid saturation SL
• Porosity εp to set up the volume fraction of the fluid.
Effective Medium
This subfeature is available only when Effective medium is selected as the material
parameter (for example, Relative permeability or Relative permittivity) in the parent
feature node when it is available with the physics interface (for example, the Wave
Equation, Electric node). Then the subnodes are made available from the context menu
(right-click the parent node) as well as from the Physics toolbar, Attributes menu.
Use the Effective Medium subfeature to specify the material properties of a domain
consisting of a porous medium using a mixture model. Depending on the specific
physics interface being used, the subfeature can be used to provide a mixture model
for the electric conductivity σ, the relative dielectric permittivity εr, or the relative
magnetic permeability μr.
EFFECTIVE MEDIUM
This section is always available and is used to define the mixture model for the domain.
For each material (Material 1, Material 2, and so on), select either Domain material, to
use the material specified for the domain, or one of the other materials specified in the
Materials node. For each material, enter a Volume fraction θ1, θ2, and so on.
When this physics interface is added, these default nodes are also added to the Model
Builder — Wave Equation, Electric, Perfect Electric Conductor, and Initial Values. Then,
from the Physics toolbar, add other nodes that implement, for example, boundary
conditions and mass sources. You can also right-click Electromagnetic Waves, Transient
to select physics features from the context menu.
Except where indicated, most of the settings are the same as for The Electromagnetic
Waves, Frequency Domain Interface.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is temw.
COMPONENTS
This section is available for 2D and 2D axisymmetric components.
DEPENDENT VARIABLES
The dependent variable (field variable) is for the Magnetic vector potential A. The name
can be changed but the names of fields and dependent variables must be unique within
a model.
DISCRETIZATION
Select the shape order for the Magnetic vector potential dependent variable — Linear,
Quadratic (the default), or Cubic. For more information about the Discretization
section, see Settings for the Discretization Sections in the COMSOL Multiphysics
Reference Manual.
The domain, boundary, edge, point, and pair nodes are available from the Physics
ribbon toolbar (Windows users), Physics context menu (Mac or Linux users), or
right-click to access the context menu (all users).
BOUNDARY CONDITIONS
With no surface currents present the boundary conditions
n2 × ( E1 – E2 ) = 0
n2 × ( H1 – H2 ) = 0
need to be fulfilled. Depending on the field being solved for, it is necessary to analyze
these conditions differently. When solving for A, the first condition can be formulated
in the following way.
∂A 2 ∂A 1 ∂
n2 × ( E1 – E2 ) = n2 × – = ( n2 × ( A2 – A1 ) )
∂t ∂t ∂t
The tangential component of the magnetic vector potential is always continuous and
thus the first condition is fulfilled. The second condition is equivalent to the natural
boundary condition.
–1 –1 –1
–n × ( μr ∇ × A1 – μr ∇ × A2 ) = –n × μr ( H1 – H2 ) = 0
These nodes and subnodes are available and described for the Electromagnetic Waves,
Frequency Domain interface (listed in alphabetical order):
In the COMSOL Multiphysics Reference Manual see Table 2-4 for links
to common sections and Table 2-5 to common feature nodes. You can
also search for information: press F1 to open the Help window or Ctrl+F1
to open the Documentation window.
∂A ∂ ∂A
+ μ 0 ε 0 ε r + ∇ × ( μ r ∇ × A ) = 0
–1
μ0 σ
∂t ∂t ∂t
for transient problems with the constitutive relations B = μ0μrH and D = ε0εrE. Other
constitutive relations can also be handled for transient problems. The Divergence
Constraint subnode is available from the context menu (right-click the parent node) or
from the Physics toolbar, Attributes menu.
Relative Permittivity
When Relative permittivity is selected, the default Relative permittivity εr
(dimensionless) takes values From material. For User defined select Isotropic, Diagonal,
Symmetric, or Full and enter values or expressions in the field or matrix. If Effective
medium is selected, the Effective medium subnode is available from the context menu
(right-click the parent node) or from the Physics toolbar, Attributes menu.
Refractive Index
When Refractive index is selected, the default Refractive index n (dimensionless) takes
the value From material. To specify the refractive index and assume a relative
permeability of unity and zero conductivity, for one or both of the options, select User
defined then choose Isotropic, Diagonal, Symmetric, or Full. Enter values or expressions
in the field or matrix.
Notice that only the real part of the refractive index is used for the
transient formulation.
Polarization
For Polarization enter coordinates for the Polarization P (SI unit: C/m2).
D = ε0 ε∞ E + Pn ,
n=1
∂2 ∂ 2 2
-------2- + Γ n ----- + ω n P n = ε 0 f n ω p E .
∂t ∂t
MAGNETIC FIELD
This section is available if Relative permittivity, Polarization, or Remanent electric
displacement are chosen as the Electric displacement field model.
Select the Constitutive relation — Relative permeability (the default), Remanent flux
density, or Magnetization.
Relative Permeability
For Relative permeability the relative permeability μr uses values From material. For User
defined select Isotropic, Diagonal, Symmetric, or Full based on the characteristics of the
magnetic field, and then enter values or expressions in the field or matrix. If Effective
medium is selected, the Effective medium subnode is available from the context menu
(right-click the parent node) or from the Physics toolbar, Attributes menu.
Magnetization
For Magnetization enter coordinates for M (SI unit: A/m).
CONDUCTION CURRENT
This section is available if Relative permittivity, Polarization, or Remanent electric
displacement are chosen as the Electric displacement field model.
By default, the Electrical conductivity σ (SI unit: S/m) uses values From material.
• For User defined select Isotropic, Diagonal, Symmetric, or Full based on the
characteristics of the current and enter values or expressions in the field or matrix.
Initial Values
The Initial Values node adds an initial value for the magnetic vector potential and its
time derivative that serves as initial conditions for the transient simulation.
INITIAL VALUES
Enter values or expressions for the initial values of the components of the magnetic
vector potential A (SI unit: Wb/m) and its time derivative ∂A/∂t (SI unit: V/m). The
default values are 0 Wb/m and 0 V/m, respectively.
Drude–Lorentz Polarization
This subfeature is available only when Drude-Lorentz Dispersion Model is selected as the
Electric displacement field model in the Wave Equation, Electric feature node. Then the
subnodes are made available from the context menu (right-click the parent node) as
well as from the Physics toolbar, Attributes menu.
D = ε0 ε∞ E + Pn ,
n=1
∂2 ∂ 2 2
-------2- + Γ n ----- + ω n P n = ε 0 f n ω p E .
∂t ∂t
Enter values or expressions for the Oscillator strength fn (SI unit: 1), the Resonance
frequency ωn (SI unit: rad/s), and the Damping in time coefficient Γn (SI unit: rad/s).
INITIAL VALUES
Enter values or expressions for the initial values of the components of the Drude–
Lorentz polarization Pn (SI unit: C/m2) and its time derivative ∂Pn/∂t (SI unit: A/
m2).
DISCRETIZATION
To display this section, click the Show More Options button ( ) and select
Discretization from the Show More Options dialog box. Select the element order from
the list box for the Drude-Lorentz polarization Pn.
The physics interface is used when solving for electromagnetic wave propagation along
one-dimensional transmission lines and is available in 1D, 2D, and 3D. The physics
interface has Eigenfrequency and Frequency Domain study types available. The
Frequency Domain study is used for source driven simulations for a single frequency
or a sequence of frequencies.
When this physics interface is added, these default nodes are also added to the Model
Builder — Transmission Line Equation, Absorbing Boundary, and Initial Values. Then,
from the Physics toolbar, add other nodes that implement, for example, boundary
conditions. You can also right-click Transmission Line to select physics features from the
context menu.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is tl.
DEPENDENT VARIABLES
The dependent variable (field variable) is the Electric potential V (SI unit: V). The name
can be changed but the names of fields and dependent variables must be unique within
a model.
DISCRETIZATION
Select the shape order for the Electric potential dependent variable — Linear, Quadratic
(the default), or Cubic. For more information about the Discretization section, see
Settings for the Discretization Sections in the COMSOL Multiphysics Reference
Manual.
• Domain, Boundary, Edge, Point, and Pair Nodes for the Transmission
Line Equation Interface
• Theory for the Transmission Line Interface
• Visualization and Selection Tools in the COMSOL Multiphysics
Reference Manual
For all space dimensions, select Points for the boundary condition.
1 - ∂V
∂ -------------------- – ( G + iωC )V = 0
∂ x R + iωL ∂ x
The default values give a characteristic impedance for the transmission line of 50 Ω.
INITIAL VALUES
Enter values or expressions for the initial values of the Electric potential V (SI unit: V).
Absorbing Boundary
The Absorbing Boundary condition is stated as
n ⋅ ∇V - ------ V
-------------------- + = 0
R + jωL Z 0
γ = ( R + iωL ) ( G + iωC )
and n is the normal pointing out of the domain. The absorbing boundary condition
prescribes that propagating waves are absorbed at the boundary and, thus, that there
is no reflection at the boundary. The Absorbing Boundary condition is only available
on external boundaries.
Incoming Wave
The Incoming Wave boundary condition
V – 2V 0
n ⋅ ∇V - --------------------
-------------------- + = 0
R + jωL Z0
lets a wave of complex amplitude Vin enter the domain. The complex propagation
constant γ and the outward-pointing normal n are defined in the section describing
the Absorbing Boundary node. The Incoming Wave boundary condition is only
available on external boundaries.
Open Circuit
The Open Circuit boundary condition is a special case of the Terminating Impedance
boundary condition, assuming an infinite impedance, and, thus, zero current at the
boundary. The condition is thus
n ⋅ ∇V = 0
Terminating Impedance
The Terminating Impedance boundary condition
n ⋅ ∇V - ------ V
-------------------- + - = 0
R + jωL Z L
specifies the terminating impedance to be ZL. Notice that the Absorbing Boundary
condition is a special case of this boundary condition for the case when
R + jωL-
ZL = Z0 = ---------------------
G + jωC
The Open Circuit and Short Circuit boundary conditions are also special cases of this
condition. The Terminating Impedance boundary condition is only available on
external boundaries.
Short Circuit
The Short Circuit node is a special case of the Terminating Impedance boundary
condition, assuming that impedance is zero and, thus, the electric potential is zero.
The constraint at this boundary is, thus, V = 0.
CONSTRAINT SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
Lumped Port
Use the Lumped Port node to apply a voltage or current excitation of a model or to
connect to a circuit. The Lumped Port node also defines S-parameters (reflection and
transmission coefficients) that can be used in later postprocessing steps.
PORT PROPERTIES
Enter a unique Port Name. It is recommended to use a numeric name as it is used to
define the elements of the S-parameter matrix and numeric port names are also
required for port sweeps and Touchstone file export.
SETTINGS
If a Circuit port type is selected under Port Properties, this section does not
require any selection.
If Cable is selected as the port type, select the Wave excitation at this port check box to
enter values or expressions for the:
This physics interface solves two first-order partial differential equations (Faraday’s law
and Maxwell-Ampère’s law) for the electric and magnetic fields using the time explicit
discontinuous Galerkin method.
When this physics interface is added, these default nodes are also added to the Model
Builder — Wave Equations, Perfect Electric Conductor, and Initial Values. Then, from the
Physics toolbar, add other nodes that implement, for example, boundary conditions.
You can also right-click Electromagnetic Waves, Time Explicit to select physics features
from the context menu.
The interface includes absorbing layers that are used to set up effective nonreflecting
like boundary conditions. These features are added from the Definitions toolbar, by
clicking Absorbing Layer. If COMSOL Multiphysics is not running in full-screen mode
nor in a large window, Absorbing Layer is accessible on the Definitions toolbar by first
clicking Coordinate Systems and then Absorbing Layer. You can also right-click
Definitions in the Model Builder and select Absorbing Layer from the context menu.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is ewte.
COMPONENTS
This section is available for 2D and 2D axisymmetric components.
• Full wave (the default) to solve using a full three-component vector for the electric
field E and the magnetic field H.
• E in plane (TM wave) to solve for the electric field vector components in the modeling
plane and one magnetic field vector component perpendicular to the plane,
assuming that there is no electric field perpendicular to the plane and no magnetic
field components in the plane.
• H in plane (TE wave) to solve for the magnetic field vector components in the
modeling plane and one electric field vector component perpendicular to the plane.
DEPENDENT VARIABLES
The dependent variables (field variables) are for the Electric field vector E and for the
Magnetic field vector H. The name can be changed but the names of fields and
dependent variables must be unique within a model.
DISCRETIZATION
Select the shape order for the Electric and magnetic fields dependent variables (the same
order for both fields) — Linear, Quadratic, or Cubic (the default). For more information
about the Discretization section, see Settings for the Discretization Sections in the
COMSOL Multiphysics Reference Manual.
Wave Equations
The Wave Equations node is the main node for the Electromagnetic Waves, Time
Explicit interface. The governing transient equations can be written in the form
∂D
∇ × H = σE + -------
∂t
∂------
B-
∇×E = –
∂t
∂E
ε 0 ε r ------- – ∇ × H + σE = 0
∂t
∂H
μ 0 μ r -------- + ∇ × E = 0
∂t
NUMERICAL PARAMETERS
The defaults for each parameter are as follows:
• Lax-Friedrichs flux parameter for E field τE (SI unit: S), the default is 0.5/Z for
Ampere’s law.
• Lax-Friedrichs flux parameter for H fieldτH (SI unit: Ω), the default is 0.5 Z for
Faraday’s law, where Z is the impedance of vacuum.
• Estimate of maximum wave speed cmax (SI unit: m/s) the default is taken from the
speed of light in a vacuum c_const.
FILTER PARAMETERS
The filter provides higher-order smoothing of nodal discontinuous Galerkin
formulations and is intended to be used for absorbing layers, but you can also use it to
stabilize linear wave problems with highly varying coefficients. The filter is constructed
by transforming the solution (in each global time step) to an orthogonal polynomial
representation, multiplying with a damping factor and then transforming back to the
(Lagrange) nodal basis. Select the Activate check box to use this filter.
–1
VΛV
1, 0 ≤ η ≤ η c
Λ mm = σ(η) = η – η c 2s
--------------
– α 1 – η c
-
e , ηc ≤ η ≤ 1
where
im
η = η ( m ) = -------
Np
When using Absorbing Layer features, add an additional Wave Equations feature for the
corresponding domain selection. Select the Activate check box and add filter
parameters. An example of a filter parameter combination that can be used for a Wave
Equations feature active on an Absorbing Layer domain selection is α = 0.1, ηc = 0.01,
and s = 2. However, other combinations could work better, depending on the
particular application.
Absorbing Layers
Reference
1. J.S. Hesthaven and T. Warburton, Nodal Discontinuous Galerkin Methods —
Algorithms, Analysis, and Applications, Springer, 2008.
Initial Values
The Initial Values node adds the initial values for the Electric field and Magnetic field
variables that serve as an initial condition for the transient simulation.
INITIAL VALUES
Enter values or expressions for the initial values of the components of the Electric field
E (SI unit: V/m) and Magnetic field H (SI unit: A/m). The default values are 0 for all
vector components.
∂E
ε 0 ε r ------- – ∇ × H + σE = – J e
∂t
∂H
μ 0 μ r -------- + ∇ × E = – J m
∂t
Electric Field
The Electric Field boundary condition
n × E = n × E0
ELECTRIC FIELD
Enter values or expressions for the components of the Electric field E0 (SI unit: V/m).
n×E = 0
is a special case of the electric field boundary condition that sets the tangential
component of the electric field to zero. It is used for the modeling of a lossless metallic
surface, for example, a ground plane or as a symmetry type boundary condition.
It imposes symmetry for magnetic fields and antisymmetry for electric fields and
electric currents. It supports induced electric surface currents and thus any prescribed
or induced electric currents (volume, surface, or edge currents) flowing into a perfect
electric conductor boundary is automatically balanced by induced surface currents.
Magnetic Field
The Magnetic Field node adds a boundary condition for specifying the tangential
component of the magnetic field at the boundary:
n × H = n × H0
MAGNETIC FIELD
Enter values or expressions for the components of the Magnetic field H0 (SI unit: A/
m).
n×H = 0
is a special case of the surface current density boundary condition that sets the
tangential component of the magnetic field and thus also the surface current density
–n × H = Js
n × ( H1 – H2 ) = Js
specifies a surface current density at both exterior and interior boundaries. The current
density is specified as a three-dimensional vector, but because it needs to flow along
the boundary surface, COMSOL Multiphysics projects it onto the boundary surface
and neglects its normal component. This makes it easier to specify the current density
and avoids unexpected results when a current density with a component normal to the
surface is given.
n × E = Z0 H
IMPEDANCE
Enter the value or expression for the medium Impedance Z0 (SI unit: Ω). By default,
the Z0 uses the value of the vacuum’s impedance. Then select Isotropic, Diagonal,
Symmetric, or Full based on the material characteristics and enter values or expressions
in the field or matrix.
n × E = E0
n × H = H0
specifies the tangential component of both electric and magnetic fields. This boundary
condition is available when Advanced Physics Options is selected in the Show More
Options dialog box on the Model Builder toolbar.
BOUNDARY FLUX/SOURCE
Enter values or expressions for the components of the tangential Electric field E0
(SI unit: V/m) and the tangential Magnetic field H0 (SI unit: A/m).
Background Field
The Background Field feature triggers the scattered field formulation, where the
dependent variable is the relative field. The same wave equations are used as in the full
field formulation, but the total field that enters the equations are written as the sum of
the relative field and the background field, E = Erelative + Ebackground, and it is the
dependent variable Erelative that is solved for. When the background field is a solution
of the wave equation, the relative field is the scattered field.
SETTINGS
Select a Background wave type — User defined (the default), or Modulated Gaussian
pulse.
User Defined
Enter the component expressions for the Background electric field Eb (SI unit: V/m)
and Background magnetic field Hb (SI unit: A/m). The entered expressions must be
differentiable in time domain since the derivative of the background field is used in the
governing equations.
For a modulated Gaussian pulse propagating in the positive x direction, the electric
field is expressed as
x + d offset 2
t – μ – ------------------------ -
vp
- sin 2πf 0 t – ------
1 x
E ( x, t ) = -------------- exp – --------------------------------------------------
τ 2π 2
v p
2τ
where τ is the pulse duration, defined as 1/2f0, μ is a time delay set to 2/f0, and vp is
the phase velocity. The time delay μ is used to excite a modulated Gaussian pulse whose
initial magnitude is very small when it is launched and gradually increases as it
propagates.
Far-Field Domain
To set up a far-field calculation, add a Far-Field Domain node and specify the far-field
domains in its Settings window. Use Far-Field Calculation subnodes (one is added by
default) to specify all other settings needed to define the far-field calculation. By
default, all of the domains are selected. The selection can be modified. In that case,
select only a homogeneous domain or domain group that is outside of all radiating and
scattering objects and which has the material settings of the far-field medium.
Far-Field Calculation
A Far-Field Calculation subnode is added by default to the Far-Field Domain node and
is used to select boundaries corresponding to a single closed surface surrounding all
FAR-FIELD CALCULATION
Enter a Far-field variable name. The default is Efar.
Select as needed the Symmetry in the x=0 plane, Symmetry in the y=0 plane, or
Symmetry in the z=0 plane check boxes to use it your model when calculating the
far-field variable. The symmetry planes have to coincide with one of the Cartesian
coordinate planes.
When a check box is selected, also choose the type of symmetry to use from the
Symmetry type list that appears — Symmetry in E (PMC) or Symmetry in H (PEC). The
selection should match the boundary condition used for the symmetry boundary.
Using these settings, include the parts of the geometry that are not in the model for
symmetry reasons in the far-field analysis.
From the Boundary relative to domain list, select Inside or Outside (the default) to define
if the selected boundaries are defining the inside or outside of the far-field domain (that
is, whether facing away from infinity or toward infinity).
A Time to Frequency FFT study step must be added after the Time Dependent
study step to generate the necessary frequency-domain data, used in the
far-field analysis.
Use this physics interface in 2D and 3D when approximating the scattered far-field of
an object configured only by a perfect electric conductor boundary condition. The
physics interface supports the Frequency Domain study type. The Frequency Domain
study can also be used for sweeping the background field by a sequence of frequencies.
When this physics interface is added, these default nodes are also added to the Model
Builder: Asymptotic Scattering, Far-Field Calculation, and Initial Values. No additional
boundary feature is needed in general. However, from the Physics toolbar, a new
Far-Field Calculation can be added to process the far-field calculation on user-defined
boundary selections. You can also right-click Electromagnetic Waves, Asymptotic
Scattering to select physics features from the context menu.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is ewas.
Physics-Controlled Mesh
The physics-controlled mesh is controlled from the Mesh node’s Settings window (if the
Sequence type is Physics-controlled mesh). In the table in the Physics-Controlled Mesh
section, find the physics interface in the Contributor column and select or clear the
check box in the Use column on the same row for enabling (the default) or disabling
contributions from the physics interface to the physics-controlled mesh.
When the Use check box is selected for the physics interface in the section for the
physics interface below the table, choose one of the four options for the Maximum mesh
element size control parameter — From study (the default), User defined, Frequency, or
Wavelength. When From study is selected, 1/5 of the vacuum wavelength from the
highest frequency defined in study step is used for the maximum mesh element size.
For the option User defined, enter a suitable Maximum element size in free space. For
example, 1/5 of the vacuum wavelength or smaller. When Frequency is selected, enter
the highest frequency intended to be used during the simulation. The maximum mesh
element size in free space is 1/5 of the vacuum wavelength for the entered frequency.
For the Wavelength option, enter the smallest vacuum wavelength intended to be used
during the simulation. The maximum mesh element size in free space is 1/5 of the
entered wavelength.
FORMULATION
For Scattered field select a Background wave type according to the following table:
TABLE 4-2: BACKGROUND WAVE TYPE BASED ON COMPONENT DIMENSION
User Defined
Enter the component expressions for the Background electric field Eb (SI unit: V/m).
The entered expressions must be differentiable.
• Enter an Electric field amplitude E0 (SI unit: V/m). The default is 1 V/m.
Figure 4-5: Schematic of the directions for the wave vector k, the electric field E0, and the
roll, pitch and yaw rotations. The top image represents an initial wave propagating in the
x direction with a polarization along the z direction.
DEPENDENT VARIABLES
The dependent variables (field variables) are for the Electric field E and its components
(in the Electric field components fields). The name can be changed but the names of
fields and dependent variables must be unique within a model.
DISCRETIZATION
Select the shape order for the Electric field dependent variable — Linear, Quadratic (the
default), or Cubic. For more information about the Discretization section, see Settings
for the Discretization Sections in the COMSOL Multiphysics Reference Manual.
Far-Field Calculation
Use a Far-Field Calculation to select boundaries corresponding to a single closed surface
surrounding all radiating and scattering objects. By default, all exterior boundaries of
the simulation domain are selected. Also use this feature to indicate symmetry planes
and symmetry cuts applied to the geometry.
FAR-FIELD CALCULATION
Enter a Far-field variable name. The default is Efar.
Select, as applicable for your model, the Symmetry in the x=0 plane, Symmetry in the
y=0 plane, or Symmetry in the z=0 plane check boxes to simplify the calculation of the
far-field variable. The symmetry planes must coincide with one of the Cartesian
coordinate planes.
When a check box is selected, also choose the type of symmetry to use from the
Symmetry type list that appears: Symmetry in E (PMC) or Symmetry in H (PEC). The
selection should match the boundary condition used for the symmetry boundary.
Using these settings, includes the parts of the geometry that are not in the model for
symmetry reasons in the far-field analysis.
From the Boundary relative to domain list, select Inside or Outside (the default) to define
if the selected boundaries are defining the inside or outside of the far-field domain (that
is, whether facing away from infinity or toward infinity).
Initial Values
The Initial Values node adds an initial value for the electric field that can be used to
perform the far-field transformation for a single frequency case without running the
complete computation.
INITIAL VALUES
Enter values or expressions for the initial values of the components of the Electric field
E (SI unit: V/m). The default values are set by the background field.
∂D
∇ × H = J + -------
∂t
∂B
∇ × E = – -------
∂t
∂εE
∇ × H = σE + ----------
∂t
∂-------
H-
∇ × E = –μ
∂t
jωt
E ( x, y, z, t ) = E ( x, y, z )e
jωt
H ( x, y, z, t ) = H ( x, y, z )e
the two laws can be combined into a time harmonic equation for the electric field, or
a similar equation for the magnetic field
∇ × ( μ – 1 ∇ × E ) – ω 2 εE = 0
–1
∇ × ( ε ∇ × H ) – ω 2 μH = 0
The first of these, based on the electric field is used in The Electromagnetic Waves,
Frequency Domain Interface.
Using the relation εr = n2, where n is the refractive index, the equation can
alternatively be written
2 2
∇ × ( ∇ × E ) – k0 n E = 0
ω
k 0 = ω ε 0 μ 0 = -----
c0
When the equation is written using the refractive index, the assumption is that μr = 1
and σ = 0 and only the constitutive relations for linear materials are available. When
solving for the scattered field the same equations are used but E = Esc + Ei and Esc is
the dependent variable.
EIGENFREQUENCY ANALYSIS
When solving the frequency domain equation as an eigenfrequency problem the
eigenvalue is the complex eigenfrequency λ = jω + δ, where δ is the damping of the
solution. The Q-factor is given from the eigenvalue by the formula
ω
Q fact = ---------
2δ
The spatial parameter, α = δz + jβ = −λ, can have a real part and an imaginary part. The
propagation constant is equal to the imaginary part, and the real part, δz, represents
the damping along the propagation direction. When solving for all three electric field
components the allowed anisotropy of the optionally complex relative permittivity and
relative permeability is limited to:
In 2D, the electric field varies with the out-of-plane wave number kz as
˜
E ( x, y, z ) = E ( x, y ) exp ( – ik z z ) .
˜ 2 ˜
( ∇ – ik z z ) × [ μ r–1 ( ∇ – ik z z ) × E ] – k 0 ε rc E = 0 ,
Similarly, in 2D axisymmetry, the electric field varies with the azimuthal mode number
m as
˜
E ( r, ϕ, z ) = E ( r, z ) exp ( – imϕ )
∇ – im ˜ ˜
----- ϕ × μ r–1 ∇ – i m
----- ϕ × E – k 20 ε rc E = 0 ,
r r
In-plane TM Waves
The TM waves polarization has only one magnetic field component in the z direction,
and the electric field lies in the modeling plane. Thus the time-harmonic fields can be
obtained by solving for the in-plane electric field components only. The equation is
formally the same as in 3D, the only difference being that the out-of-plane electric field
component is zero everywhere and that out-of-plane spatial derivatives are evaluated
for the prescribed out-of-plane wave vector dependence of the electric field.
2
– ∇ ⋅ ( μ̃ r ∇E z ) – ε rzz k 0 E z = 0
where
T
μr
μ̃ r = -------------------
det ( μ r )
To be able to write the fields in this form, it is also required that εr, σ, and μr are
nondiagonal only in the xy-plane. μr denotes a 2-by-2 tensor, and εrzz and σzz are the
relative permittivity and conductivity in the z direction.
Axisymmetric TM Waves
A TM wave has a magnetic field with only a ϕ component and thus an electric field
with components in the rz-plane only. The equation is formally the same as in 3D, the
only difference being that the ϕ component is zero everywhere and that spatial
derivatives with respect to ϕ are evaluated for the prescribed azimuthal mode number
dependence of the electric field.
Axisymmetric TE Waves
A TE wave has only an electric field component in the ϕ direction, and the magnetic
field lies in the modeling plane. Given these constraints, the 3D equation can be
simplified to a scalar equation for E ϕ . To write the fields in this form, it is also required
that εr and μr are nondiagonal only in the rz-plane. μr denotes a 2-by-2 tensor, while
ε rϕϕ and σ ϕϕ are the relative permittivity and conductivity in the ϕ direction.
Electric Losses
The frequency domain equations allow for several ways of introducing electric losses.
Finite conductivity results in a complex permittivity,
σ
ε c = ε – j ----
ω
ε c = ε 0 ( ε' – jε'' )
where ε' is the real part of εr, and all losses (dielectric and conduction losses) are given
by ε''. The dielectric loss model can also single out the losses from finite conductivity
(so that ε'' only represents dielectric losses) resulting in:
σ
ε c = ε 0 ε' – j --------- + ε''
ωε 0
ε c = ε 0 ε' ( 1 – j tan δ )
In optics and photonics applications, the refractive index is often used instead of the
permittivity. In materials where μr is 1, the relation between the complex refractive
index
n = n – jκ
2
ε rc = n
that is
2 1 2 2
n = --- ( ε' r + ε' r + ε'' r )
2
2 1 2 2
κ = --- ( – ε' r + ε' r + ε'' r )
2
In the physics and optics literature, the time harmonic form is often written with a
minus sign (and “i” instead of “j”):
– i ωt
E ( x, y, z, t ) = E ( x, y, z )e
Magnetic Losses
The frequency domain equations allow for magnetic losses to be introduced as a
complex relative permeability.
μ r = ( μ' – jμ'' )
The complex relative permeability can be combined with any electric loss model except
refractive index.
∂A ∂ ∂A –1
μ0 σ + μ0 ε + ∇ × ( μr ∇ × A ) = 0
∂t ∂t ∂t
Using the relation εr = n2, where n is the refractive index, the equations can
alternatively be written
∂ n 2 ∂A + ∇ × ( ∇ × A ) = 0
μ0 ε0
∂t ∂t
WAVES IN 2D
In 2D, different polarizations can be chosen by selecting to solve for a subset of the
3D vector components. When selecting all three components, the 3D equation applies
with the addition that out-of-plane spatial derivatives are set to zero.
In-plane TM Waves
The TM waves polarization has only one magnetic field component in the z direction,
and thus the electric field and vector potential lie in the modeling plane. Hence it is
obtained by solving only for the in-plane vector potential components. The equation
is formally the same as in 3D, the only difference being that the out-of-plane vector
potential component is zero everywhere and that out-of-plane spatial derivatives are
set to zero.
In-plane TE Waves
As the field propagates in the modeling xy-plane a TE wave has only one nonzero
vector potential component, namely in the z direction. The magnetic field lies in the
modeling plane. Thus the equation in the time domain can be simplified to a scalar
equation for Az:
∂A z ∂ ∂A z
+ μ 0 ε 0 ε r
–1
μ0 σ + ∇ ⋅ ( μ r ( ∇A z ) ) = 0
∂t ∂t ∂t
Using the relation εr = n2, where n is the refractive index, the equation can
alternatively be written
2 ∂A z
μ 0 ε 0 ∂ n + ∇ ⋅ ( ∇A z ) = 0
∂t ∂t
Axisymmetric TM Waves
TM waves have a magnetic field with only a ϕ component and thus an electric field
and a magnetic vector potential with components in the rz-plane only. The equation
is formally the same as in 3D, the only difference being that the ϕ component is zero
everywhere and that spatial derivatives with respect to ϕ are set to zero.
Axisymmetric TE Waves
A TE wave has only a vector potential component in the ϕ direction, and the magnetic
field lies in the modeling plane. Given these constraints, the 3D equation can be
simplified to a scalar equation for A ϕ . To write the fields in this form, it is also required
that εr and μr are nondiagonal only in the rz-plane. μr denotes a 2-by-2 tensor, while
ε rϕϕ and σ ϕϕ are the relative permittivity and conductivity in the ϕ direction.
Vector Elements
Whenever solving for more than a single vector component, it is not possible to use
Lagrange elements for electromagnetic wave modeling. The reason is that they force
the fields to be continuous everywhere. This implies that the physics interface
conditions, which specify that the normal components of the electric and magnetic
fields are discontinuous across interior boundaries between media with different
permittivity and permeability, cannot be fulfilled. To overcome this problem, the
Electromagnetic Waves, Frequency Domain interface uses vector elements, which do
not have this limitation.
The solution obtained when using vector elements also better fulfills the divergence
conditions ∇ · D = 0 and ∇ · B = 0 than when using Lagrange elements.
Using the default parameters for the eigenfrequency study, it might find a large
number of false eigenfrequencies, which are almost zero. This is a known consequence
of using vector elements. To avoid these eigenfrequencies, change the parameters for
the eigenvalue solver in the Study Settings. Adjust the settings so that the solver
searches for eigenfrequencies closer to the lowest eigenfrequency than to zero.
In the paraxial approximation the field for a Gaussian beam propagating along the
z-axis is defined below,
w0 2 2
ρ - ρ
E G ( x, y, z ) = E G0 ------------ exp – -------------- – jkz – jk ---------------- + j η ( z ) ,
w(z) w (z)
2 2R ( z )
where w0 is the beam radius, p0 is the focal plane on the z-axis, EG0 is the Gaussian
beam electric field amplitude and the spot radius for different positions along the
propagation axis is given by
z – p0 2
w ( z ) = w 0 1 + --------------- .
z0
defines the radius of curvature for the phase of the field and the so-called Gouy phase
shift is given by
z–p
η ( z ) = atan --------------0- .
z0
The equations above are expressed using the Rayleigh range z0 and the transverse
coordinate ρ, defined by
2
k0 w0 2 2 2
z 0 = --------------, ρ = x + y .
2
Note that the time-harmonic ansatz in COMSOL is ejωt, and with this convention, the
beam above propagates in the +z direction. The equations are modified accordingly for
beams propagating along the other coordinate axes.
The field for a Gaussian beam is defined in a similar way for 2D components. In the
particular case where the beam propagates along the x-axis, the field is defined as
2 2
w0 y - η(x)
y - -----------
E G ( x, y ) = E G0 ------------ exp – -------------- – jkx – jk --------------- +j .
w(x) w (x)
2 2R ( x ) 2
For a beam propagating along the y-axis, the coordinates x and y are interchanged.
Notice that the expressions above for Gaussian beams are not solutions to the
Helmholtz equation, but to the so called paraxial approximation of the Helmholtz
equation. This means that these equations become less accurate the smaller the spot
radius is and should not be used when the spot radius is of the same size as or smaller
than the wavelength.
To circumvent the problem that the paraxial approximation formula is not a solution
to the Helmholtz equation, a plane wave expansion can be used to approximate a
Gaussian beam background field. Since each plane wave is a solution to Helmholtz
equation, also the expansion is a solution to Helmholtz equation.
The plane wave expansion approximates the Gaussian distribution in the focal plane
where the beam is assumed to be propagating in the z direction, the focal plane is
spanned by the x and y coordinates, e is the unit magnitude transverse polarization in
the focal plane, l and m denote the indices for the wave vectors, the index n accounts
for the two polarizations per wave vector klm, almn is the amplitude, un(klm) is the
unit magnitude polarization, and r is the position vector.
Multiplying with the conjugate of the exponential factor above and the polarization
factor un(klm) and applying a surface integral over the entire focal plane allows us to
extract the amplitudes as
2 2 2
E 0 w 0 ( e ⋅ u n ( k lm ) ) k t, lm w 0
a lmn = -------------------------------------------------- exp – --------------------- ,
4π 4
There are several possible approaches to compute an average material property starting
from the material properties and the volume fraction of each material.
The following sections illustrate the different formulas available to compute the
effective electrical conductivity, the effective relative permittivity and the effective
relative permeability of a homogenized medium. In the following, volume fractions
of the materials are indicated with θi, where i is the material index, and they are
assumed to be fractional (between 0 and 1). Up to five different materials can be
specified as phases of the mixture. Typically, their volume fractions should add up to 1.
σ = θi σi = θ1 σ1 + θ2 σ 2 + …
i=1
If the conductivities are defined by second order tensors (such as for anisotropic
materials), the volume average is applied element by element.
n
θi θ1 θ2
--1- =
σ ----σ-i = ------ + ------ + …
σ1 σ2
i=0
If the conductivities are defined by second order tensors, the inverse of the tensors are
used.
POWER LAW
A power law gives the following expression for the equivalent conductivity:
n
θi θ θ
σ = ∏ σi = σ 11 σ 22 …
i=0
ε = θi εi = θ1 ε1 + θ2 ε2 + …
i=1
n
1 θi θ1 θ2
--- =
ε ----εi = ----- + ----- + …
ε1 ε2
i=0
If the permittivity is defined by a second-order tensor, the inverse of the tensor is used.
POWER LAW
A power law gives the following expression for the equivalent permittivity:
n
θi θ θ
ε = ∏ εi = ε 11 ε 22 …
i=0
μ = θi μi = θ1 μ1 + θ2 μ2 + …
i=1
n
θi θ1 θ2
--1- =
μ ----μ-i = ------ + ------ + …
μ1 μ2
i=0
If the permeability is defined by a second-order tensor, the inverse of the tensor is used.
POWER LAW
A power law gives the following expression for the equivalent permeability:
n
θi θ θ
μ = ∏ μi = μ 11 μ 22 …
i=0
The porous rocks and sediments can have even lower conductivities. In variably
saturated soils, the conductivity of air is roughly ten orders of magnitude lower than
the ground water. A simple volume average (of either conductivity or resistivity) in
rocks or soils might give different results compared to experimental data.
Since most crustal rocks, sedimentary rocks, and soils are formed by nonconducting
materials, Archie (Ref. 1) assumed that electric current are mainly caused by ion fluxes
through the pore network. Originally, Archie’s law is an empirical law for the effective
conductivity of a fully saturated rock or soil, but it can be extended to variably
saturated porous media.
Archie’s law relates the effective conductivity to the fluid conductivity σL, fluid
saturation sL, and porosity εp:
n m
σ = sL εp σL
here, m is the cementation exponent, a parameter that describes the connectivity of the
pores. The cementation exponent normally varies between 1.3 and 2.5 for most
sedimentary rocks and is close to 2 for sandstones. The lower limit m = 1 represents a
volume average of the conductivities of a fully saturated, insulating (zero conductivity)
porous matrix, and a conducting fluid. The saturation coefficient n is normally close
to 2. The ratio F = σL/σ is called the formation factor.
Archie’s law does not take care of the relative permittivity of either fluids or solids, so
the effective relative permittivity of the porous medium is normally consider as εr = 1.
The distribution of the electric potential V and the current I describes the propagation
of the signal wave along the line. The following equations relate the current and the
electric potential
∂V
= – ( R + jωL )I (4-1)
∂x
∂I
= – ( G + jωC )V (4-2)
∂x
Equation 4-1 and Equation 4-2 can be combined to the second-order partial
differential equation
where
γ = ( R + jωL ) ( G + jωC ) = α + jβ
Here γ, α, and β are called the complex propagation constant, the attenuation
constant, and the (real) propagation constant, respectively.
– γx γx (4-4)
V ( x ) = V+ e + V- e
γ – γx γx
I ( x ) = --------------------- ( V + e – V - e )
R + jωL
V R + jωL R + jωL-
Z 0 = ---- = --------------------- = ---------------------
I γ G + jωC
To make sure that the current is conserved across interior boundaries, COMSOL
Multiphysics solves the following wave equation (instead of Equation 4-3)
1 - ∂V
∂ -------------------- – ( G + jωC )V = 0 (4-5)
∂ x R + jωL ∂ x
V1 = V2 (4-6)
I1 = I2 (4-7)
In Equation 4-6 and Equation 4-7, the indices 1 and 2 denote the domains on the two
sides of the boundary. The currents flowing out of a boundary are given by
n i ⋅ ∇V i
I i = – ------------------------- , i = 1, 2
R i + jωL i
Because V is solved for, the electric potential is always continuous, and thus
Equation 4-6 is automatically fulfilled. Equation 4-7 is equivalent to the natural
boundary condition
1 - ∂V – --------------------------- ∂V
-------------------------- 1
= 0
R 2 + jωL 2 ∂ x 2 R 1 jωL 1 ∂ x
+ 1
When the transmission line is terminated by a load impedance, as Figure 4-6 shows,
the current through the load impedance is given by
V(L )
I ( L ) = ------------- (4-8)
ZL
1 - ∂V V
-------------------- + ------- = 0 (4-9)
R + jωL ∂ x Z L
The Open Circuit boundary condition is obtained by letting the load impedance
become infinitely large, that is, no current flows through the load impedance.
On the other hand, the Short Circuit boundary condition specifies that the voltage at
the load is zero. In COMSOL Multiphysics this is implemented as a constraint on the
electric potential.
V ( 0 ) = V = V0 + V-
1 ∂V 1 2V 0 – V
I ( 0 ) = – --------------------- = ------ ( V 0 – V - ) = --------------------
R + jωL ∂ x x=0
Z0 Z0
1 ∂V V – 2V 0
– --------------------- + -------------------- = 0
R + jωL ∂ x Z0
For the Lumped Port boundary condition, the port current (positive when entering
the transmission line) defines the boundary condition as
1 ∂V
– --------------------- – I port = 0
R + jωL ∂ x
2V 0 – V
I port = --------------------
Z0
for a Cable lumped port (see the Lumped Port section for a description of the lumped
port settings).
• The Equations
• In-Plane E Field or In-Plane H Field
• Fluxes as Dirichlet Boundary Conditions
• Absorbing Layers
The Equations
Maxwell’s equations are a set of equations, written in differential or integral form,
stating the relationships between the fundamental electromagnetic quantities. These
quantities are the:
For general time-varying fields, the differential form of Maxwell’s equations can be
written as
∂D
∇ × H = J + -------
∂t
∂B
∇ × E = – ------- (4-10)
∂t
∇⋅D = ρ
∇⋅B = 0
The first two equations are also called Maxwell-Ampère’s law and Faraday’s law,
respectively. Equation three and four are two forms of Gauss’ law, the electric and
magnetic form, respectively.
D = ε0 E + P
B = μ0 ( H + M ) (4-11)
J = σE
Here ε0 is the permittivity of a vacuum, μ0 is the permeability of a vacuum, and σ the
electric conductivity of the medium. In the SI system, the permeability of a vacuum is
chosen to be 4π·10−7 H/m. The velocity of an electromagnetic wave in a vacuum is
given as c0 and the permittivity of a vacuum is derived from the relation
1 – 12 1 –9
ε 0 = ----------
2
= 8.854 ⋅ 10 F/m ≈ --------- ⋅ 10 F/m
c0 μ0 36π
The electric polarization vector P describes how the material is polarized when an
electric field E is present. It can be interpreted as the volume density of electric dipole
moments. P is generally a function of E. Some materials might have a nonzero P also
when there is no electric field present.
The magnetization vector M similarly describes how the material is magnetized when
a magnetic field H is present. It can be interpreted as the volume density of magnetic
dipole moments. M is generally a function of H. Permanent magnets, for example,
have a nonzero M also when there is no magnetic field present.
To get a wave equation for the E field, for example, take the curl of the second
equation in Equation 4-10 (previously divided by μ0), and insert it into the time
derivative of the first row in Equation 4-10
2 2
∂M ∂E ∂ E- ---------
∂ P
– ∇ × ----- ∇ × E + -------- = σ ------- + ε 0 ---------
1
+ 2-
μ0 ∂t ∂t ∂t 2 ∂t
this is referred as curl-curl formulation in the literature (second order time derivatives
and second order space derivatives).
LINEAR MATERIALS
In the simplest case linear materials, the polarization is directly proportional to the
electric field, that is
∂P ⁄ ∂E = ε 0 χ e and P = ε 0 χ e E
∂M ⁄ ∂H = χ m and M = χ m H
As a consequence, for linear materials, the constitutive relations in Equation 4-11 can
be written as
D = ε 0 E + P = ε 0 ( 1 + χ e )E = ε 0 ε r E
B = μ 0 ( H + M ) = μ 0 ( 1 + χ m )H = μ 0 μ r H
Here, ε = ε0εr and μ = μ0μr are the permittivity and permeability of the material. The
relative permittivity εr and the relative permeability μr are usually scalar properties but
these can be second-rank symmetric (Hermitian) tensors for a general anisotropic
material.
∂E
∇ × H = σE + ε 0 ε r -------
∂t
(4-12)
∂H
∇ × E = – μ 0 μ r --------
∂t
∂u
da + ∇ ⋅ Γ(u) = f
∂t
Maxwell’s equations in 3D
∂E
ε 0 ε r ------- – ∇ × H = – σE
∂t
∂H
μ 0 μ r -------- + ∇ × E = 0
∂t
∂E
d E ------- + ∇ ⋅ Γ E ( H ) = f
∂t
∂-------
H
dH - + ∇ ⋅ ΓH ( E ) = 0
∂t
d E = ε 0 ε r and d H = μ 0 μ r
0 h3 –h2 0 e3 –e2
ΓE ( H ) = – –h3 0 h1 and Γ H ( E ) = – e 3 0 e1
h2 –h1 0 e2 –e1 0
When using SI units (or other) for the electromagnetic fields and material
properties, the Lax-Friedrichs flux parameter is not dimensionless and
must have units of τE = 1/(2Z) for Ampère’s law and τH = Z/2 for
Faraday’s law, where Z is the impedance of the medium.
TM WAVES IN 2D
For TM waves in 2D, solve for an in-plane electric field vector and one out-of-plane
variable for the magnetic field. Maxwell’s equations then read
∂E
ε 0 ε r ------- + ∇ ⋅ Γ E ( H ) = – σ ⋅ E
∂t
(4-14)
∂H
μ 0 μ r -------- + ∇ ⋅ Γ H ( E ) = 0
∂t
0 –h3
ΓE ( H ) = and Γ H ( E ) = e 2 – e 1 (4-15)
h3 0
The default Lax-Friedrichs flux parameters are τE = 1/(2Z) for Ampère’s law, and the
scalar τH = Z/2 for Faraday’s law, where Z is the impedance of a vacuum.
TE WAVES IN 2D
For TE waves in 2D, solve for an in-plane magnetic field vector and one out-of-plane
variable for the electric field. Maxwell’s equations then read
∂E
ε 0 ε r ------- + ∇ ⋅ Γ E ( H ) = – σE
∂t
(4-16)
∂H
μ 0 μ r -------- + ∇ ⋅ Γ H ( E ) = 0
∂t
0 e3
Γ E ( H ) = – h 2 h 1 and Γ H ( E ) = (4-17)
–e3 0
The default Lax-Friedrichs flux parameters are τE = 1/(2Z) for Ampère’s law, and two
scalar τH = Z/2 for Faraday’s law, where Z is the impedance of a vacuum.
∂E
ε 0 ε r ------- + ∇ ⋅ Γ E ( H ) = – σE
∂t
∂H
μ 0 μ r -------- + ∇ ⋅ Γ H ( E ) = 0
∂t
0 –h3 h2 0 e3 –e2
ΓE ( H ) = h3 0 – h 1 and Γ H ( E ) = – e 3 0 e1
–h2 h1 0 e2 –e1 0
For Ampère’s law, the normal to the flux term on exterior boundaries reads
n ⋅ ΓE ( H ) = –n × H
n ⋅ ΓH ( E ) = n × E
which means that normal fluxes on external boundaries can only prescribe tangential
components for the fields.
BOUNDARY CONDITIONS
The boundary conditions for outer boundaries are computed from the normal fluxes
n · ΓH(E) and n · ΓE(H).
Absorbing Layers
The Electromagnetic Waves, Time Explicit Interface includes so-called absorbing
layers, also often referred to as sponge layers. The layers work by combining three
techniques: a scaling system, filtering, and simple nonreflecting conditions. For a
review of the method see, for example, Ref. 1.
The layers are set up by adding the Absorbing Layer under the Definitions node. This
adds a special scaled system. The scaling effectively slows down the propagating waves
and ensures that they hit the outer boundary in the normal direction. For the
Absorbing Layer domain selection, add an additional Wave Equations feature, mark
the Activate check box under the Filter Parameters section, and enter filter parameters.
Filtering attenuates and filters out high-frequency components of the wave. Finally, at
the outer boundary of the layer add a simple Scattering Boundary Condition
condition, which will work well to remove all remaining waves as normal incidence has
been ensured.
For more detailed information about the filter see the Filter Parameters
section under Wave Form PDE in the COMSOL Multiphysics Reference
Manual.
For more detailed on the Geometry and Scaling see the Infinite Elements,
Perfectly Matched Layers, and Absorbing Layers in the COMSOL
Multiphysics Reference Manual.
For the layers to work optimally the filter should not be too aggressive. Moreover, the
scaled coordinates in the layer domain should also vary smoothly. To inspect the scaled
system you can plot the coordinate variables x_absorb_ab1, y_absorb_ab1, and
z_absorb_ab1. Using the absorbing layers with the three combined techniques will
enable the reduction of spurious reflections by a factor between 100 and 1000
compared to the incident amplitude.
Reference
1. P.G. Petropoulos, L. Zhao, and A.C. Cangellaris, “A Reflectionless Sponge Layer
Absorbing Boundary Condition for the Solution of Maxwell’s Equations with
High-Order Staggered Finite Difference Schemes,” J. Comp. Phys., vol. 139, pp. 184–
208, 1998.
AC/DC Interfaces
This chapter summarizes the functionality of the electrical circuit interface found
under the AC/DC branch ( ) when adding a physics interface.
In this chapter:
209
The Electrical Circuit Interface
The Electrical Circuit (cir) interface ( ), found under the AC/DC branch ( ) when
adding a physics interface, is used to model currents and voltages in circuits including
voltage and current sources, resistors, capacitors, inductors, and semiconductor
devices. Models created with the Electrical Circuit interface can include connections
to distributed field models. The physics interface supports stationary,
frequency-domain and time-domain modeling and solves Kirchhoff’s conservation
laws for the voltages, currents and charges associated with the circuit elements.
When this physics interface is added, it adds a default Ground Node feature and
associates that with node zero in the electrical circuit.
Circuit nodes are nodes in the electrical circuit (electrical nodes) and
should not be confused with nodes in the Model Builder tree of the
COMSOL Multiphysics software. Circuit node names are not restricted
to numerical values but can contain alphanumeric characters.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is cir.
GROUND CONNECTION
Set the Node name for the ground node in the circuit. The convention is to use 0 (zero)
for the ground node. If adding more ground nodes, each must have a unique node
name (number).
Voltmeter
The Voltmeter ( ) feature connects a voltmeter (voltage measurement device)
between two nodes in the electrical circuit. A voltmeter behaves electrically as an open
circuit. The voltmeter node adds a Probe sampling the voltage across it.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
Ampère Meter
The Ammeter ( ) feature connects an ammeter (current measurement device)
between two nodes in the electrical circuit. An ammeter behaves electrically as a short
circuit. The ammeter node adds a Probe sampling the current through it.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
DEVICE PARAMETERS
Enter the Resistance of the resistor.
Resistor
The Resistor ( ) feature connects a resistor between two nodes in the electrical
circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
Capacitor
The Capacitor ( ) feature connects a capacitor between two nodes in the electrical
circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the capacitor.
DEVICE PARAMETERS
Enter the Capacitance of the capacitor.
Inductor
The Inductor ( ) feature connects an inductor between two nodes in the electrical
circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the inductor.
DEVICE PARAMETERS
Enter the Inductance of the inductor.
Voltage Source
The Voltage Source ( ) feature connects a voltage source between two nodes in the
electrical circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the voltage source. The first node
represents the positive reference terminal.
• For a DC-source, the Voltage Vsrc (default value: 1 V). DC-sources are active in
Stationary and Time-Dependent studies.
• For an AC-source: the Voltage Vsrc (default value: 1 V) and the Phase Θ (default
value: 0 rad). AC-sources are active in Frequency Domain studies only.
• For a sine source: the Voltage Vsrc (default value: 1 V), the Offset Voff (default value:
0 V), the Frequency (default value: 1 kHz), and the Phase Θ (default value: 0 rad).
The sine sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
• For a pulse source: the Voltage Vsrc (default value: 1 V), the Offset Voff (default value:
0 V), the Delay td (default value: 0s), the Rise time tr and Fall time tf (default values:
0 s), the Pulse width pw (default value: 1 μs), and the Period Tper (default value: 2 μs).
The pulse sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
For the AC source, the frequency is a global input set by the solver. AC
sources should be used in Frequency-domain studies only. Do not use the
Sine source unless the model is time-dependent.
Current Source
The Current Source ( ) feature connects a current source between two nodes in the
electrical circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the current source. The first node
represents the positive reference terminal from where the current flows through the
source to the second node.
• For a DC-source, the Current isrc (default value: 1 A). DC-sources are active in
Stationary and Time-Dependent studies.
• For an AC-source: the Current isrc (default value: 1 A) and the Phase Θ (default
value: 0 rad). AC-sources are active in Frequency Domain studies only.
• For a sine source: the Current isrc (default value: 1 A), the Offset ioff (default value:
0 A), the Frequency (default value: 1 kHz), and the Phase Θ (default value: 0 rad).
The sine sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
• For a pulse source: the Current isrc (default value: 1 A), the Offset ioff (default value:
0 A), the Delay td (default value: 0 s), the Rise time tr and Fall time tf (default values:
0 s), the Pulse width pw (default value: 1 μs), and the Period Tper (default value: 2 μs).
The pulse sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
For the AC source, the frequency is a global input set by the solver. AC
sources should be used in frequency-domain studies only. Do not use the
Sine source unless the model is time-dependent.
NODE CONNECTIONS
Specify four Node names: the first pair for the connection nodes for the voltage source
and the second pair defining the input control voltage. The first node in a pair
represents the positive reference terminal.
NODE CONNECTIONS
Specify four Node names: the first pair for the connection nodes for the current source
and the second pair defining the input control voltage. The first node in a pair
represents the positive voltage reference terminal or the one from where the current
flows through the source to the second node.
DEVICE PARAMETERS
There are two options to define the relationship between the control voltage and
resulting current. The Use gain method defines the resulting current to be the control
voltage multiplied by the gain (SI units: S). The Custom expression method can define
the relationship with an arbitrary expression.
NODE CONNECTIONS
Set two Node names for the connection nodes for the voltage source. The first node in
a pair represents the positive reference terminal.
DEVICE PARAMETERS
There are two options to define the relationship between the control current and
resulting voltage. The Use gain method defines the resulting voltage to be the control
current multiplied by the gain (SI units: Ω). The Custom expression method can define
the relationship with an arbitrary expression.
NODE CONNECTIONS
Specify two Node names for the connection nodes for the current source. The first node
in a pair represents the positive reference terminal from where the current flows
through the source to the second node.
DEVICE PARAMETERS
There are two options to define the relationship between the control current and
resulting current. The Use gain method defines the resulting current to be the control
current multiplied by the gain. The Custom expression method can define the
relationship with an arbitrary expression.
Switch
The Switch ( ) feature is used to connect or disconnect the conducting path in a
circuit under specific conditions.
NODE CONNECTIONS
Specify two Node names for the connection nodes for the current source. The first node
in a pair represents the positive reference terminal from where the current flows
through the source to the second node.
SWITHCH CONDITIONS
There are three types of conditions, Voltage controlled, Current controlled, and Custom
expressions. For each type of condition there are two conditions, one for turn on and
one for turn off. The on condition is true if the On condition expression is larger than
zero, while the off condition is true if the Off condition is less than zero.
The Initial state list has three options, Use on condition, Use off condition, and Boolean
expression. The two former options mean that the switch will have an initial state
matching to the on or off condition. The third option makes the switch's initial state
match a custom Boolean expression. Separating on, off, and initial states makes the
switch more flexible and can support Schmitt-trigger style switches and various latches.
SWITCH PARAMETERS
When the switch is in the on state it has a non-zero resistance specified by the On
resistance expression. For the off state no current flows through the switch (infinite
resistance). There is also a transition time for the switch to turn on and off set by the
Switching time expression. The switch triggers an implicit event that updates a discrete
state variable (with suffix _state).
Subcircuit Definition
The Subcircuit Definition ( ) feature is used to define subcircuits, which can be
inserted as devices into the main circuit using Subcircuit Instance nodes. Create the
subcircuit by adding subnodes to the Subcircuit Definition node, either by using the
Physics toolbar, or by right-clicking the Subcircuit Definition.
SUBCIRCUIT PINS
Define the Pin names at which the subcircuit connects to the main circuit or to other
subcircuits when referenced by a Subcircuit Instance node. The Pin names refer to
circuit nodes in the subcircuit. The order in which the Pin names are defined is the
order in which they are referenced by a Subcircuit Instance node. The devices
constituting the subcircuit should be connected only to the subcircuit’s pins and to
themselves.
INPUT PARAMETERS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options. Specify input parameters to a subcircuit that can be changed from a
subcircuit instance. These input parameters can be used in all expression-style edit
fields that affect the parameters of a device, for example, resistance, capacitance, and
current gain. In this way, a subcircuit can represent a parameterized custom device
model.
NODE CONNECTIONS
Select the Name of subcircuit link from the list of defined subcircuits in the circuit model
and the circuit Node names at which the subcircuit instance connects to the main circuit
or to another subcircuit if used therein.
INPUT PARAMETERS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options. Specify input parameters to a subcircuit that can be changed from a
subcircuit instance. These input parameters can be used in all expression-style edit
fields that affect the parameters of a device, for example, resistance, capacitance, and
current gain. In this way, a subcircuit can represent a parameterized custom device
model.
NODE CONNECTIONS
Specify three Node names for the connection nodes for the BJT device. These represent
the collector, base, and emitter nodes for the NPN transistor, and the emitter, base, and
collector nodes for the PNP transistor.
NODE CONNECTIONS
Specify four Node names for the connection nodes for the n-Channel MOSFET or
p-Channel MOSFET device. These represent the drain, gate, source, and bulk nodes,
respectively.
MODEL PARAMETERS
Specify the Model Parameters. Reasonable defaults are provided but for any particular
MOSFET, the device manufacturer should be the primary source of information.
Mutual Inductance
The Mutual Inductance allows specifying a coupling between two existing Inductor
features in the circuit. The mutual inductance of the coupling is
where k is the coupling factor and L1 and L2 are the inductances of the inductors.
DEVICE PARAMETERS
Enter values or expressions for the:
• Coupling factor k (dimensionless). The value must be between 0 and 1, and the
default is 0.98.
• First inductance L1 (SI unit: H) and Second inductance L2 (SI unit: H). These must
be set to two different Inductor features in the circuit.
Transformer
The Transformer feature represents either a combination of two Inductor and a Mutual
Inductance features, or an ideal transformer.
NODE CONNECTIONS
Enter or edit the table in the Node names column for the primary and secondary node
connections.
DEVICE PARAMETERS
Choose a Transformer model — Specify inductors (the default) or Ideal transformer.
For Ideal transformer enter values or expressions for the Winding ratio N1/N2
(dimensionless). The default is 10.
Diode
The Diode device model ( ) is a large-signal model for a diode. It is an advanced
device model and no thorough description and motivation of the many input
parameters are attempted here. The interested reader is referred to Ref. 1 for more
details on semiconductor modeling within circuits. Many device manufacturers
NODE CONNECTIONS
Specify two Node names for the positive and negative nodes for the Diode device.
MODEL PARAMETERS
Specify the Model Parameters. Reasonable defaults are provided but for any particular
diode, the device manufacturer should be the primary source of information.
External I vs. U
The External I vs. U ( ) feature connects an arbitrary voltage measurement (for
example, a circuit terminal or circuit port boundary or a coil domain from another
physics interface) as a voltage source between two nodes in the electrical circuit. The
resulting circuit current from the first node to the second node is typically coupled
back as a prescribed current source in the context of the voltage measurement.
NODE CONNECTIONS
Specify the two Node names for the connecting nodes for the voltage source. The first
node represents the positive reference terminal.
EXTERNAL DEVICE
Enter the source of the Voltage. If circuit or current excited terminals or circuit ports
are defined on boundaries or domains or a multiturn coil domains is defined in other
physics interfaces, these display as options in the Voltage list. Also select the User defined
option and enter your own voltage variable, for example, using a suitable coupling
operator. For inductive or electromagnetic wave propagation models, the voltage
measurement must be performed as an integral of the electric field because the electric
External U vs. I
The External U vs. I ( ) feature connects an arbitrary current measurement (for
example, a coil domain from another physics interface) as a current source between
two nodes in the electrical circuit. The resulting circuit voltage between the first node
and the second node is typically coupled back as a prescribed voltage source in the
context of the current measurement.
NODE CONNECTIONS
Specify the two Node names for the connecting nodes for the current source. The
current flows from the first node to the second node.
EXTERNAL DEVICE
Enter the source of the Current. Voltage excited terminals or lumped ports defined on
boundaries in other physics interfaces are natural candidates but do not appear as
options in the Voltage list because those do not have an accurate built-in current
External I-Terminal
The External I-Terminal ( ) feature connects an arbitrary voltage-to-ground
measurement (for example, a circuit terminal from another physics interface) as a
voltage-to-ground assignment to a node in the electrical circuit. The resulting circuit
current from the node is typically coupled back as a prescribed current source in the
context of the voltage measurement. This node does not apply when coupling to
inductive or electromagnetic wave propagation models because then voltage must be
defined as a line integral between two points rather than a single point measurement
of electric potential. For such couplings, use the External I vs. U node instead.
NODE CONNECTIONS
Set the Node name for the connecting node for the voltage assignment.
EXTERNAL TERMINAL
Enter the source of the Voltage. If circuit- or current-excited terminals are defined on
boundaries in other physics interfaces, these display as options in the Voltage list. Also
Except when coupling to a circuit terminal, the current flow variable must
be manually coupled back in the electrical circuit to the context of the
voltage measurement. This applies also when coupling to a current
excited terminal. The name of this current variable follows the convention
cirn.termIm_i, where cirn is the tag of the Electrical Circuit interface
node and termIm is the tag of the External I-Terminal node. The tags are
typically displayed within curly brackets {} in the Model Builder.
See SPICE Import and Export about the supported SPICE commands.
See SPICE Export for more details on the supported SPICE commands.
• The field model is used to get a better, more accurate description of a single device
in the electrical circuit model.
• The electrical circuit is used to drive or terminate the device in the field model in
such a way that it makes more sense to simulate both as a tightly coupled system.
The Electrical Circuit interface makes it possible to add nodes representing circuit
elements directly to the Model Builder tree in a COMSOL Multiphysics model. The
circuit variables can then be connected to a physical device model to perform
co-simulations of circuits and multiphysics. The model acts as a device connected to
the circuit so that its behavior is analyzed in larger systems.
The fundamental equations solved by the Electrical Circuit interface are Kirchhoff’s
circuit laws, which in turn can be deduced from Maxwell’s equations. The supported
study types are Stationary, Frequency Domain, and Time Dependent.
Bipolar Transistors
Figure 5-1 illustrates the equivalent circuit for the npn bipolar junction transistor.
The pnp transistor model is similar in all regards to the npn transistor, with the
difference that the polarities of the currents and voltages involved are reversed. The
following equations are used to compute the relations between currents and voltages
in the circuit.
v be
--------------
v bc
--------------
– 1
NF VT NR VT
= ----------------------------------------------- 1 + 1 + 4I S ----------------------- + ------------------------
1 e –1 e
f bq
v v KF I A I A
2 1 – ----------- – -----------
bc be KR
V AF V AR
v be v be
I S -------------
N V
-
--------------
N V
i be = A ------- e F T – 1 + I SE e E T – 1
B
F
v bc v bc
I S --------------
N V --------------
N V
i bc = A -------- e R T – 1 + I SC e C T – 1
B
R
v be v bc
I S -------------
N V
-
N V
--------------
i ce = A ------- e F T + e C T
f
bq
k B T NOM
V T = ------------------------
q
There are also two capacitances that use the same formula as the junction capacitance
of the diode model. In the parameter names below, replace x with C for the
base-collector capacitance and E for the base-emitter capacitance.
v bx – MJx
1 – ---------
-
V Jx v bx < F C V Jx
C jbx = AC Jx ×
v bx v bx ≥ F C V Jx
( 1 – F ) – 1 – M Jx 1 – F ( 1 + M ) + M ---------
-
C C Jx Jx V
Jx
The following equations are used to compute the relations between currents and
voltages in the circuit.
C gd = C gd0 W
C gs = C gs0 W
1 – v
–M J
bd
--------
-
PB v bx < F C P B
C jbd = C BD ×
v bx v bx ≥ F C P B
( 1 – F ) – 1 – M J 1 – F ( 1 + M ) + M -------
-
C C J JP
B
The following equations are used to compute the relations between currents and
voltages in the circuit.
vd –M
1 – ------
- vd < FC VJ
V J
C j = C J0 ×
– 1 – M vd
( 1 – FC ) 1 – F C ( 1 + M ) + M ------- v d ≥ F C V J
V J
k B T NOM
V T = ------------------------
q
This chapter describes The Microwave Heating Interface found under the Heat
Transfer>Electromagnetic Heating branch ( ) when adding a physics interface.
See The Heat Transfer Interfaces and The Joule Heating Interface in the COMSOL
Multiphysics Reference Manual for other Heat Transfer interface and feature node
settings.
237
The Microwave Heating Interface
The Microwave Heating interface ( ) is used to model electromagnetic heating for
systems and devices that are on a scale ranging from 1/10 of a wavelength up to,
depending on available computer memory, about 10 wavelengths. This multiphysics
interface adds an Electromagnetic Waves, Frequency Domain interface and a Heat
Transfer in Solids interface. The multiphysics couplings add the electromagnetic losses
from the electromagnetic waves as a heat source, and the electromagnetic material
properties can depend on the temperature. The modeling approach is based on the
assumption that the electromagnetic cycle time is short compared to the thermal time
scale.
The following table shows what study step combinations of Electromagnetic Waves,
Frequency Domain interface and Heat Transfer in Solids interface are supported from
the Microwave Heating interface.
TABLE 6-1: STUDY STEP COMBINATIONS OF ELECTROMAGNETIC WAVES, FREQUENCY DOMAIN INTERFACE
AND HEAT TRANSFER IN SOLIDS INTERFACE
The Heat Transfer in Solids interface provides features for modeling heat transfer by
conduction, convection, and radiation. A Heat Transfer in Solids model is active by
default on all domains. All functionality for including other domain types, such as a
fluid domain, is also available. The temperature equation defined in solid domains
corresponds to the differential form of Fourier's law that may contain additional
contributions like heat sources.
However, if physics interfaces are added one at a time, followed by the coupling
features, these modified settings are not automatically included.
For example, if single Electromagnetic Waves, Frequency Domain and Heat Transfer in
Solids interfaces are added, the COMSOL adds an empty Multiphysics node. You can
choose Electromagnetic Heating from the available coupling features, but the modified
settings are not included.
Coupling features are available from the context menu (right-click the
Multiphysics node) or from the Physics toolbar, Multiphysics menu.
A side effect of adding physics interfaces one at a time is that four study
types — Frequency-Stationary; Frequency-Transient;
Frequency-Stationary, One-Way Electromagnetic Heating; and
Frequency-Transient, One-Way Electromagnetic Heating — are not
available for selection until after at least one coupling feature is added. In
this case, it is better to initially not add any study at all, then add the
coupling features to the Multiphysics node, and lastly, open the Add Study
window and add a study sequence below the Preset Studies for Selected
Multiphysics heading.
Use the online help in COMSOL Multiphysics to locate and search all the
documentation. All these links also work directly in COMSOL
Multiphysics when using the Help system.
Coupling Feature
The Electromagnetic Heating coupling feature node is described in this section.
• The available physics features for The Electromagnetic Waves, Frequency Domain
Interface are listed in the section Domain, Boundary, Edge, Point, and Pair Nodes
for the Electromagnetic Waves, Frequency Domain Interface.
• See The Heat Transfer Interfaces in the COMSOL Multiphysics Reference Manual
for information about the available physics features for heat transfer.
If you have an add-on module, such as the Heat Transfer Module, there
are additional specialized physics nodes available and described in the
individual module documentation.
Electromagnetic Heating
The Electromagnetic Heating node represents the electromagnetic losses, Qe (SI unit:
W/m3), as a heat source in the heat transfer part of the model. It is given by
Q e = Q rh + Q ml
1 *
Q rh = --- Re ( J ⋅ E )
2
1 *
Q ml = --- Re ( iωB ⋅ H )
2
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is emh.
DOMAIN SELECTION
When nodes are added from the context menu, you can select Manual (the default)
from the Selection list to choose specific domains to define the electromagnetic heat
source or select All domains as needed.
Only domains that are active in the physics interfaces selected in the Coupled Interfaces
section can be selected.
BOUNDARY SELECTION
When nodes are added from the context menu, you can select Manual (the default)
from the Selection list to choose specific boundaries to define the electromagnetic
boundary heat source or select All boundaries as needed.
Only boundaries that are active in the physics interfaces selected in the Coupled
Interfaces section can be selected.
COUPLED INTERFACES
This section defines the physics involved in the Electromagnetic Heating multiphysics
coupling. By default, the applicable physics interface is selected in the Electromagnetic
list to apply the Heat transfer to its physics interface to establish the coupling.
If a physics interface is deleted and then added to the model again, and in
order to re-establish the coupling, you need to choose the physics
interface again from the Heat transfer or Electromagnetic lists. This is
applicable to all multiphysics coupling nodes that would normally default
to the once present physics interface. See Multiphysics Modeling
Workflow in the COMSOL Multiphysics Reference Manual.
Glossary
245
Glossary of Terms
absorbing boundary A boundary that lets an electromagnetic wave propagate through
the boundary without reflections.
constitutive relation The relation between the D and E fields and between the B and
H fields. These relations depend on the material properties.
cutoff frequency The lowest frequency for which a given mode can propagate
through, for example, a waveguide or optical fiber.
eigenmode A possible propagating mode of, for example, a waveguide or optical fiber.
electric dipole Two equal and opposite charges +q and −q separated a short distance
d. The electric dipole moment is given by p = qd, where d is a vector going from −q
to +q.
lumped port A type of port feature. Use the lumped port to excite the model with a
voltage, current, or circuit input. The lumped port must be applied between two
metallic objects, separated by much less than a wavelength.
magnetic dipole A small circular loop carrying a current. The magnetic dipole
moment is m = IAe, where I is the current carried by the loop, A its area, and e a unit
vector along the central axis of the loop.
perfect electric conductor (PEC) A material with high electrical conductivity, modeled
as a boundary where the electric field is zero.
surface current density Current density defined on the surface. The component
normal to the surface is zero. The unit is A/m.
vector element A finite element often used for electromagnetic vector fields. The
tangential component of the vector field at the mesh edges is used as a degree of
freedom. Also called Nedelec’s edge element or just edge element.
INDEX| 249
boundary nodes dielectric medium theory 81
electromagnetic waves, frequency do- diffraction order (node) 123
main interface 100 diode (node) 221
electromagnetic waves, time explicit diode transistor model 233
166 dispersive materials 78
electromagnetic waves, transient 152 divergence constraint (node) 108
transmission line 159 documentation 19
bulk node 220 domain nodes
electromagnetic waves, frequency do-
C calculating
main interface 100
S-parameters 46
electromagnetic waves, time explicit
capacitor (node) 213
166
Cartesian coordinates 28
drain node 220
cementation exponent 148, 195
Drude-Lorentz dispersion model 106
circuit import, SPICE 225
Drude-Lorentz polarization (node) 156
circular port reference axis (node) 123
collector node 219 E E (PMC) symmetry 37
common settings 16 edge current (node) 146
complex permittivity, electric losses and effective medium (node) 148
185 eigenfrequency analysis 55
complex propagation constant 197 eigenfrequency calculations theory 189
complex relative permeability, magnetic eigenfrequency study 181
losses and 186 eigenmode analysis 83
constitutive relations 201 eigenvalue (node) 57
constitutive relations, theory 75 electric current density (node) 170
continuity, periodic boundaries and 33 electric field (node) 131, 170
coupling, to the electrical circuits inter- electric losses theory 185
face 61 electric point dipole (node) 146
curl-curl formulation 201 electric potential (node) 122
current source (node) 214 electric scalar potential 77
current-controlled current source electric susceptibility 202
(node) 217 electrical circuit interface 210
current-controlled voltage source theory 226
(node) 216 electrical circuits
cutoff frequency 84 modeling techniques 60
cylindrical coordinates 29 electrical conductivity 75
cylindrical waves 133 electrical conductivity, porous media 195
electrical size, modeling 13
D Debye dispersion model 106
electromagnetic energy theory 77
device models, electrical circuits 227
250 | I N D E X
electromagnetic heating (node) 241 Frequency-Domain Modal Method 68
electromagnetic quantities 88
G gate node 220
electromagnetic sources, applying 31
Gauss’ law 200
electromagnetic waves asymptotic scat-
geometry, simplifying 28
tering interface 176
ground (node) 122
electromagnetic waves, frequency do-
ground node (node) 212
main interface 92, 94
H H (PEC) symmetry 37
theory 180
high-frequency modeling 13
electromagnetic waves, time explicit in-
hybrid-mode waves
terface 165
axisymmetric, frequency domain 184
theory 200
axisymmetric, time domain 188
electromagnetic waves, transient inter-
in-plane, frequency domain 183
face 150
in-plane, time domain 187
theory 180
perpendicular 182
emailing COMSOL 21
emitter node 219 I impedance boundary condition (node)
error message, electrical circuits 61 137
exponential filter, for wave problems 168 importing
exporting SPICE netlists 64, 225
SPICE netlists 65 incoming wave (node) 161
external current density (node) 108 inductor (node) 213
external I vs. U (node) 222 inhomogeneous materials 78
external I-terminal (node) 224 initial values (node) 179
external U vs. I (node) 223 electromagnetic waves, frequency do-
main interface 108
F far field variables 39
electromagnetic waves, time explicit
Faraday’s law 200
interface 169
far-field calculation (node) 109, 174, 179
electromagnetic waves, transient 156
far-field calculations 86
transmission line 161
far-field domain (node) 109–110, 174
in-plane TE waves theory
far-field variables 37
frequency domain 184
file, Touchstone 99, 159
time domain 187
Floquet periodicity 33, 143
in-plane TM waves theory
fluid saturation 148
frequency domain 183
flux/source (node) 173
time domain 187
formation factor 195
inports 115
forward-propagating wave 197
integration line for current (node) 122
free-space variables 102
integration line for voltage (node) 122
frequency domain equation 181
INDEX| 251
internet resources 18 nonlinear materials 78
NPN bipolar junction transistor 227
K Kirchhoff’s circuit laws 226
NPN BJT (node) 219
knowledge base, COMSOL 22
numeric modes 114
L line current (out-of-plane) (node) 147
O open circuit (node) 162
linearization point 57
orthogonal polarization (node) 125
listener ports 115
losses, electric 185 P p-Channel MOS transistor 220
losses, magnetic 186 p-Channel MOSFET (node) 220
lossy eigenvalue calculations 55 PEC. see perfect electric conductor
lumped element (node) 130 perfect conductors theory 80
lumped port (node) 127, 163 perfect electric conductor (node) 171
lumped ports 48–49 boundaries 111
perfect magnetic conductor (node) 112,
M magnetic current (node) 144
171
magnetic current density (node) 170
periodic boundary conditions 33
magnetic field (node) 132, 171
periodic condition (node) 142
magnetic losses theory 186
periodic port reference point (node) 126
magnetic point dipole (node) 146
permeability
magnetic susceptibility 76, 202
anisotropic 182
material properties 78
permeability of vacuum 75
materials 79
permittivity
Maxwell’s equations 74
anisotropic 182
electrical circuits and 226
permittivity of vacuum 75
Maxwell-Ampere’s law 200
phasors theory 81
mesh resolution 32
physics interfaces, common settings 16
microwave heating interface 238
PMC. see perfect magnetic conductor
mode analysis 57, 182
PNP BJT (node) 219
mode phase
polarization, 2D and 2D axisymmetry 30
for Port and Diffraction Order 117
port (node) 113
modeling tips 27
port boundary conditions 46
MPH-files 20
ports, lumped 48–49
multiphysics couplings
potentials theory 77
microwave heating 238
power law, porous media
mutual inductance (node) 220
conductivity 192
N n-Channel MOS transistor 220, 230
permeability 194
n-Channel MOSFET (node) 220
permittivity 193
netlists, SPICE 64, 225
Poynting’s theorem 77
nodes, common settings 16
predefined couplings, electrical circuits
252 | I N D E X
61 port node and 113
propagating waves 197 theory 83
propagation constant 197 spherical waves 133
SPICE
Q quality factor (Q-factor) 55, 181
exporting 65
quasistatic modeling 13
SPICE netlists 64, 225
R reciprocal permeability, volume average
standard settings 16
194
study types 13
reciprocal permittivity, volume average
boundary mode analysis 57, 114
193
eigenfrequency 55, 181
reference point (node) 136
frequency domain 181
refractive index 103
mode analysis 57, 182
refractive index theory 185
subcircuit definition (node) 218
relative electric field 36
subcircuit instance (node) 219
relative permeability 76
surface charge density 80
relative permittivity 76
surface current density (node) 139, 172
remanent displacement 76
surface magnetic current density (node)
resistor (node) 212
139
S saturation coefficient 195 surface roughness (node) 140
saturation exponent 148 switch (node) 217
scattered fields, definition 35 symbols for electromagnetic quantities
scattering boundary condition (node) 88
132, 172 symmetry in E (PMC) or H (PEC) 37
scattering parameters. see S-parameters symmetry planes, far-field calculations 37
selecting symmetry, axial 29
mesh resolution 32
T TE axisymmetric waves theory
solver sequences 32
frequency domain 184
study types 13, 16
time domain 188
semiconductor device models 227
TE waves theory 85
short circuit (node) 163
technical support, COMSOL 21
SI units 88
TEM waves theory 85
simplifying geometries 28
terminating impedance (node) 162
skin effect, meshes and 32
theory
solver sequences, selecting 32
constitutive relations 75
source node 220
dielectrics and perfect conductors 80
space dimensions 14, 28
electric and magnetic potentials 77
S-parameter calculations
electrical circuit interface 226
electric field, and 45
electromagnetic energy 77
INDEX| 253
electromagnetic waves 180 voltage input, ports 48
electromagnetic waves, time explicit voltage source (node) 213
interface 200 voltage-controlled current source
far-field calculations 86 (node) 216
lumped ports 49 voltage-controlled voltage source (node)
Maxwell equations 74 215
phasors 81 voltmeter (node) 212
S-parameters 83 volume average, permeability 194
surface charges 80 volume average, permittivity 193
transmission line 196 volume averages, porous media 192
time domain equation, theory 186
W wave equation, electric (node) 102, 153
TM waves
wave equations (node) 167
axisymmetric 181
wave excitation 115
TM waves theory 85
wave impedance theory 85
Touchstone file 99, 159
wave number, free-space 102
transformer (node) 221
wavelength, meshes and 32
transition boundary condition (node)
websites, COMSOL 22
141
transmission line equation (node) 160
transmission line interface 158
theory 196
TW axisymmetric waves theory
frequency domain 184
time domain 188
two-port network (node) 144
two-port network port (node) 145
V variables
boundary mode analysis 58
eigenfrequency analysis and 56
far-field 37
for far fields 39
Jones vectors 53
lumped ports 51
mode analysis 58, 182
S-parameters 46
vector elements theory 188
254 | I N D E X