ACDCModuleUsersGuide PDF
ACDCModuleUsersGuide PDF
ACDCModuleUsersGuide PDF
User’s Guide
AC/DC Module User’s Guide
© 1998–2018 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; 8,954,302; 9,098,106; 9,146,652; 9,323,503; 9,372,673; and 9,454,625. 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 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.4
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
Fundamentals of Electromagnetics 32
Maxwell’s Equations. . . . . . . . . . . . . . . . . . . . . . 32
Constitutive Relations . . . . . . . . . . . . . . . . . . . . . 33
Potentials. . . . . . . . . . . . . . . . . . . . . . . . . . 35
Reduced Potential PDE Formulations . . . . . . . . . . . . . . . 35
Electromagnetic Energy . . . . . . . . . . . . . . . . . . . . 36
The Quasi-Static Approximation and the Lorentz Term . . . . . . . . 37
Material Properties . . . . . . . . . . . . . . . . . . . . . . 38
About the Boundary and Physics Interface Conditions . . . . . . . . . 39
Phasors . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Nonlinear Permanent Magnet Constitutive Relations . . . . . . . . . 41
Effective Nonlinear Magnetic Constitutive Relations . . . . . . . . . 41
The Jiles-Atherton Hysteresis Model . . . . . . . . . . . . . . . 42
References for Electromagnetic Theory . . . . . . . . . . . . . . 44
Electromagnetic Forces 45
Overview of Forces in Continuum Mechanics . . . . . . . . . . . . 45
Forces on an Elastic Solid Surrounded by Vacuum or Air. . . . . . . . 47
Torque . . . . . . . . . . . . . . . . . . . . . . . . . . 48
CONTENTS |3
Forces in Stationary Fields . . . . . . . . . . . . . . . . . . . 48
Forces in a Moving Body . . . . . . . . . . . . . . . . . . . . 52
Electromagnetic Energy, Coenergy and Virtual Work . . . . . . . . . 54
Electromagnetic Quantities 57
Coil Features 75
About the Coil Features . . . . . . . . . . . . . . . . . . . . 75
About the Coil Name . . . . . . . . . . . . . . . . . . . . . 78
Coil Excitation . . . . . . . . . . . . . . . . . . . . . . . . 78
Coil Groups. . . . . . . . . . . . . . . . . . . . . . . . . 82
Lumped Parameter Calculations . . . . . . . . . . . . . . . . . 82
The Single-Turn Coil Feature . . . . . . . . . . . . . . . . . . 84
Using Coils in 3D Models . . . . . . . . . . . . . . . . . . . 85
Automatic Coil Geometry Analysis . . . . . . . . . . . . . . . . 90
Accurate Coil Voltage Calculation . . . . . . . . . . . . . . . . 91
4 | CONTENTS
Lumped Parameters 93
Calculating Lumped Parameters with Ohm’s Law. . . . . . . . . . . 93
Calculating Lumped Parameters Using the Energy Method . . . . . . . 95
Studying Lumped Parameters . . . . . . . . . . . . . . . . . . 96
Lumped Parameter Conversion . . . . . . . . . . . . . . . . . 97
Stationary Source Sweep . . . . . . . . . . . . . . . . . . . . 99
CONTENTS |5
Space Charge Density . . . . . . . . . . . . . . . . . . . . 120
Force Calculation. . . . . . . . . . . . . . . . . . . . . . 120
Zero Charge . . . . . . . . . . . . . . . . . . . . . . . 121
Ground . . . . . . . . . . . . . . . . . . . . . . . . . 122
Electric Potential . . . . . . . . . . . . . . . . . . . . . . 122
Surface Charge Density . . . . . . . . . . . . . . . . . . . 123
External Surface Charge Accumulation . . . . . . . . . . . . . 123
Electric Displacement Field . . . . . . . . . . . . . . . . . . 124
Periodic Condition . . . . . . . . . . . . . . . . . . . . . 124
Thin Low Permittivity Gap . . . . . . . . . . . . . . . . . . 125
Dielectric Shielding . . . . . . . . . . . . . . . . . . . . . 126
Terminal . . . . . . . . . . . . . . . . . . . . . . . . . 127
Floating Potential . . . . . . . . . . . . . . . . . . . . . . 129
Distributed Capacitance . . . . . . . . . . . . . . . . . . . 131
Line Charge . . . . . . . . . . . . . . . . . . . . . . . . 131
Line Charge (on Axis) . . . . . . . . . . . . . . . . . . . . 132
Line Charge (Out-of-Plane) . . . . . . . . . . . . . . . . . . 132
Point Charge . . . . . . . . . . . . . . . . . . . . . . . 133
Point Charge (on Axis) . . . . . . . . . . . . . . . . . . . 134
Change Cross Section . . . . . . . . . . . . . . . . . . . . 134
Change Thickness (Out-of-Plane). . . . . . . . . . . . . . . . 135
Charge Conservation, Piezoelectric . . . . . . . . . . . . . . . 135
Electrostatic Point Dipole . . . . . . . . . . . . . . . . . . 136
Archie’s Law . . . . . . . . . . . . . . . . . . . . . . . 137
Effective Medium . . . . . . . . . . . . . . . . . . . . . . 137
6 | CONTENTS
The Electric Currents Interface 149
Domain, Boundary, Edge, Point, and Pair Nodes for the Electric
Currents Interface . . . . . . . . . . . . . . . . . . . . 151
Current Conservation . . . . . . . . . . . . . . . . . . . . 153
Initial Values . . . . . . . . . . . . . . . . . . . . . . . 155
External Current Density . . . . . . . . . . . . . . . . . . 155
Current Source . . . . . . . . . . . . . . . . . . . . . . 156
Electric Insulation . . . . . . . . . . . . . . . . . . . . . 157
Floating Potential . . . . . . . . . . . . . . . . . . . . . . 157
Boundary Current Source . . . . . . . . . . . . . . . . . . 159
Normal Current Density . . . . . . . . . . . . . . . . . . . 159
Distributed Impedance. . . . . . . . . . . . . . . . . . . . 160
Terminal . . . . . . . . . . . . . . . . . . . . . . . . . 162
Electric Shielding . . . . . . . . . . . . . . . . . . . . . . 164
Contact Impedance . . . . . . . . . . . . . . . . . . . . . 165
Electrical Contact . . . . . . . . . . . . . . . . . . . . . 167
Sector Symmetry . . . . . . . . . . . . . . . . . . . . . . 168
Line Current Source . . . . . . . . . . . . . . . . . . . . 169
Line Current Source (on Axis). . . . . . . . . . . . . . . . . 169
Point Current Source . . . . . . . . . . . . . . . . . . . . 170
Point Current Source (on Axis) . . . . . . . . . . . . . . . . 170
Electric Point Dipole . . . . . . . . . . . . . . . . . . . . 171
Electric Point Dipole (on Axis) . . . . . . . . . . . . . . . . 171
Piezoresistive Material . . . . . . . . . . . . . . . . . . . . 172
CONTENTS |7
Using the Extra Dimension Coordinates . . . . . . . . . . . . . 187
Boundary and Edge nodes for the Electric Currents, Layered Shell
Interface. . . . . . . . . . . . . . . . . . . . . . . . 188
Layered Shell . . . . . . . . . . . . . . . . . . . . . . . 189
Interlayer Connection . . . . . . . . . . . . . . . . . . . . 190
Homogenized Shell Connection . . . . . . . . . . . . . . . . 191
Homogenized Shell . . . . . . . . . . . . . . . . . . . . . 191
Varying Thickness . . . . . . . . . . . . . . . . . . . . . 193
Layer Current Source . . . . . . . . . . . . . . . . . . . . 193
Interface Current Source. . . . . . . . . . . . . . . . . . . 194
Electric Currents Solid to Shell Via Boundary . . . . . . . . . . . 195
8 | CONTENTS
Chapter 5: Magnetic Field Interfaces
CONTENTS |9
Domain Group/Boundary Group . . . . . . . . . . . . . . . . 251
Reversed Current Direction . . . . . . . . . . . . . . . . . 252
Lumped Port . . . . . . . . . . . . . . . . . . . . . . . 252
Lumped Element . . . . . . . . . . . . . . . . . . . . . . 255
Edge Current . . . . . . . . . . . . . . . . . . . . . . . 256
External Magnetic Vector Potential . . . . . . . . . . . . . . . 256
Impedance Boundary Condition . . . . . . . . . . . . . . . . 256
Transition Boundary Condition . . . . . . . . . . . . . . . . 258
Thin Low Permeability Gap . . . . . . . . . . . . . . . . . . 259
Magnetic Point Dipole . . . . . . . . . . . . . . . . . . . . 259
Magnetic Point Dipole (on Axis) . . . . . . . . . . . . . . . . 260
Magnetic Shielding . . . . . . . . . . . . . . . . . . . . . 260
10 | C O N T E N T S
The Magnetic Fields, No Currents, Boundary Elements
Interface 278
Domain and Boundary Nodes for the Magnetic Fields, No Currents,
Boundary Elements Interface . . . . . . . . . . . . . . . . 281
Magnetic Flux Conservation. . . . . . . . . . . . . . . . . . 281
Magnetic Insulation . . . . . . . . . . . . . . . . . . . . . 282
Initial Values . . . . . . . . . . . . . . . . . . . . . . . 282
Zero Magnetic Scalar Potential . . . . . . . . . . . . . . . . 282
Magnetic Scalar Potential . . . . . . . . . . . . . . . . . . . 283
Magnetic Flux Density . . . . . . . . . . . . . . . . . . . . 283
Magnetic Scalar-Scalar Potential Coupling . . . . . . . . . . . . 283
Magnetic Vector-Scalar Potential Coupling . . . . . . . . . . . . 285
CONTENTS | 11
Transient Equation . . . . . . . . . . . . . . . . . . . . . 309
Maxwell’s Equations. . . . . . . . . . . . . . . . . . . . . 310
Magnetic and Electric Potentials . . . . . . . . . . . . . . . . 310
Gauge Transformations . . . . . . . . . . . . . . . . . . . 311
Selecting a Particular Gauge . . . . . . . . . . . . . . . . . . 311
The Gauge and Equation of Continuity for Dynamic Fields . . . . . . 312
Explicit Gauge Fixing/Divergence Constraint . . . . . . . . . . . 312
Ungauged Formulations and Current Conservation . . . . . . . . . 313
Time-Harmonic Magnetic Fields . . . . . . . . . . . . . . . . 313
12 | C O N T E N T S
NPN BJT and PNP BJT. . . . . . . . . . . . . . . . . . . . 331
n-Channel MOSFET and p-Channel MOSFET . . . . . . . . . . . 332
Mutual Inductance . . . . . . . . . . . . . . . . . . . . . 333
Transformer . . . . . . . . . . . . . . . . . . . . . . . 333
Diode . . . . . . . . . . . . . . . . . . . . . . . . . . 334
External I vs. U . . . . . . . . . . . . . . . . . . . . . . 334
External U vs. I . . . . . . . . . . . . . . . . . . . . . . 335
External I-Terminal . . . . . . . . . . . . . . . . . . . . . 336
SPICE Circuit Import . . . . . . . . . . . . . . . . . . . . 337
SPICE Circuit Export . . . . . . . . . . . . . . . . . . . . 337
Chapter 8: Glossary
CONTENTS | 13
14 | C O N T E N T S
1
Introduction
This guide describes the AC/DC Module, an optional add-on package for
COMSOL Multiphysics® designed to assist you to solve and model low-frequency
electromagnetics.
In this chapter:
15
About the AC/DC Module
In this section:
The available physics interfaces cover the following types of electromagnetics field
simulations:
• Electrostatics
• Electric currents in conductive media
• Magnetostatics
• Low-frequency electromagnetics
Material properties include inhomogeneous and fully anisotropic materials, media with
gains or losses, and complex-valued material properties. Infinite elements makes it
possible to model unbounded domains. In addition to the standard results and
visualization functionality, the module supports direct computation of lumped
parameters such as capacitances and inductances as well as electromagnetic forces.
16 | CHAPTER 1: INTRODUCTION
Like all COMSOL modules, there is a library of ready-to-run models that make it
quicker and easier to analyze discipline-specific problems. In addition, any model you
develop is described in terms of the underlying partial differential equations, offering
a unique way to see the underlying physical laws of a simulation.
The AC/DC interfaces are fully multiphysics enabled — couple them to any other
physics interface in COMSOL Multiphysics or the other modules. For example, to find
the heat distribution in a motor, first find the current in the s using one of the
quasi-static interfaces in this module and then couple it to a heat equation in the main
COMSOL Multiphysics package or the Heat Transfer Module. This forms a powerful
multiphysics model that solves all the equations simultaneously.
The AC/DC Module also provides interfaces for modeling electrical circuits.
Each physics interface has a tag which is of special importance when performing
multiphysics simulations. This tag helps distinguish between physics interfaces and the
variables defined by the physics interface have an underscore plus the physics interface
tag appended to their names.
The Model Wizard is an easy way to select the physics interface and study type when
creating a model for the first time, and physics interfaces can be added to an existing
model at any time. Full instructions for selecting interfaces and setting up a model are
in the COMSOL Multiphysics Reference Manual.
See What Problems Can You Solve? and Table 1-1 for information about the available
study types and variables. See also Overview of the User’s Guide for links to the
chapters in this guide.
AC/DC
18 | CHAPTER 1: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
Electromagnetic Heating
20 | CHAPTER 1: INTRODUCTION
AC/DC Module Study Availability
TABLE 1-1: AC/DC MODULE DEPENDENT VARIABLES, FIELD COMPONENTS, AND PRESET STUDY AVAILABILITY
**
PHYSICS INTERFACE TAG DEPENDENT FIELD PRESET STUDIES
*
VARIABLES COMPONENTS
FREQUENCY-STATIONARY
SMALL-SIGNAL ANALYSIS,
FREQUENCY-TRANSIENT
MAGNETIC POTENTIAL
FREQUENCY DOMAIN
FREQUENCY DOMAIN
CURRENT DENSITY
TIME DEPENDENT
MAGNETIC FIELD
ELECTRIC FIELD
STATIONARY
Electrostatics es V x √ √ √ √
y
z
Electrostatics, esbe V x √
Boundary Elements y
z
Electric Currents ec V x x √ √ √ √
y y
z z
Electric Currents, ecs V x x √ √ √ √
Shell y y
z z
Magnetic Fields mf A x x x x √ √ √ √
y y y y
z z z z
Magnetic Field mfh H x √ √ √ √
Formulation y
z
Magnetic and mef V, A x x x x √ √ √
Electric Fields y y y y
z z z z
Magnetic Fields, No mfnc Vm x √ √
Currents y
z
Magnetic Fields, No mfncbe Vm x √
Currents, Boundary y
Elements z
FREQUENCY-STATIONARY
SMALL-SIGNAL ANALYSIS,
FREQUENCY-TRANSIENT
MAGNETIC POTENTIAL
FREQUENCY DOMAIN
FREQUENCY DOMAIN
CURRENT DENSITY
TIME DEPENDENT
MAGNETIC FIELD
ELECTRIC FIELD
STATIONARY
Rotating Machinery, rmm A, Vm x z z z √ √
Magnetic y
Electrical Circuit cir not √ √ √
applicable
Induction Heating — A,T, J x x x x √ √ √ √ √
y y y y
z z z z
*
These are the nonzero field components. For Cartesian coordinates, these are indexed by x, y, and z;
for cylindrical coordinates, r, ϕ, and z are used.
**
Custom studies are also available based on the physics interface, for example, Eigenfrequency and
Eigenvalue.
22 | CHAPTER 1: INTRODUCTION
In each module’s documentation, only unique or extra information is included;
standard information and procedures are centralized in the COMSOL Multiphysics
Reference Manual.
• 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.
24 | CHAPTER 1: INTRODUCTION
THE APPLICATION LIBRARIES WINDOW
Each model or application includes documentation with the theoretical background
and step-by-step instructions to create a model or app. The models and applications
are available in COMSOL Multiphysics as MPH files that you can open for further
investigation. You can use the step-by-step instructions and the actual models as
templates for your own modeling. In most models, SI units are used to describe the
relevant properties, parameters, and dimensions, but other unit systems are available.
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.
To include the latest versions of model examples, from the Help menu
select ( ) Update COMSOL Application Library.
26 | CHAPTER 1: INTRODUCTION
Overview of the User’s Guide
The AC/DC 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.
THEORY OF ELECTROMAGNETICS
In the Review of Electromagnetics chapter contains an overview of the theory behind
the AC/DC Module. It is intended for readers that wish to understand what goes on
in the background when using the physics interfaces and discusses the Fundamentals
of Electromagnetics, Electromagnetic Forces, and Electromagnetic Quantities.
ELECTRIC FIELDS
Electric Field Interfaces chapter describes these interfaces and includes the underlying
theory for each physics interface at the end of the chapter:
MAGNETIC FIELDS
Magnetic Field Interfaces chapter describes these interfaces and includes the
underlying theory for each physics interface at the end of the chapter:
• The Magnetic Fields Interface, which handles problems for magnetic fields with
source currents, solving for the magnetic vector potential and, for voltage driven s,
the scalar electric potential. Preset Stationary, Time Dependent, Frequency
Domain, and Small-Signal Analysis study types are available. This is the
recommended primary choice for modeling of magnetic fields involving source
currents.
• The Magnetic Field Formulation Interface has the equations, boundary conditions,
and currents for modeling magnetic fields, solving for the magnetic field. It is
especially suitable for modeling involving nonlinear conductivity effects, for
example in superconductors. Preset Stationary, Time Dependent, Frequency
Domain, and Small-Signal Analysis study types are available.
• The Magnetic Fields, No Currents Interface, which handles magnetic fields without
currents. When no currents are present, the problem is easier to solve using the
magnetic scalar potential. Thus it is the primary choice for modeling permanent
magnets. Preset Stationary and Time Dependent study types are available.
• The Magnetic Fields, No Currents, Boundary Elements Interface, which handles
magnetic fields without currents using the boundary element method. It can be
combined with the Magnetic Fields or Magnetic Fields, No Currents interfaces for
easy modeling of open boundary conditions. The preset Stationary study type is
available.
• The Rotating Machinery, Magnetic Interface combines a Magnetic Fields
formulation (magnetic vector potential) and Magnetic Fields, No Currents
formulation (magnetic scalar potential) with a selection of predefined frames for
prescribed rotation or rotational velocity - most of its features are taken either from
28 | CHAPTER 1: INTRODUCTION
the Magnetic Fields or the Magnetic Fields, No Currents interfaces. Preset
Stationary and Time Dependent study types are available.
• The Magnetic and Electric Fields Interface computes magnetic field and current
distributions when the exciting current is driven by an applied voltage. In most
cases, using the Magnetic Fields interface with its dedicated modeling features is the
preferred choice over using the Magnetic and Electric Fields interface. Preset
Stationary and Frequency Domain study types are available.
ELECTRICAL CIRCUIT
The Electrical Circuit Interface chapter describes the physics interface, which has the
equations for modeling electrical circuits with or without connections to a distributed
fields model, solving for the voltages, currents, and charges associated with the circuit
elements. The underlying theory for the physics interface is included at the end of the
chapter. Preset Stationary, Time Dependent and Frequency Domain study types are
available.
HEAT TRANSFER
Heat Transfer Interfaces chapter describes the physics interface, which combines all
physics features from the Magnetic Fields interface in the time harmonic formulation
with the Heat Transfer interface for modeling of induction and eddy current heating.
Heat transfer through conduction and convection in solids and free media (fluids) is
supported by physics interfaces shipped with the basic COMSOL Multiphysics license.
The Heat Transfer Interfaces, The Joule Heating Interface, and Theory
for Heat Transfer in the COMSOL Multiphysics Reference Manual.
Review of Electromagnetics
This chapter contains an overview of the theory behind the AC/DC Module. It is
intended for readers that wish to understand what goes on in the background when
using the physics interfaces.
In this chapter:
• Fundamentals of Electromagnetics
• Electromagnetic Forces
• Electromagnetic Quantities
• References for the AC/DC Interfaces
31
Fundamentals of Electromagnetics
In this section:
• Maxwell’s Equations
• Constitutive Relations
• Potentials
• Reduced Potential PDE Formulations
• Electromagnetic Energy
• The Quasi-Static Approximation and the Lorentz Term
• Material Properties
• About the Boundary and Physics Interface Conditions
• Phasors
• Nonlinear Permanent Magnet Constitutive Relations
• Effective Nonlinear Magnetic Constitutive Relations
• The Jiles-Atherton Hysteresis Model
• References for Electromagnetic Theory
Maxwell’s Equations
The problem of electromagnetic analysis on a macroscopic level is that of solving
Maxwell’s equations subject to certain boundary conditions. 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 equations can be formulated in differential form or integral form. The differential
form is presented here because it leads to differential equations that the finite element
∂D
∇ × H = J + -------
∂t
∂------
B-
∇×E = –
∂t
∇⋅D = ρ
∇⋅B = 0
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 equations include constitutive relations that describe
the macroscopic properties of the medium. They are given as:
D = ε0 E + P
B = μ0 ( H + M ) (2-1)
J = σE
where ε0 is the permittivity of vacuum, μ0 is the permeability of vacuum, and σ is the
electrical conductivity. In the SI system, the permeability of 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 electromagnetic constants ε0, μ0, and c0 are available in COMSOL Multiphysics
as predefined physical constants.
FUNDAMENTALS OF ELECTROMAGNETICS | 33
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.
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 instance,
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 χe E, where χe is the electric susceptibility. Similarly in linear materials, the
magnetization is directly proportional to the magnetic field, M = χm H, where χm is
the magnetic susceptibility. For such materials, the constitutive relations are:
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. Usually these are scalar properties but can, in the general case, be 3-by-3
tensors when the material is anisotropic. The properties ε and μ (without subscripts)
are the permittivity and permeability of the material, respectively.
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.
For some materials, there is a nonlinear relationship between B and H such that
Potentials
Under certain circumstances, it can be helpful to formulate the problems in terms of
the electric scalar potential V and the magnetic vector potential A. They are given by
the equalities:
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. In the
magnetostatic case where there are no currents present, Maxwell-Ampère’s law
reduces to ∇ × H = 0. When this holds, it is also possible to define a magnetic scalar
potential by the relation H = −∇Vm.
–1 dD
∇×( μ ∇×A ) = J +
dt
DOMAIN EQUATIONS
Time-Harmonic
For time-harmonic quasi-static systems solving for an A formulation, the reduced
potential formulation results in the following PDE:
FUNDAMENTALS OF ELECTROMAGNETICS | 35
2 –1
( jωσ – ω ε ) ( A ext + A red ) + ∇×( μ ∇×( A ext + A red ) ) = J e
Transient
Similarly to the time-harmonic formulation, in the transient formulation, the above
substitution results in the reduced equation
Static
In static formulations, the induced current is zero. Maxwell-Ampère’s law reduces to:
–1
∇×( μ ∇×( A ext + A red ) ) = J e
In this case it is also possible to express the external field through a known external
magnetic flux density, Bext. The domain equation in reduced form then reads:
–1
∇×( μ ( ∇×A red + B ext ) ) = J e
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 = V 0 H ⋅ dB dV = V 0 H ⋅ ------
∂t
- dt dV
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. 1)
The first term on the right-hand side represents the resistive losses,
Ph = V J ⋅ E dV
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 ∂ 1
H ⋅ ------- = --- B ⋅ ------- = ------- B ⋅ B
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), the
result is:
∂
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μ
FUNDAMENTALS OF ELECTROMAGNETICS | 37
fields are always delayed relative to the changes of the sources, reflecting the finite
speed of propagation of electromagnetic waves. Under the assumption that this effect
can be ignored, it is possible to obtain the electromagnetic fields by considering
stationary currents at every instant. This is called the quasi-static approximation. The
approximation is valid provided that the variations in time are small and that the
studied geometries are considerably smaller than the wavelength (Ref. 5).
The quasi-static approximation implies that the equation of continuity can be written
as ∇ ⋅ J = 0 and that the time derivative of the electric displacement ∂D/∂t can be
disregarded in Maxwell-Ampère’s law.
There are also effects of the motion of the geometries. Consider a geometry moving
with velocity v relative to the reference system. The force per unit charge, F/q, is then
given by the Lorentz force equation:
F
---- = E + v × B
q
This means that to an observer traveling with the geometry, the force on a charged
particle can be interpreted as caused by an electric field E' = E + v × B. In a conductive
medium, the observer accordingly sees the current density
J = σ ( E + v × B ) + Je
∇ × H = σ ( E + v × B ) + Je
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 for which they require some consideration. A given material
can belong to one or more of these groups.
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 two. Of these four conditions, only two
are independent. This is an overdetermined system of equations, so it needs to be
reduced. First select either equation one or equation four. Then select either equation
two or equation three. Together these selections form a set of two independent
conditions.
From these relationships, the interface condition is derived for the current density,
FUNDAMENTALS OF ELECTROMAGNETICS | 39
∂ρ 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, 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
If the fields satisfy a linear time-dependent equation, then the corresponding phasors
must satisfy a similar equation in which the time derivatives are replaced by a factor jω.
All time-harmonic equations in the AC/DC Module are expressed as equations for the
phasors. (The tilde is dropped from the variable denoting the phasor.)
FUNDAMENTALS OF ELECTROMAGNETICS | 41
The formulation still applies the assumption of harmonic time dependency of the fields
(which are still described as phasors): the fields still oscillate at the given frequency;
higher-order harmonic effects, or other similar effects due to the nonlinearity, are not
accounted for.
B = μ0 ( H + M )
• Interdomain coupling: α
• Maximum (saturation) magnetization: Ms
• Langevin slope: a
These quantities are second-order tensors, something that is not reflected by the
equation notation here and in the graphical user interface.
The interdomain coupling parameter is used for defining the effective magnetic field
He, which is the field that is assumed to affect the magnetic moments on a
microscopical level.
He = H + α ⋅ M
The maximum magnetization and Langevin slope determines the anhysteretic part of
the material behavior, which is supposed to follow the Langevin function
Lan ( x ) = cot h ( x ) – 1 ⁄ x
–1 He
M an = M s ⋅ Lan ( a ⋅ H e ) ----------
He
Ms and α enter the expression with the nonlinear Langevin function and consistent
with [4], they are available as diagonal matrices. The use of the Langevin function
should be interpreted as applied component wise.
The remaining parameter, k, is controlling the pinning and hysteresis of the magnetic
moments.
Given the parameters above, the key equation in the Jiles-Atherton model,
determining the magnetization M, is
dH e χ dM an
= max χ ⋅
dM
, 0 ⋅ ------ + c ⋅
dt dt χ dt
–1
χ = k ⋅ ( M an – M )
The Jiles-Atherton hysteresis model is available mainly for time-dependent studies, but
it also supports stationary simulation. The latter is especially useful if the material is
FUNDAMENTALS OF ELECTROMAGNETICS | 43
starting from a magnetized configuration that is not the solution of a previous
transient study. In the graphical user interface, it is possible to select the initial value of
magnetization, in which case the initial magnetic field is set accordingly.
The internal state variables are discretized directly in the Gauss points, and an auxiliary
variable is added to the same domain also for aiding support to simulations with
moving domains.
The order of Gauss point discretization can be selected in the Discretization section
under the Magnetic Field section. For the Gauss point discretization order, the default
is Automatic. This default sets the Gauss points discretization order depending on the
order of the parent physics. If the parent physics variable has order n, the Gauss point
order is 2*(n–1). It is good practice to test that the results do not change significantly
when using a Gauss point discretization of order 2*n since a good discretization order
depends on the kind of magnetic potential, mesh element shape, and space dimension,
rather than the pure mesh size. For an example, see.
2. J. Jin, The Finite Element Method in Electromagnetics, John Wiley & Sons, New
York, 1993.
The derivation of the expressions for the electromagnetic stress tensor utilizes
thermodynamic potential (energy) principles (Ref. 1 and Ref. 3). The distribution of
electromagnetic forces in a system depends on the material. Accordingly, the
techniques and expressions used when calculating electromagnetic forces are different
for different types of materials.
Another technique for calculating forces using the method of virtual work is described
in the section Electromagnetic Energy, Coenergy and Virtual Work.
In this section:
2
d r
ρ 2
= ∇ ⋅ T + f ext
dt
where ρ is the density, r denotes the coordinates of a material point, T is the stress
tensor, and fext is an external volume force such as gravity (fext = ρg). This is the
equation solved in the structural mechanics interfaces for the special case of a linear
elastic material, neglecting the electromagnetic contributions.
ELECTROMAGNETIC FORCES | 45
In the stationary case there is no acceleration, and the equation representing the force
balance is
0 = ∇ ⋅ T + f ext
The stress tensor must be continuous across a stationary boundary between two
materials. This corresponds to the equation
n1 ( T2 – T1 ) = 0
where T1 and T2 represent the stress tensor in Materials 1 and 2, respectively, and n1
is the normal pointing out from the domain containing Material 1. This relation gives
rise to a surface force acting on the boundary between Material 1 and 2.
Material 2
Material 1
n1
In certain cases, the stress tensor T can be divided into one part that depends on the
electromagnetic field quantities and one part that is the mechanical stress tensor,
T = T EM + σ M
For the special case of an elastic body, the mechanical stress tensor is proportional only
to the strain and the temperature gradient. The exact nature of this split of the stress
tensor into an electromagnetic and a mechanical part depends on the material model,
if it can be made at all.
For more information on the mechanical stress tensor for elastic materials,
see the documentation for the interfaces. For example, The Solid
Mechanics Interface in the COMSOL Multiphysics Reference Manual.
It is sometimes convenient to use a volume force instead of the stress tensor. This force
is obtained from the relation
f em = ∇ ⋅ T EM
0 = ∇ ⋅ σ M + f em + f ext
–∇ ⋅ σM = f where f = f em + f ext
In practice, the equation for the force balance also needs to include an external
boundary force gext. It is nonzero on those parts of the boundary where it is necessary
to compensate for the contributions to the stress tensor that you are not interested in
or do not have enough information on. These contributions come from the influence
of the adjacent domains. By approximating the surroundings by vacuum or air, the
influence of these boundaries and their adjacent domains (that are not part of our
model) on the electromagnetic fields are neglected.
˜
n1 ( T2 – T1 ) = 0
˜
n 1 T 2 = n 1 T 2 + g ext
The external boundary force gext can represent the reaction force from another body
that the solid is attached to.
The equations for the balance of forces on the solid now become
∇ ⋅ T 1 + f ext = 0
n 1 ( T 2 – T 1 ) + g ext = 0
For calculating the total force F on the solid these equations need to be integrated over
the entire solid and the solid/vacuum boundary
( ∇ ⋅ T1 + fext ) dV + ° ( n1 ( T2 – T1 ) + gext ) dS = 0
Ω1 ∂Ω 1
∇ ⋅ T1 dV – ° n1 T1 dS = 0
Ω1 ∂Ω 1
ELECTROMAGNETIC FORCES | 47
this means that the external force
is needed to balance the term for the boundary integral of the stress tensor in the
surrounding vacuum
F = n1 T2 dS
°
∂Ω 1
to keep the solid stationary. That is Fext + F = 0. If the external forces are suddenly
removed, the solid is no longer stationary, but F causes the solid to begin to move with
an initial acceleration according to
2
d r
ma = ρ ---------
dt
2
dV = F
Ω1
Torque
The torque in the case of Forces on an Elastic Solid Surrounded by Vacuum or Air is
given by
MO =
° ( r – rO ) × ( n1 T2 ) dS
∂Ω 1
where rO is a point on the axis of rotation. This follows from a derivation similar to
the one made for forces.
that is, if the fields vary so slowly that the contributions from induced currents and
displacement currents can be neglected.
Also assume that the objects modeled are not moving v = 0 so that there is no
contributions from Lorentz forces. These are treated later on.
To apply the stress tensor in air to calculate the total force and torque on
a magnetizable rod close to a permanent magnet, see Permanent
Magnet: Application Library path ACDC_Module/Magnetostatics/
permanent_magnet.
The expressions for the stress tensor in a general electromagnetic context stems from
a fusion of material theory, thermodynamics, continuum mechanics, and
electromagnetic field theory. With the introduction of thermodynamic potentials for
mechanical, thermal, and electromagnetic effects, explicit expressions for the stress
tensor can be derived in a convenient way by forming the formal derivatives with
respect to the different physical fields (Ref. 1 and Ref. 3). Alternative derivations can
be made for a vacuum (Ref. 4) but these cannot easily be generalized to polarized and
magnetized materials.
ε0
T 2 = – pI – ----- E ⋅ E + --------- B ⋅ B I + ε 0 EE + ------ BB
1 T 1 T
2 2μ 0 μ0
where p is the air pressure, I is the identity 3-by-3 tensor (or matrix), and E and B are
3-by-1 vectors. In this expression of the stress tensor, air is considered to be
nonpolarizable and nonmagnetizable. When air is approximated by vacuum, p = 0.
This expression, with p = 0, of the stress tensor is also known as the Maxwell stress
tensor.
Using the fact that, for air, D = ε0E and B = μ0H the expression for the stress tensor
can be written as
ELECTROMAGNETIC FORCES | 49
T 2 = – pI – --- E ⋅ D + --- H ⋅ B I + ED + HB
1 1 T T
2 2
Maxwell’s equations in free space give that the contribution of the electromagnetic
part of the stress tensor is zero, and the resulting expression is
0 = – ∇p + f ext
Thus, using the same terminology as earlier, fem = 0 for air, with σM = −pI. In the
derivation of the total force on an elastic solid surrounded by vacuum or air, the
approximation ∇p = 0 has been used.
When operating with the divergence operator on the stress tensor, the relation
∇ ⋅ EE – --- E ⋅ EI = E ( ∇ ⋅ E ) – E × ( ∇ × E )
T 1
2
is useful (and similarly for B). From the right-hand side it is clear (using Maxwell’s
equations) that this is zero for stationary fields in free space.
Consider again the case of a solid surrounded by air. To compute the total force, the
projection of the stress tensor on the outside of the solid surface is needed,
n 1 T 2 = – pn 1 – --- E ⋅ D + --- H ⋅ B n 1 + ( n 1 ⋅ E )D + ( n 1 ⋅ H )B
1 1 T T
2 2
where n1 is the surface normal, a 1-by-3 vector, pointing out from the solid. This
expression can be used directly in the boundary integral of the stress tensor for
calculating the total force F on the solid.
For an example of how to compute the total force on two parallel wires
either by integrating the volume force or by integrating the stress tensor
on the surrounding surface, see Electromagnetic Forces on Parallel
Current-Carrying Wires: Application Library path ACDC_Module/
Verification_Examples/parallel_wires.
T 1 = σ M – --- E ⋅ D + --- H ⋅ B I + ED + HB
1 1 T T
2 2
The situation is slightly different from the case of air because there can be currents and
volume charges in the conductor. The current density is
1
J = ∇ × H = ------ ∇ × B
μ0
ρ = ∇ ⋅ D = ε0 ∇ ⋅ E
0 = ∇ ⋅ σ M + ρE + J × B + f ext
f em = ρE + J × B
ε0
T 1 = σ ( E, B ) – ----- E ⋅ E + --------- B ⋅ B – M ⋅ B I
1
2 2μ 0
T 1 T T T
+ ε 0 EE + ------ BB + EP – MB
μ0
where in σ(E, B) the dependence of E and B has not been separated out. Thus σ is
not a purely mechanical stress tensor in this general case. Different material models
give different appearances of σ(E, B). The electromagnetic contributions to σ(E, B)
typically represent pyroelectric, pyromagnetic, piezoelectric, piezomagnetic, dielectric,
ELECTROMAGNETIC FORCES | 51
and magnetization effects. The expression for the stress tensor in vacuum, air, and pure
conductors can be derived from this general expression by setting M = P = 0.
T1 must be symmetric. The terms EPT and −MBT are symmetric in the case of a linear
dielectric and magnetic material because
P = ε0 χe E
M = χB B
Here, the magnetic susceptibility χB differs slightly from the classical χm. The other
explicit terms are all symmetric, as is σ(E, B). In the general case this imposes
constraints on the properties of σ(E, B). For a nonlinear material σ(E, B) might need
to include terms such as −EPT or +MBT to compensate for asymmetric EPT or
−MBT.
To instantiate the stress tensor for the general elastic case, an explicit material model
including the magnetization and polarization effects is needed. Such material models
can easily be found for piezoelectric materials (Ref. 3).
r' = r + v 0 t
The electromotive intensity is the most important of these invariants. The Lorentz
magnetization is significant only in materials for which neither the magnetization M
nor the polarization P is negligible. Such materials are rare in practical applications.
The same holds for the magnetization term of the magnetomotive intensity. Notice
that the term ε0v × E is very small compared to B/μ0 except for cases when v and E
are both very large. Thus in many practical cases this term can be neglected.
T 2 = – pI – --- E ⋅ D + --- H ⋅ B I + ED + HB + ( D × B )v
1 1 T T T
2 2
T 1 = σ M – --- E ⋅ D + --- H ⋅ B I + ED + HB + ( D × B )v
1 1 T T T
2 2
To get the equation for the balance of forces the divergence of this expression needs
to be computed. This requires an introduction of an extra term in Cauchy’s equation
corresponding to an additional electromagnetic contribution to the linear momentum.
Cauchy’s equation with this extra term is
2
d r
ρ 2
+ D × B = ∇ ⋅ T + f ext
dt
The extra term is canceled out by the additional term in the stress tensor, and the final
result is
ELECTROMAGNETIC FORCES | 53
2
d r ˜ ˜
ρ 2
= ∇ ⋅ σ M + ρE + J × B + f ext
dt
For the case of no acceleration, with the explicit appearance of the transformed
quantities,
0 = ∇ ⋅ σ M + ρ ( E + v × B ) + ( J – ρv ) × B + f ext
The terms containing v × B cancel out, which yields the following equation:
0 = ∇ ⋅ σ M + ρE + J × B + f ext
˜ ε0 ˜
T 1 = σ ( E, B ) – ----- E ⋅ E + --------- B ⋅ B – M ⋅ B I +
1
2 2μ 0
T 1 T ˜ T ˜ T T
+ ε 0 EE + ------ BB + E P – M B + ε 0 ( E × B )v
μ0
To instantiate the stress tensor for the general elastic case a material model explicitly
including the magnetization and polarization effects is needed.
The method of virtual work is used for the electric energy and magnetic energy
separately for calculating the total electric or magnetic force as follows.
Knowing either of these, the other can be computed through the relation
dW m = H ⋅ B – dW' m
For the case of linear materials, the magnetic energy and coenergy densities are equal.
The method of virtual work utilizes that, under constant magnetic flux conditions
(Ref. 5), the total magnetic force on a system is computed as
F Φ = – ∇W m
∂W m
TΦ = –
∂ϕ
Under the condition of constant currents and linear materials, the total force and
torque can be computed in the same way but with opposite signs,
F I = ∇W m
∂W m
TI =
∂ϕ
Under the condition of constant currents and nonlinear materials, the total force
and torque must be computed from the magnetic coenergy and with opposite signs,
F I = ∇W' m
∂W' m
TI =
∂ϕ
ELECTROMAGNETIC FORCES | 55
D
dW e = E ⋅ dD'
0
E
dW' e = D ⋅ dE'
0
Knowing either of these, the other can be computed through the relation
dW e = E ⋅ D – dW' e
For the case of linear materials, the electric energy and coenergy densities are equal.
Under the condition of constant charges, the total electric force and torque on a
system are computed as
F Q = – ∇W e
∂W e
TQ = –
∂ϕ
Under the condition of constant potentials and linear materials, the total electric
force and torque on a system can be computed as
F V = ∇W e
∂W e
TV =
∂ϕ
Under the condition of constant potentials and nonlinear materials, the total
electric force and torque on a system must be computed from the electric coenergy as
F V = ∇W' e
∂W' e
TV =
∂ϕ
The method of virtual work can be employed by using the features for
deformed mesh and sensitivity analysis in COMSOL Multiphysics.
ELECTROMAGNETIC QUANTITIES | 57
TABLE 2-1: ELECTROMAGNETIC QUANTITIES
4. R.K. Wangsness, Electromagnetic Fields, 2nd ed., John Wiley & Sons, 1986.
5. D.K. Cheng, Field and Wave Electromagnetics, 2nd ed., Addison-Wesley, 1991.
The goal of this chapter is to familiarize you with the modeling procedure in the
AC/DC Module. Because this module is fully integrated with COMSOL
Multiphysics®, the modeling process is similar.
In this chapter:
61
Preparing for Modeling
This section is intended a guide you through the selection process among the
interfaces in the AC/DC Module and does not contain detailed physics interface
descriptions. Several topics about the art of modeling are covered here that might not
be in ordinary textbooks on electromagnetic theory. This section discusses these
topics:
• What Problems Can You Solve? — Can I use the quasi-static physics interfaces or
do I need wave propagation?
• Selecting the Space Dimension for the Model Geometry — Is a 2D, 3D, or
axisymmetric geometry best for my model?
• Simplifying the Geometry Using Boundary Conditions — When do I need to
resolve the thickness of thin shells?
• Applying Electromagnetic Sources — What sources can I use to excite the fields?
• Selecting a Study Type — Is my problem suited for time-dependent or
time-harmonic (frequency domain) formulations?
• 2D Field Variables
• About the AC/DC and Nonlinear Magnetic Material Databases
• Meshing and Solving — What issues might arise with respect to meshing and
solving?
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.
It can be more accurate and efficient to use several simple models instead
of a single, complex one.
One major difference between quasi-static and high-frequency modeling is that the
formulations depend on the electrical size of the structure. This dimensionless
measure is the ratio between the largest distance between two points in the structure
divided by the wavelength of the electromagnetic fields.
The quasi-static physics interfaces in this module are suitable for simulations of
structures with an electrical size in the range up to 1/10. The physical assumption of
these situations is that the currents and charges generating the electromagnetic fields
vary so slowly in time that the electromagnetic fields are practically the same at every
instant as if they had been generated by stationary sources.
When the variations in time of the sources of the electromagnetic fields are more rapid,
it is necessary to solve the full Maxwell equations for high-frequency electromagnetic
waves. They are appropriate for structures of electrical size 1/100 and larger. Thus, an
Independent of the structure size, the AC/DC 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.
Examples of applications that successfully simulate with this module include electric
motors, generators, permanent magnets, induction heating devices, and dielectric
heating. For a more detailed description of some of these applications, refer to the
Application Libraries included with the module.
This is because 2D models are easier to modify and solve much faster. Thus, modeling
mistakes are much easier to find when working in 2D. Once the 2D model is verified,
you are in a much better position to build a 3D model.
2D PROBLEMS
Cartesian Coordinates
In this case you view a cross section in the xy-plane of the actual 3D geometry. The
geometry is mathematically extended to infinity in both directions along the z-axis,
There are usually two approaches that lead to a 2D cross-section view of a problem:
The geometry has a finite width but the model neglects the (end) effects from the faces
parallel to the cross section because the strongest forces are between the perpendicular
faces (those seen as lines in the cross section).
Figure 3-1: The cross sections and their real geometry for Cartesian coordinates and
cylindrical coordinates (axial symmetry).
The spatial coordinates are called r and z, where r is the radius. The flow at the
boundaries is given per unit length along the third dimension. Because this dimension
3D PROBLEMS
Is it 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.
See 2D Problems if this is applicable.
Are there symmetries in the geometry and model? Many problems have planes
where the solution on either side of the plane looks the same. A good way to check this
is to flip the geometry around the plane, for example, by turning it upside down
around the horizontal plane. You can then remove the geometry below the plane if you
do not see any differences 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.
Can you replace the thin layers with boundary conditions? There are several types of
boundary conditions in COMSOL Multiphysics suitable for such replacements. You
can, for example, replace materials with high conductivity with the shielding boundary
condition, which assumes a constant potential through the thickness of the layer. If you
have a magnetic material with a high relative permeability, you can also model it using
the shielding boundary condition.
Use boundary conditions for known solutions. A body with a high conductivity at high
frequency has the current density confined to a thin region beneath the surface of the
wire. You can often replace the current in the body by either a surface current
boundary condition or an impedance boundary condition.
For nonlinear problems, use a frequency domain study after a linearization of the
problem, which assumes that the distortion of the sinusoidal signal is small.
Specify a Time Dependent study when you think that the nonlinear influence is very
strong, or if you are interested in the harmonic distortion of a sinusoidal signal. It
might also be more efficient to use a Time Dependent study if there is a periodic input
with many harmonics, like a square-shaped signal.
There are some special predefined study types for the Induction Heating interface.
This physics interface is based on the assumption that the magnetic cycle time is short
compared to the thermal time scale (adiabatic assumption). Thus, it is associated with
four predefined study sequences.
• Frequency-Stationary
- Time-harmonic magnetic fields
- Stationary heat transfer
• Sequential Frequency-Stationary
- Time-harmonic magnetic fields
- Stationary heat transfer
• Sequential Frequency-Transient
- Time-harmonic magnetic fields
- Transient heat transfer
The difference between the first two and the last two study sequences is that in the
Frequency-Stationary and the Frequency-Transient studies, the coupling between the
electromagnetic and the heat transfer problem is bidirectional. The electromagnetic
problem defines a heat source for the heat transfer problem and the material properties
in the electromagnetic problem are temperature dependent. For the Sequential
Frequency-Stationary and Sequential Frequency-Transient study sequences, it is
assumed that the material properties for the electromagnetic problem is temperature
independent, so a Frequency Domain study is added for the electromagnetic problem
followed by a Stationary or Time Dependent study for the heat transfer problem.
2D Field Variables
When solving for a vector field in 2D, the physics interface has three options: to solve
for the out-of-plane vector, the in-plane vector, or the three-component vector.
Depending on the choice, the available source specification options on the domain,
boundary, edge, and point levels change accordingly.
Nonlinear magnetic materials are generally provided as interpolation functions for the
magnetization curve without hysteresis effects. In the main equation nodes, that is in
any of the nodes; Ampère’s Law, Ampère’s Law and Current Conservation, Faraday’s Law,
or Magnetic Flux Conservation and, depending on which of the AC/DC interfaces is
being used, one needs to select the magnetic field constitutive relation B-H curve or H-B
curve.
MESH RESOLUTION
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:
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.
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. You need at least two linear
elements per skin depth to capture the variation of the fields. If you do not study the
skin depth, you can replace regions with a small skin depth with a boundary condition,
thereby saving elements.
SELECTING A SOLVER
In most cases the solver that COMSOL Multiphysics suggests can be used. The choice
of solver is optimized for the typical case for each physics interface and study type in
the AC/DC Module. However, in special cases the solver settings might need fine
tuning. This is especially important for 3D problems because they use a large amount
of memory. For large 3D problems, a 64-bit platform might be required.
• Meshing
• Studies and Solvers
1 T
n 1 T 2 = – --- n 1 ( E ⋅ D ) + ( n 1 ⋅ E )D (3-1)
2
1 T
n 1 T 2 = – --- n 1 ( H ⋅ B ) + ( n 1 ⋅ H )B
2
is integrated on the surface to obtain the force. In the Magnetic and Electric Fields
interface, both expressions are included. E is the electric field, D the electric
displacement, H the magnetic field, B the magnetic flux density, and n1 the outward
normal from the object. For a theoretical discussion about the stress tensor, see
Electromagnetic Forces.
LORENTZ FORCES
The Lorentz force is defined as F = J × B. The Lorentz force is very accurate for
electromagnetic force calculations in electrically conducting domains. The Lorentz
CYCLE-AVERAGED QUANTITIES
In Frequency Domain studies, the formulas presented above can be used to compute
the instantaneous force and torque by taking the real part of all the variables involved.
It is often more interesting to compute the average of these quantities over a cycle to,
for example, couple the electromagnetic forces to another physics. The cycle-averaged
Maxwell’s stress tensor (including both electric and magnetic forces) is computed as:
The Permanent Magnet model demonstrates how to compute the total force on a
magnetizable rod close to a permanent magnet by integrating the Maxwell stress
tensor in the air on the outside of the rod. This is the most important method for
accurately calculating the total force on magnetic devices for which the exact
distribution of volume forces is not known. To retrieve the exact distribution of
volume forces requires a material that describes the interactions of the magnetizations
and strains. Such materials are not always available. Therefore you are often limited to
compute the total force by integrating the stress tensor or using the method of virtual
work.
The Coil is available both as a domain and as a boundary feature. The latter can be
used to model flat coils for which the thickness is negligible compared to the other two
dimensions.
COIL FEATURES | 75
The Coil feature supports two different Conductor model options, which differ by the
physical system represented and the modeling details.
• The Single conductor option models a single, solid region of a conducting material
(for example, metal) in which the current flows. The lumped voltage and current of
the coil correspond respectively to the integral of the electric field along the coil
length and to the integral of the current density on a cross section. Use this option
in a domain feature to model a single wire with a non-negligible cross section; in a
boundary feature it can be used for a thin conductive layer whose thickness is
negligible (for example, in frequency domain, it is small with respect to the skin
depth).
• The Homogenized multi-turn option implements a homogenized model of a coil
consisting of numerous tightly-wound conducting wires, separated by an electrical
insulator. The computation of the voltage and current of the coil is performed in a
similar way as for the Single conductor model, but it also takes into account
parameters such as the number of wires and the cross-section area. Use this feature
to model a coil containing a large number of wires without the need to model each
wire individually. The boundary feature can be used when the wires are arranged in
a thin region whose thickness is negligible. Note that capacitive coupling and skin
effect are neglected in the homogenized model.
For 3D model geometries, the direction of the current flow is not easily
determined. The coil domains have settings and subnodes to solve this
problem. See Using Coils in 3D Models for more information.
• When using the Single conductor option, the coil models a solid conductive domain
(typically a metal). The material properties of the metal (electric permitting,
magnetic permeability, and electrical conductivity) should be used in this case.
• When using the Homogenized multi-turn option, the conduction current (induced
current) is assumed to flow only in the wires. To prevent induced current density
from flowing in the domain, the region’s electrical conductivity σ is set to zero. The
material used in the domain, only specifies the constitutive relations for the electric
displacement field and the magnetic flux density. Therefore, choose or create a
material that provides constitutive relations representative of the entire coil domain,
not just of the material making up the metallic wires.
COIL FEATURES | 77
• The RLC Coil Group feature, that models coils for which capacitive effects between
coil turns and with other parts of the models are comparable with the inductive
effects.
• The Single-Turn Coil feature, to introduce an externally applied current or voltage
on domains where the Ampère’s Law and Current Conservation feature is active.
This feature is contributing, that is, it simply adds an external current density to the
domain model represented by Ampère’s Law and Current Conservation.
Coil Excitation
The current density flowing in the coil domain is computed from a lumped quantity
that constitutes the coil excitation. The choice of this quantity can be done by setting
the coil excitation parameter, while the value is specified in the associated text box. All
the coil features can be excited either with current excitation, voltage excitation or,
in certain cases, power excitation. The supplied value or expression is translated to a
current density or electric field applied to the domain according to the coil model used.
Coil features can also retrieve the value of the current or the voltage from a circuit
connection, if the parameter coil excitation is set to Circuit (current) or Circuit
(voltage), respectively.
CURRENT EXCITATION
where L is equal to the physics interface thickness d for 2D models, 2πr for 2D axially
symmetric models, and the coil length in 3D models. V is an unknown applied
potential. The potential V is solved for using an additional equation, which constrains
the total integrated current to be equal to the current value Icoil specified, acting as a
controller.
NI coil
J e = --------------- (3-3)
A
σV
J e = -------- (3-4)
L
where L is equal to the physics interface thickness d for 2D models and 2πr for 2D
axially symmetric models, and V is an unknown applied potential. The potential V is
solved for using an additional algebraic equation, which constrains the total integrated
current to be equal to the current value Icoil specified.
VOLTAGE EXCITATION
V coil
E e = ----------- (3-5)
L
COIL FEATURES | 79
where Vcoil is the applied voltage specified, and L is equal to the physics interface’s
thickness d for 2D models, to 2πr for 2D axially symmetric models, and the coil length
in 3D models.
N ( V coil + V ind )
J e = --------------------------------------- (3-6)
AR coil
where Vcoil is the applied voltage which is specified, A is the total cross-sectional area
of the coil domain, N is the number of turns specified, Rcoil is the total resistance of
the coil calculated as
NL
R coil = σ--------------------------
coil a coil A
- (3-7)
A
where L is equal to the physics interface’s thickness d for 2D models, 2πr for 2D axially
symmetric models, and the coil length in 3D models. The expression σcoilacoil is the
product of the wire bulk conductivity and wire cross-section area. Vind is the induced
voltage calculated by integrating the electric field along the coil.
V coil = Vi
i=1
where Vcoil is the user-defined voltage drop across the coil, Vi is the individual applied
potentials, and N is the number of turns in the coil group.
σV coil
J e = --------------- (3-8)
L
POWER EXCITATION
In 2D and 2D axisymmetric models, Power is an option available for the coil excitation
parameter. When Power is selected, the physics interface sets up the equation as
specified in the Current Excitation section, plus a constraint between the coil current
and voltage:
1
--- real ( V coil I coil∗ ) = P coil
2
in frequency domain (cycle average). Select this option to specify the input power for
the coils. The problem becomes nonlinear with these settings.
The Power option is not available for 3D models due to the increased
complexity of the problem compared to 2D models.
Be aware that, in general, the values of the voltage and the current might
not uniquely be determined by this constraint. Particularly in the
frequency domain, the absolute phase of the quantities can be arbitrary.
This indeterminacy can have an impact on the solution process. See
Compile Equations in the COMSOL Multiphysics Reference Manual for
a possible solution to this problem.
• Coil
• Single-Turn Coil
• Studies and Solvers in the COMSOL Multiphysics Reference Manual
COIL FEATURES | 81
Coil Groups
The Coil nodes, in 2D and 2D axisymmetry, present an additional Coil Group check box
that allows to easily model one or more coil features connected in series. Each domain
in the coil feature selection is considered a separate turn of the coil. Use this feature in
2D and 2D axisymmetric models for coils that cross the modeling plane more than
once.
The Coil Group mode automatically sets up the relations between the variables
associated to each domain. Current and voltage variables for each domain (or turn) are
computed the same as for the base coil feature, while the total variables for the coil are
computed using the rules for series-connection—the total voltage is the sum of the
voltages, while the total current is equal to the current flowing in each domain.
Special subnodes can be added to the Coil nodes once the Coil Group check box has
been selected in order to better control the set-up of the connections:
Domain Group
When the Coil Group mode is active, each domain in the selection will be considered
individually as a series-connected entity. The Domain Group subnode can be used to
mark a group of domains which should be treated as one domain for the purpose of
the coil set up. Multiple instances of this subnodes can be added to mark different
groups. Domains that do not belong to a group are connected individually in series.
V coil
R coil = -----------
I coil
VARIABLE EQUATION
Impedance V coil
Z coil = -----------
I coil
Admittance
Y coil = Z coil -1
COIL FEATURES | 83
lindev ( V coil )
Z coil = -----------------------------------
lindev ( I coil )
In COMSOL Multiphysics version 5.2 and earlier, the Single-Turn Coil was
available in The Magnetic Fields Interface and The Rotating Machinery,
Magnetic Interface. The functionality provided by this feature has been
incorporated in the Coil feature in the form of the Single conductor option.
COIL FEATURE
The options available in a Coil feature depend on the chosen Conductor model:
• When Single conductor is used, the coil direction is computed in the dedicated
preprocessing study step Coil Geometry Analysis.
• When Homogenized multi-turn is used, the Coil type input determines how the
direction of the current flow is specified. The available alternatives are: Linear,
Circular, Numeric and User defined. The Numeric option computes the current flow
in the dedicated preprocessing study step Coil Geometry Analysis. The other three
options do not require any preprocessing.
In the Homogenized multi-turn case, the coil type determines how the direction of the
wires constituting the coil is specified, as well as the method used to compute the
average length and the average cross-section area of the domain or boundary. The coil
length and coil area are used to compute lumped variables, such as the induced voltage
or the total resistance.
COIL FEATURES | 85
For models representing only a part of the geometry (due to symmetry), the total
length and cross-section area of the coil (referred to as coil length and coil area) are
computed by multiplying the domain or boundary length and area by the appropriate
multiplication factors specified under Symmetry specification in the subfeatures.
The following sections detail the possible options available for the Coil feature.
The geometry analysis algorithm can determine the axis and direction only if there is
a geometry in the model. If the model does not have a geometry, for example when
using imported meshes, use the alternative analysis method by selecting the Use robust
geometry analysis method check box. This method can be used even without a
geometry, but it requires that the selected edges form a complete circle to work
correctly. If neither the default method nor the alternative method work, an alternative
is to set the Coil type to User defined and manually specify the direction of the wires
using a Cylindrical System.
By default, the Input (for Geometry Analysis) boundary condition is available under
the Geometry Analysis subnode. For open coils (whose ends are on exterior
boundaries), apply this condition on the input boundary, that is, the boundary at
which the wires enter the coil domain. Right-click the Geometry Analysis node, add
an Output (for Geometry Analysis) subnode, and apply it on the exterior boundary
where the wires exit the domain.
For closed-loop coils, apply the Input boundary condition on an interior boundary
orthogonal to the direction of the wires.
COIL FEATURES | 87
To complete the set up, add a Coil Geometry Analysis study step to the study, before the
main study step.
The Coil Geometry Analysis can process multiple coils in the same step.
In order to analyze only certain coils, select the Solve only specific coils
check box and enter the coil names as a comma-separated list.
Since the coil feature does not have any information on the geometry, additional
settings are required:
• Enter the Domain length (for domain features) or the Boundary length (for boundary
feature), which should be the average length of the coil from the input to the
output—the average length of a wire, not accounting for symmetry nor the number
of turns.
• For domain features, use the subnodes to User Defined Coil Geometry — Input and
Output — to select the coil input and output boundaries, where the current flow
enters or leaves the simulation domain. If the coil is closed in a loop, no selection is
required.
SINGLE-TURN COIL
The 3D Single-Turn Coil feature was available in The Magnetic Fields Interface in
previous versions of COMSOL Multiphysics. This feature is obsolete — the
recommended approach to model the same physical system is to use a Coil feature with
the Single conductor model, or The Magnetic and Electric Fields Interface.
The feature is still available for backward compatibility when opening models created
in previous versions of the software.
This node represents a solid conducting domain, typically a wire or a coil, with a
non-negligible cross section. The boundary node represents a conducting thin layer
whose thickness is small (also compared to the skin depth). To enforce the current
conservation in the domain, an additional dependent variable with the dimension of
an electric potential (SI unit: V) is added to the problem and the continuity equation
for the current is introduced in the system of equations. This variable is referred to as
the coil potential, but it is only loosely related to the electrostatic potential and it
should be considered a help variable rather than representing a tangible physical
quantity. In the Single-Turn Coil node, it is possible to specify the material properties
that are used in the continuity equation.
A Gap Feed subnode models a thin gap in the conductive domain across which a
difference of potential or a current is applied. This feature should be applied on
internal boundaries to the conductive domain and is useful for modeling closed loops.
For the Single-Turn Coil feature, only one active feed feature can be used
at a time. That is, either a Boundary Feed or Gap Feed feature can be
used, not both.
COIL FEATURES | 89
The Gap Feed should be used with care in high frequency modeling as
there is a mesh dependent displacement current density flowing in the
elements just outside the bounding edges of the Gap Feed boundary. If
these elements have nonzero conductivity, this problem can appear even
at lower frequencies.
The study step can be used to compute both the current flow when using the Single
conductor model and the direction of the wires when using the Homogenized multi-turn
model.
In the case of a Single conductor coil, the study step solves a simple current conservation
problem to determine the flow of a current applied to the metallic domain.
In the case of a Homogenized multi-turn coil, the feature computes a vector field e
which represents the local wire density in the coil, as well as the length and average
cross-section of the wires. The vector variable eCoil can be plotted (for example, in a
Streamline or Arrow Volume plot) to visualize the computed direction of the wires.
The Coil Geometry Analysis study step must precede the main study step
(for example, a Stationary study step) in which the Coil is used.
The coil geometry analysis method tries to construct a coil wires density that is
plausible for an actual bundle of thin wires. The computed wire density is constant over
arbitrary cross section, without violating current conservation for the current flowing
The Segregated solver uses Geometric Multigrid to solve the linear system associated
with the analysis. If the Geometric Multigrid causes problems during the solution (for
example, in presence of very coarse meshes), the direct solver can be used instead.
Change the Linear solver in all Segregated Steps nodes to Direct.
COIL FEATURES | 91
as the coil impedance, inductance, and power, also benefit from the improved accuracy.
The functionality is enabled by default in frequency domain studies, and can be
manually controlled using the corresponding check box in the Coil Settings window,
Stabilization section. The Stabilization section can be shown by clicking the Show
button ( ) and selecting Stabilization, and is only available when the Conductor model
is set to Homogenized multi-turn.
The Accurate coil voltage calculation functionality performs a filtering of the external
current density applied by the coil to obtain a current density vector field that is
numerically divergence free. The filtering is performed by introducing an additional
scalar degree of freedom Vf (the filtering variable) and solving the equation
∇ ⋅ ( ∇V f – J u ) = 0
in the coil domain, where Ju is the unfiltered external current density applied by the
Coil feature. The filtered current density
J e = ∇V f – J u
is then applied as right-hand side in the magnetic equations. This filtering reduces the
spurious electric fields in the solution, in turn improving significantly the accuracy of
the induced voltage (which is computed by integrating the electric field in the coil
domain).
The filtering equation requires appropriate boundary conditions for the filtering
variable, in particular, on the inlet and outlet boundaries of the coil, where the external
current density is orthogonal. The approach used to set up these boundary conditions
vary according to the coil geometry:
• For Linear and Circular coils, the geometry analysis algorithm automatically detects
the inlet/outlet boundaries from the geometry. If this detection fails, use the
Override Coil Inlet/Outlet subnode to manually select the inlet/outlet boundaries.
Right-click the Coil node to add this subnode.
• For Numeric coils, the subnodes of the Geometry Analysis node are used to identify
the inlet/outlet boundaries.
• For User defined coils, use the subfeatures of the User Defined Coil Geometry
subnode — Input and Output — to select the coil inlet/outlet boundaries.
In this section:
There are several available techniques to extract the lumped parameters. Which one to
use depends on the physics interface, the parameter of interest, and how the model is
solved. The overview of the techniques in this section use a 4-by-4 matrix example for
the lumped parameter matrix. This represents a system of at least five electrodes, where
four are used as terminals and the rest are grounded, as illustrated in Figure 3-2.
V1 V3
Ground
V2 V4
Figure 3-2: A five-electrode system with 4 terminals and one ground electrode.
LUMPED PARAMETERS | 93
If a system specifies that all electrodes are terminals, the results are redundant matrix
elements. This is better understood by considering a two-electrode system. If both
electrodes are declared as terminals, a 2-by-2 matrix is obtained for the system. This is
clearly too many elements because there is only one unique lumped parameter between
the terminals. If in addition one or more ground electrodes are declared, the system
has three unique electrodes and the lumped parameter matrix becomes a 2-by-2
matrix.
FORCED VOLTAGE
If voltages are applied to the terminals, the extracted currents represent elements in the
admittance matrix, Y. This matrix determines the relation between the applied voltages
and the corresponding currents with the formula
I1 Y 11 Y 12 Y 13 Y 14 V 1
I2 Y 21 Y 22 Y 23 Y 24 V 2
=
I4 Y 31 Y 32 Y 33 Y 34 V 3
I4 Y 41 Y 42 Y 43 Y 44 V 4
so when V1 is nonzero and all other voltages are zero, the vector I is proportional to
the first column of Y.
In electrostatics the current is replaced with charge and the admittance matrix is
replaced with the capacitance matrix
Q1 C 11 C 12 C 13 C 14 V 1
Q2 C 21 C 22 C 23 C 24 V 2
=
Q4 C 31 C 32 C 33 C 34 V 3
Q4 C 41 C 42 C 43 C 44 V 4
FIXED CURRENT
It might be necessary to calculate the Z-matrix in a more direct way. Similar to the Y
calculation, the Z calculation can be done by forcing the current through one terminal
at the time to a nonzero value while the others are set to zero. Then, the columns of
the impedance matrix are proportional to the voltage values on all terminals:
In magnetostatics this option means that the energy method is used; see Calculating
Lumped Parameters Using the Energy Method below.
FIXED CHARGE
The Electrostatics interface can use total charge instead of total current. This gives the
inverted capacitance matrix in a similar manner as the Z and Y matrices.
–1
V1 C 11 C 12 C 13 C 14 Q1
V2 C 21 C 22 C 23 C 24 Q2
=
V3 C 31 C 32 C 33 C 34 Q4
V4 C 41 C 42 C 43 C 44 Q4
2 0 j≠i
Vi Ω
C ii = ------2- W e dΩ Vj =
Vi j = i
0 k ≠ i, j
1 Vi Vj
= ------------ W e dΩ – --- ------ C + ------ C jj
1
C ij
Vi Vj 2 V j ii V i
Vk = Vi
k = i
Ω
Vj k = j
Calculate the inductance matrix in the same way from the magnetic energy density:
LUMPED PARAMETERS | 95
2 0 j≠i
Ii Ω
L ii = ----2- W m dΩ Ij =
Ii j = i
0 k ≠ i, j
1 Ii Ij
= -------- W m dΩ – --- ---- L + ---- L
1
L ij
Ii Ij 2 I j ii I i jj
Ik = Ii
k = i
Ω
Ij k = j
The lumped parameters are defined as global variables. Evaluate these from the Derived
Values node under Results in the Model Builder or define 1D plot groups.
ACCURACY
Use reaction terms to be accurate when calculating the total current over the
boundary. This is necessary for the forced voltage input property. The reaction terms
(representing current or charge density) come from default information stored in the
–1 –1
S = G ref ⋅ ( E – ( Z ref∗ ⋅ Y ) ) ⋅ ( E + Z ref ⋅ Y ) ⋅ G ref ,
–1 Im ( Z ) Im ( Y )
Z = Y , L = ----------------- , C = ----------------- ,
ω ω
R = Re ( Z ), G = Re ( Y )
1 0 0 0
E = 0 1 0 0
0 0 1 0
0 0 0 1
Z ref = E ⋅ Z 0
1
G ref = E ⋅ ------------------------------
2 Re ( Z 0 )
For capacitance matrix calculations special transformations are available to obtain the
mutual or SPICE capacitance matrix from the basic output from the solver. The latter
is, for voltage excitation the so called Maxwell capacitance matrix and, for charge
excitation the inverse Maxwell capacitance matrix. The Maxwell capacitance values,
terminal charges, and terminal voltages are linked by the following matrix relation:
LUMPED PARAMETERS | 97
Q1 C 11 C 12 C 13 V 1
Q 2 = C 21 C 22 C 23 V 2 (3-9)
Q3 C 31 C 32 C 33 V 3
In order to analyze the response of a device, it is more common to use the form of a
mutual capacitance matrix. The mutual capacitance matrix can be converted from the
Maxwell capacitance matrix:
where CMij are elements of the Maxwell capacitance matrix and Cmij are elements of
the mutual capacitance matrix.
1 Set each diagonal element of the mutual capacitance matrix to the sum of the
corresponding row elements of the Maxwell capacitance matrix (that is, to get Cm11,
sum CM1j for all j).
2 Reverse the sign of all of the off-diagonal elements in the Maxwell capacitance
matrix to get the corresponding mutual capacitance matrix elements.
This conversion can also be used to transform resistance and conductance matrices
obtained from the Electric Currents physics. In that context the output resistance matrix
corresponds to the inverse Maxwell capacitance matrix, electrode currents correspond
to charges and the mutual or SPICE conductance matrix corresponds to the mutual or
SPICE capacitance matrix.
Matrix Output
The resulting potential values at the terminals are used to extract one of the lumped
matrices according to the following table.
During postprocessing, other matrices can be produced via the Global Matrix
Evaluation. The settings for Global Matrix Evaluation has a Transformation section that
allows for the output of the Maxwell capacitance matrix or the mutual capacitance
matrix.
If the Generate default plots check box on the Study node is enabled,
default Global Matrix Evaluation nodes for matrix output are generated.
LUMPED PARAMETERS | 99
Solver Selection
The Stationary Source Sweep is optimized for solution speed. It works best with a
direct solver as the LU decomposition has to be performed only once, followed by
back substitution for the multiple source excitations (multiple right-hand sides). For
3D models, the Stationary Source Sweep generates both an iterative and a direct solver.
The former will be selected if the model is having a large number of elements; that is,
a number that is above a built-in threshold. For hardware with sufficiently large RAM,
it is recommended to manually switch to the direct solver.
Voltage Output
With unit charge feeding the absolute value of the electric potential is difficult to
predict. For this reason there is an auxiliary global variable — for example,
es.VexcTerm — that represents the value of potential of the excited terminal. The
default potential plots will typically show the electric potential scaled by this value.
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.
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
= – ( J s ⋅ a h ) dl
w w
where the integration is taken in the direction of ah × n. This gives the following
relation between the transmission line impedance and the surface impedance
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.
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 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 lumped ports with the numbers 1 and 2 and Lumped 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.
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
• A choice is made in the Settings window for the non-circuit 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 non-circuit 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.
Some modeling errors lead to the error message The DAE is structurally
inconsistent, being displayed when solving. This 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.
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/
This chapter summarizes the functionality of the electric field interfaces found
under the AC/DC branch ( ). The AC/DC Module enhances the Electrostatics
and Electric Currents interfaces included with the basic COMSOL Multiphysics®
license.
In this chapter:
113
The Electrostatics Interface
The Electrostatics (es) interface ( ), found under the AC/DC branch ( ) when
adding a physics interface, is used to compute the electric field, electric displacement
field, and potential distributions in dielectrics under conditions where the electric
charge distribution is explicitly prescribed. The formulation is stationary except for use
together with other physics interfaces. Eigenfrequency, frequency-domain,
small-signal analysis, and time-domain modeling are supported in all space dimensions.
The physics interface solves Gauss’ Law for the electric field using the scalar electric
potential as the dependent variable.
Charge Conservation is the main node, which adds the equation for the electric
potential and has a Settings window for defining the constitutive relation for the
electric displacement field and its associated properties such as the relative permittivity.
When this physics interface is added, these default nodes are also added to the Model
Builder — Charge Conservation, Zero Charge (the default boundary condition), and
Initial Values. Then, from the Physics toolbar, add other nodes that implement, for
example, boundary conditions and space charges. You can also right-click Electrostatics
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 es.
PHYSICS-CONTROLLED MESH
Select the Enable check box to allow the physics interface to control the meshing
process. Information from the physics, such as the presence of an infinite elements
domain or periodic condition, will be used to set up automatically an appropriate mesh
sequence.
Select the Activate manual terminal sweep check box to switch on the sweep and invoke
a parametric sweep over the terminals. Enter a Sweep parameter name to assign a
specific name to the variable that controls the terminal number solved for during the
sweep. The Sweep parameter name must also be declared as a model parameter. The
default is PortName.
The generated lumped parameters are in the form of capacitance matrix elements. The
terminal settings must consistently be of either fixed voltage or fixed charge type.
The lumped parameters are subject to Touchstone file export. Enter a file path or Browse
for a file. Select a Parameter format (value pairs) for the Touchstone export—Magnitude
and angle (MA) (the default), Magnitude in dB and angle (DB), or Real and imaginary
parts (RI). Select an option from the If file exists list — Overwrite or Create new. Create
new is useful when the model is solved multiple times with different settings. Select a
Parameter to export: Z (the default), Y, or S.
From the Home toolbar add a Parameters node and enter the chosen name and assign
to it a temporary Expression of unity into the Parameters table. You can also right-click
Global Definitions to add the node.
DEPENDENT VARIABLES
The dependent variable is the Electric potential V. You can change its name, which
changes both the field name and the variable name. If the new name coincides with the
name of another electric potential field in the model, the physics interfaces shares
degrees of freedom. The new name must not coincide with the name of a field of
another type or with a component name belonging to some other field.
DISCRETIZATION
In the COMSOL Multiphysics Reference Manual, see Table 2-3 for links
to common sections and Table 2-4 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.
• Domain, Boundary, Edge, Point, and Pair Nodes for the Electrostatics
Interface
• Theory for the Electrostatics Interface
In the absence of surface charges, this condition is fulfilled by the natural boundary
condition
n ⋅ [ ( ε 0 ∇V – P ) 1 – ( ε 0 ∇V – P ) 2 ] = – n ⋅ ( D 1 – D 2 ) = 0
AVAILABLE NODES
These nodes, listed in alphabetical order, 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). Also see Table 4-1 for a list of interior and exterior
boundary conditions, including edge, point, and pair availability.
In the COMSOL Multiphysics Reference Manual, see Table 2-3 for links
to common sections and Table 2-4 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.
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation, and other variables
evaluated in a material reference configuration (material frame). Select Non-solid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
ELECTRIC FIELD
Select a Constitutive relation to describe the macroscopic properties of the medium
(relating the electric displacement D with the electric field E) and the applicable
material properties, such as the relative permittivity. Select:
• Relative permittivity (the default) to use the constitutive relation D = ε0εrE. Then
the default is to take the Relative permittivity εr (dimensionless) values From material.
For User defined, select Isotropic, Diagonal, Symmetric, or Anisotropic and enter values
or expressions in the field or matrix. The default is 1. If Effective medium is selected,
the Effective Medium subnode is available from the context menu (right-click the
parent node) as well as from the Physics toolbar, Attributes menu.
• Polarization to use the constitutive relation D = ε0E + P. Then enter the components
based on space dimension for the Polarization vector P (SI unit: C/m2). The
defaults are 0 C/m2.
• Remanent electric displacement to use constitutive relation D = ε0εrE + Dr, where Dr
is the remanent displacement (the displacement when no electric field is present).
Then the default is to take the Relative permittivity εr (dimensionless) values From
material. For User defined, select Isotropic, Diagonal, Symmetric, or Anisotropic and
enter values or expressions in the field or matrix. Then enter the components based
on space dimension for the Remanent electric displacement Dr (SI unit: C/m2). The
defaults are 0 C/m2.
INITIAL VALUES
Enter a value or expression for the initial value of the Electric potential V (SI unit: V).
The default value is 0 V.
Force Calculation
Use the Force Calculation node to define globally available force and torque variables
for the selected domains. Force Calculation is also available and defined here for the
Magnetic Fields, Magnetic and Electric Fields, and Magnetic Field Formulation
interfaces.
FORCE CALCULATION
Enter a Force name, which is then appended to global variables. The method used to
compute forces and torques is integration of the Maxwell’s stress tensor over the
exterior surfaces of the set of domains. This feature also gives access to the normal
component of the Maxwell stress tensor on the external surfaces. For the Magnetic and
Electric Fields and Magnetic Fields interfaces, the force calculation includes both
electric and magnetic forces. The computed force is made available as a global vector
variable with name es.Force<component>_<force name>.
Enter a direction vector for the Torque axis rax and coordinates for the Torque rotation
point r0. A torque calculation about a given point (Torque rotation point) is made and
defined as a global vector variable, es.T<component>_<force name>. The resulting
torque component parallel to the given Torque axis is given as a global scalar variable,
typically es.Tax_<force name>.
For more information on calculation of forces and torques, see Force and
Torque Computations.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
The Compute both instantaneous and cycle-averaged variables in frequency domain check
box is selected by default. The check box controls which variables are defined by the
Force Calculation feature in frequency-domain studies. When selected (the default
setting in new models), the variables are defined as described above. This check box
selected by default for backward compatibility. If the check box is cleared, only the
cycle-averaged variables are computed, replicating the behavior of the Force
Calculation feature in version 4.3 and earlier. The check box is cleared in existing Force
Calculation features when opening models created in version 4.3 and earlier.
Zero Charge
The Zero Charge node adds the condition that there is zero charge on the boundary so
that n ⋅ D = 0. This boundary condition is also applicable at symmetry boundaries
where the potential is known to be symmetric with respect to the boundary. This is the
default boundary condition at exterior boundaries. At interior boundaries, it means
that no displacement field can penetrate the boundary and that the electric potential is
discontinuous across the boundary.
Ground means that there is a zero potential on the boundary. This boundary condition
is also applicable at symmetry boundaries where the potential is known to be
antisymmetric with respect to the boundary.
For some physics interfaces, also select additional Ground nodes from the Edges (3D
components) or Points (2D and 3D components) submenus. For 2D axisymmetric
components, it can be applied on the Symmetry axis.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Electric Potential
The Electric Potential node provides an electric potential V0 as the boundary condition
V = V0.
Because the electric potential is being solved for in the physics interface, the value of
the potential is typically defined at some part of the geometry. For some physics
interfaces, also select additional Electric Potential nodes from the Edges (3D
components) or Points (2D and 3D components) submenus. For 2D axisymmetric
components, it can be applied on the symmetry axis.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
–n ⋅ D = ρs , n ⋅ ( D1 – D2 ) = ρs
–n ⋅ D = ρs
dρ s
= n ⋅ Ji + n ⋅ Je
dt
where n·Ji is the normal component of the total ion current density on the wall and
n·Je is the normal component of the total electron current density on the wall, which
are feature inputs.
n ⋅ D = n ⋅ D0
Periodic Condition
The Periodic Condition node defines periodicity or antiperiodicity between two
boundaries. If required, activate periodic conditions on more than two boundaries, in
which case the Periodic Condition tries to identify two separate surfaces that can each
consist of several connected boundaries. For more complex geometries, it might be
necessary to use the Destination Selection subnode. With this subnode the boundaries
which constitute the source and destination surfaces can be manually specified. The
BOUNDARY SELECTION
When using nonconforming meshes on the source and destination of a periodic
boundary pair, for numerical stability, a finer mesh should be applied on the
destination side. Use conforming meshes if possible.
PERIODIC CONDITION
Select a Type of periodicity — Continuity (the default), Antiperiodicity, or Floquet
periodicity. Select:
• Continuity to make the electric potential periodic (equal on the source and
destination),
• Antiperiodicity to make it antiperiodic, or
• Floquet periodicity (only available with products supporting piezoelectric modeling).
Specify the components of the k-vector for Floquet periodicity kF (SI unit: rad/m).
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
For information about the Orientation of Source section, see Orientation of Source
and Destination in the COMSOL Multiphysics Reference Manual.
ε0 εr
n ⋅ D 1 = ---------- ( V 1 – V 2 )
ds
ε0 εr
n ⋅ D 2 = ---------- ( V 2 – V 1 )
ds
Dielectric Shielding
The Dielectric Shielding node adds dielectric shielding as a boundary condition. It
describes a thin layer with thickness ds and a bulk relative permittivity; εr that shields
the electric field:
n ⋅ D = – ∇ t ⋅ ε 0 ε rs d s ∇tV
Use this boundary condition when approximating a thin domain with a boundary to
reduce the number of mesh elements.
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation and other variables
evaluated in a material reference configuration (material frame). Select Non-solid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
ELECTRIC SHIELDING
The default is to take the Relative permittivity εr (dimensionless) values From material.
It takes it from the adjacent domains if not explicitly defined. For User defined select
Terminal
The Terminal node provides a boundary or domain condition for connection to
external circuits, to transmission lines, or with a specified voltage or charge. By
specifying zero charge, a floating potential condition is obtained.
The Terminal node also enables the computation of the lumped parameters of the
system, such as capacitance. See Lumped Parameters.
TERMINAL
Specify the terminal properties. Enter a unique Terminal name that will be used to
identify the global variables (such as the terminal voltage or current) that are defined
by the terminal. The Terminal name can contain alphanumeric characters, but it should
be numeric when performing terminal sweeps.
Select a Terminal type — Charge (the default), Voltage, Circuit, or Terminated. Select:
• Charge to prescribe the total charge deposited on the terminal boundaries. Enter the
total charge Q0. The default is 0 C, which corresponds to a uncharged floating
electrode.
The Initial value for voltage is available for all Terminal types except Voltage. Enter an
Initial value for voltage Vinit.
The Initial value for charge is available for the Circuit excitation. It specifies the initial
charge on the terminal when the circuit is connected. Enter an Initial value for charge
Qinit.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Floating Potential
The Floating Potential node is used when modeling a metallic electrode at floating
potential. It applies a constant voltage V0 on the boundary (for domain features, this
is the boundary enclosing the selected domain), such that the total normal electric
displacement field D equals a specific charge Q0:
V = V0
( D ⋅ n ) dS = Q0 ,
∂Ω
where Ω represents the boundary and n refers to the surface normal. The constant
boundary voltage implies the tangential electric field equals zero; the electric field will
be perpendicular to the boundary:
If the floating potential touches a point, boundary or domain feature that is not
floating (a Terminal or Ground feature), the floating potential will acquire that feature's
potential. If the floating potential is set to a certain charge, or connected to a circuit,
it behaves like a terminal.
FLOATING POTENTIAL
The Floating potential identifier text area shows the unique identifier for the floating
potential feature. It is used to identify the global variables created by the node, such as
the voltage. The Floating potential group check box on the Settings window for Floating
Potential controls how potentials are assigned to boundary segments in the feature’s
boundary selection. If this check box is not selected, a single potential is applied to all
boundaries in the selection. If the check box is selected, each group of contiguous
boundaries in the selection is given a unique potential. This simplifies the setup of
models with many floating electrodes. The values of the potential at each group of
boundaries are then made available in postprocessing, collected in a vector variable.
The following options are not available if the Floating potential group check box is
selected. Select a specification for the Electric charge — choose User defined to specify
a total Charge Q0 deposited on the surface. Select Circuit to connect the floating
potential to an Electrical Circuit.
Initial Values
The initial values subsection is used to provide initial values for the internal state
variables in the Floating potential. These settings can be useful to provide an initial
guess to the nonlinear solver or the initial value in a time-dependent study.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Distributed Capacitance
The Distributed Capacitance node adds a distributed capacitance boundary condition
according to the following equations for exterior boundaries (left) and interior
boundaries (right):
V ref – V V ref – V
– n ⋅ D = ε 0 ε r -------------------- n ⋅ ( D 1 – D 2 ) = ε 0 ε r --------------------
ds ds
Use this boundary condition to model a thin sheet or film of a dielectric material. The
sheet has the relative permittivity εr and the surface thickness ds, and it is connected
to the reference potential Vref.
DISTRIBUTED CAPACITANCE
Enter the values or expressions for Relative permittivity εr (dimensionless), Surface
thickness ds (SI unit: m), and Reference potential Vref. The default value for the surface
thickness is 10−3 m (1 mm) and 0 V for the reference potential.
Line Charge
For 3D components, use the Line Charge node to specify line charges along the edges
of a geometry. Add a contribution as a Harmonic Perturbation by right-clicking the
parent node or clicking Harmonic Perturbation on the Physics toolbar. For more
information see Harmonic Perturbation — Exclusive and Contributing Nodes in the
COMSOL Multiphysics Reference Manual.
LINE CHARGE
Enter a value or expression to apply a Line charge QL (SI unit: C/m). This source
represents electric charge per unit length and the default is 0 C/m.
Use the Line Charge (Out-of-Plane) node to specify line charges along the points of a
geometry for 2D and 2D axisymmetric components.
POINT SELECTION
Point Charge
The Point Charge node adds a point source to 3D components. The point charge
represents an electric displacement field flowing out of the point.
POINT SELECTION
POINT SELECTION
Use the Change Thickness (Out-of-Plane) node to set the out-of-plane thickness for
specific geometric entities.
ELECTRIC DISPLACEMENT
If the node is used together with an active Piezoelectric Effect multiphysics coupling
node, then these settings are locked. Note that if they are unlocked, then the material
behaves like a dielectric and not a piezoelectric. In this case, the default is to take the
Relative permittivity εrS (dimensionless) values From material. For User defined, select
Isotropic, Diagonal, Symmetric, or Anisotropic and enter values or expressions in the field
or matrix.
The Electrostatic Point Dipole represents the limiting case of zero separation distance
between two equally strong point sources of opposing signs while maintaining the
product between separation distance and source strength at a fixed value (p). The
dipole moment is a vector entity with positive direction from the negative charge to
the positive one.
• For Magnitude and direction enter coordinates for the Electric dipole moment
direction np (dimensionless) and the Electric dipole moment, magnitude p.
• For Electric dipole moment enter coordinates for the Electric dipole moment p.
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.
For User defined enter a value or expression. If another type of temperature dependence
is used other than a linear temperature relation, enter any expression for the
conductivity as a function of temperature.
• 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 Charge
Conservation, Current Conservation, or Ampere’s Law nodes). Then the subnodes are
made available from the context menu (right-click the parent node) as well as from the
Physics toolbar, Attributes menu.
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.
The Volume fractions specified for the materials should be fractional (between 0 and
1) and should add to 1 in normal cases.
The physics interface solves Laplace’s equation for the electric potential using the scalar
electric potential as the dependent variable
– ∇ ⋅ ( ε 0 ε r ∇V ) = 0 .
For a discussion about the boundary element method, see Theory for the
Boundary Elements PDE in the COMSOL Multiphysics Reference
Manual.
When this physics interface is added, these default nodes are also added to the Model
Builder— Charge Conservation, Zero Charge (on exterior boundaries and in 3D also on
interior edges) and Initial Values. Then, from the Physics toolbar, add other nodes that
implement, for example, electric potential and surface charge density conditions. You
can also right-click Electrostatics, Boundary Elements to select physics features from the
context menu.
Selections can also be entered using the Selection List window, available from the
Windows menu on the Home toolbar.
PHYSICS SYMBOLS
Select the Enable physics symbols check box to display symmetry planes (in 3D) and
lines (in 2D) in the Graphics window, as specified in the Symmetry settings.
SYMMETRY
For more information about the Far Field Approximation settings, see
Far-Field Approximation Settings in the COMSOL Multiphysics
Reference Manual.
INFINITY CONDITION
QUADRATURE
To display this section, click the Show button ( ) and select Advanced Physics Options.
MESH CONTROL
The Enable check box is selected by default. This allow the physics interface to control
the meshing process. Information from the physics will be used to automatically set up
an appropriate mesh sequence
DISCRETIZATION
From the Electric potential/Surface charge density list, choose from predefined options
for the boundary element discretization order for the electric potential variable and the
surface charge density variable, respectively. The predefined options represent the
suitable combinations of element orders such as Quadratic/Linear (the default).
For 3D components, the Line charge variable will have the same
discretization order as the Electric potential variable.
The settings under Value types when using splitting of complex variables are important
for sensitivity and optimization computations. See the description of the built-in
operators fsens and fsensimag.
DEPENDENT VARIABLES
The dependent variable (field variable) is for the Electric potential V. The name can be
changed but the names of fields and dependent variables must be unique within a
model.
AVAILABLE NODES
These nodes and subnodes are described in this section:
The Zero Charge node is the default boundary and edge condition available for this
physics interface and described for the Electrostatics interface.
Initial Values
Initial Values adds an initial value for the electric potential that can serve as an initial
condition for a transient simulation or as an initial guess for a nonlinear solver. Add
more Initial Values nodes from the Physics toolbar, Global menu.
INITIAL VALUES
Enter a value or expression for the initial value of the Electric potential V (SI unit: V).
ELECTRIC FIELD
Enter the Relative permittivity εr (dimensionless) to describe the macroscopic
properties of the medium. The default is to take the Relative permittivity values From
material. For User defined, enter a value or expression in the field.
Electric Potential
The Electric Potential node provides an electric potential V0 as the boundary condition
V = V0.
For more information about the Electric Potential condition, see Electric
Potential in the documentation for the Electrostatics interface.
For 3D components, enter the Edge radius (SI unit: m) for the virtual tube
surrounding the edge. The default value is 0.1 m. It is generally assumed
that the edge radius is small in comparison to other geometrical measures
of the model, such as the edge length and the distances between different
edges in the geometry.
For 3D components, enter the Edge radius (SI unit: m) for the virtual tube
surrounding the edge. The default value is 0.1 m. It is generally assumed
that the edge radius is small in comparison to other geometrical measures
of the model, such as the edge length and the distances between different
edges in the geometry.
Ground
The Ground node implements ground (zero potential) as the boundary condition
V = 0.
For more information about the Ground condition, see Ground in the
documentation for the Electrostatics interface.
For 3D components, enter the Edge radius (SI unit: m) for the virtual tube
surrounding the edge. The default value is 0.1 m. It is generally assumed
that the edge radius is small in comparison to other geometrical measures
of the model, such as the edge length and the distances between different
edges in the geometry.
–n ⋅ D = ρs , n ⋅ ( D1 – D2 ) = ρs
For more information about the Surface Charge Density condition, see
Surface Charge Density in the documentation for the Electrostatics
interface.
For 3D components, enter the Edge radius (SI unit: m) for the virtual tube
surrounding the edge. The default value is 0.1 m. It is generally assumed
that the edge radius is small in comparison to other geometrical measures
of the model, such as the edge length and the distances between different
edges in the geometry. Enable Compensate for tube volume to allow for a
potential gradient across the cross section of the virtual tube. This can be
useful for establishing a zero surface charge density condition on the
virtual tube periphery.
Terminal
The Terminal node provides a boundary or domain condition for connection to
external circuits, to transmission lines, or with a specified voltage or charge.
For more information about the Terminal condition, see Terminal in the
documentation for the Electrostatics interface.
For 3D components, enter the Edge radius (SI unit: m) for the virtual tube
surrounding the edge. The default value is 0.1 m. It is generally assumed
that the edge radius is small in comparison to other geometrical measures
of the model, such as the edge length and the distances between different
edges in the geometry.
Zero Charge
The Zero Charge node adds the condition that there is zero charge on the boundary so
that n ⋅ D = 0. This is the default boundary condition on exterior boundaries and on
edges embedded within a Charge Conservation feature.
BOUNDARY SELECTION
Select Manual or All boundaries from the Selection list. Make additional edits to the list
of boundary entity numbers using the Selection toolbar buttons. When All boundaries
is selected from the Selection list, the boundaries exterior to the Electrostatics interface
that intersect the exterior boundaries to the Electrostatics, Boundary Elements
interface are available in the boundary entity number list.
COUPLED INTERFACES
• Select Electrostatics, Boundary Elements as Primary interface and Electrostatics and
Secondary interface.
The physics interface solves a current conservation equation based on Ohm’s law using
the scalar electric potential as the dependent variable.
Current Conservation is the main node, which adds the equation for the electric
potential and provides a Settings window for defining the electrical conductivity as well
as the constitutive relation for the electric displacement field and its associated material
properties, such as the relative permittivity.
When this physics interface is added, these default nodes are also added to the Model
Builder — Current Conservation, Electric Insulation (the default boundary condition),
and Initial Values. Then, from the Physics toolbar, add other nodes that implement, for
example, boundary conditions and current sources. You can also right-click Electric
Currents 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 ec.
PHYSICS-CONTROLLED MESH
Select the Enable check box to allow the physics interface to control the meshing
process. Information from the physics, such as the presence of an infinite elements
THICKNESS (2D)
Enter a default value for the Out-of-plane thickness d (SI unit: m) (see Equation 4-1).
The default value of 1 m is typically not representative for a thin dielectric medium, for
example. Instead it describes a unit thickness that makes the 2D equation identical to
the equation used for 3D components. See also Change Thickness (Out-of-Plane)
(described for the Electrostatics interface).
Select the Activate manual terminal sweep check box to switch on the sweep and invoke
a parametric sweep over the terminals. Enter a Sweep parameter name to assign a
specific name to the variable that controls the terminal number solved for during the
sweep. The default is PortName. The generated lumped parameters are in the form of
capacitance matrix elements. The terminal settings must consistently be of either fixed
voltage or fixed charge type.
The lumped parameters are subject to Touchstone file export. Enter a file path or Browse
for a file. Select a Parameter format (value pairs) for the Touchstone export —
Magnitude and angle (MA) (the default), Magnitude in dB and angle (DB), or Real and
imaginary parts (RI). Select an option from the If file exists list — Overwrite or Create
new. Create new is useful when the model is solved multiple times with different
settings. Select a Parameter to export — Z (the default), Y, or S.
DISCRETIZATION
• Lumped Parameters
• Domain, Boundary, Edge, Point, and Pair Nodes for the Electric
Currents Interface
• Theory for the Electric Currents Interface
Domain, Boundary, Edge, Point, and Pair Nodes for the Electric
Currents Interface
The Electric Currents interface has these domain, boundary, edge, point, and pair
nodes available from the Physics ribbon toolbar (Windows users) or Physics context
menu (Mac or Linux users). You can also right-click to access the context menu (all
users).
n2 ⋅ ( J1 – J2 ) = 0
AVAILABLE NODES
These nodes are available for this physics interface, listed in alphabetical order. Also see
Table 4-1 for a list of interior and exterior boundary conditions, including edge, point,
and pair availability.
1This feature is available with the Piezoresistivity, Domain Currents interface, which
In the COMSOL Multiphysics Reference Manual, see Table 2-3 for links
to common sections and Table 2-4 for common feature nodes. You can
also search for information: press F1 to open the Help window or Ctrl+F1
to open the Documentation window.
Table 4-1 lists the interior and exterior boundary conditions available with this physics
interface. It also includes edge, point, and pair availability.
TABLE 4-2: INTERIOR AND EXTERIOR BOUNDARY CONDITIONS (INCLUDING EDGE, POINT, AND PAIR
AVAILABILITY) FOR THE ELECTRIC CURRENTS INTERFACE
Current Conservation
The Current Conservation node adds the continuity equation for the electrical potential
and provides an interface for defining the electric conductivity as well as the
constitutive relation and the relative permittivity for the displacement current.
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation, and other variables
evaluated in a material reference configuration (material frame). Select Non-solid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
User Defined
For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic depending on the
characteristics of the electrical conductivity, and then enter values or expressions for the
electrical conductivity σ in the field or matrix. The default is 0 S/m. If type of
temperature dependence is used other than a linear temperature relation, enter any
expression for the conductivity as a function of temperature.
Linearized Resistivity
Select Linearized resistivity for a temperature-dependent conductivity (this occurs in,
for example, Joule heating, and is also called resistive heating). The equation
describing the conductivity:
1
σ = ---------------------------------------------------
ρ 0 ( 1 + α ( T – T ref ) )
where ρ0 is the resistivity at the reference temperature Tref, and α is the temperature
coefficient of resistance, which describes how the resistivity varies with temperature.
The default Reference resistivity ρ0 (SI unit: Ω⋅m), Reference temperature Tref
(SI unit: K), and Resistivity temperature coefficient α (SI unit: 1/K) are taken From
material, which means that the values are taken from the domain (or boundary)
material. T is the current temperature, which can be a value that is specified as a model
input or the temperature from a heat transfer interface. The definition of the
temperature field is in the Model Inputs section.
To specify other values for any of these properties, select User defined from the list and
then enter a value or expression for each. The default values are:
ELECTRIC FIELD
Select a Constitutive relation to describe the macroscopic properties of the medium
(relating the electric displacement D with the electric field E) and the applicable
material properties, such as the relative permittivity. For a description of the
constitutive relations Relative permittivity, Polarization, and Remanent electric
displacement, see Electric Field as described for the Charge Conservation node for the
Electrostatics interface. The constitutive relations specific to Electric Currents are:
• Dielectric losses: uses the constitutive relation D = ε0(ε'+ε")E. Specify that the
Relative permittivity (real part) ε' (dimensionless) and the Relative permittivity
(imaginary part) ε" (dimensionless) must be taken From material or be User defined.
For User defined, select Isotropic, Diagonal, Symmetric, or Anisotropic and enter values
or expressions in the field or matrix. The default is 1.
• Loss tangent, loss angle: uses the constitutive relation D = ε0ε'(1+jtanδ)E. Specify the
Relative permittivity (real part) ε' (dimensionless) and Loss angle δ (SI unit: rad).
• Loss tangent, dissipation factor: uses the constitutive relation D = ε0ε'(1+jtanδ)E.
Specify the Relative permittivity (real part) ε' (dimensionless) and the Dissipation
factor tanδ (dimensionless).
Initial Values
The Initial Values node adds an initial value for the electric potential that can serve as
an initial condition for a transient simulation or as an initial guess for a nonlinear solver.
If more than one set of initial values is required, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions and current sources. Add
more Initial Values nodes from the Physics toolbar.
INITIAL VALUES
Enter a value or expression for the initial value of the Electric potential V (SI unit: V).
The default value is 0 V.
The external current density does not contribute to the losses (due to Joule heating),
since there is no electric field associated with it. To include the contribution to the
losses from the external current density, select the Add contribution of the external
current density to the losses check box. Then select an option from the External losses
list—From domain conductivity (the default) or User defined. If From domain
conductivity is selected, the heat source is computed using the conductivity specified in
the material model feature (such as Current Conservation) that is applied in the domain.
For User defined, enter a value for Qe (SI unit: W/m3) to specify a user-defined heat
source.
Current Source
The Current Source node adds a distributed current source Qj in the equation that the
physics interface defines. Use this node with caution as it can violate the current
conservation law that is inherent in Maxwell-Ampère’s law.
CURRENT SOURCE
Enter a value or expression for the Current source Qj (SI unit: A/m3). The default is
0 A/m3.
n⋅J = 0
This boundary condition means that no electric current flows into the boundary. At
interior boundaries, it means that no current can flow through the boundary and that
the electric potential is discontinuous across the boundary. It is also applicable at
symmetric boundaries where the potential is known to be symmetric with respect to
the boundary.
Floating Potential
The Floating Potential node is used when modeling a metallic electrode at floating
potential. It applies a constant voltage V0 on the boundary (for domain features, this
is the boundary enclosing the selected domain), such that the total normal electric
current density J equals a specific current I0:
V = V0
( –n ⋅ J ) dS = I0 ,
∂Ω
where Ω represents the boundary and n refers to the surface normal. The constant
boundary voltage implies the tangential electric field equals zero; the electric field will
be perpendicular to the boundary:
–∇t V = Et = 0
n × E = 0.
In case of I0 = 0 (the default case), the boundary will behave as an unconnected perfect
conductor (a floating equipotential). This is a good approximation when the
conductivity of the electrode is many orders of magnitude larger than that of the
surrounding medium. Although locally the current density may vary, the total current
entering or leaving the boundary equals zero.
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or on the Physics toolbar, click the
Attributes menu and select Harmonic Perturbation. For more information see Harmonic
Perturbation — Exclusive and Contributing Nodes in the COMSOL Multiphysics
Reference Manual.
FLOATING POTENTIAL
The Floating potential identifier text area shows the unique identifier for the floating
potential feature. It is used to identify the global variables created by the node, such as
the voltage. The Floating potential group check box on the Settings window for Floating
Potential controls how potentials are assigned to boundary segments in the feature’s
boundary selection. If this check box is not selected, a single potential is applied to all
boundaries in the selection. If the check box is selected, each group of contiguous
boundaries in the selection is given a unique potential. This simplifies the setup of
models with many floating electrodes. The values of the potential at each group of
boundaries are then made available in postprocessing, collected in a vector variable.
The following options are not available if the Floating potential group check box is
selected. Select a specification for the Floating potential current—choose User defined
(the default) to specify a total Current I0 (SI unit: A) that flows from the electrode. The
default is 0 A, corresponding to an unconnected electrode. Select Circuit to connect
the floating potential to an Electrical Circuit interface.
Initial values
The initial values subsection is used to provide initial values for the internal state
variables in the Floating Potential. These settings can be useful to provide an initial
guess to the nonlinear solver or the initial value in a time-dependent study when the
initial values of the other dependent variables (such as the electric potential) are
computed from a separate physics or study. In most cases, the default values are
appropriate.
Enter an Initial value for voltage Vinit (SI unit: V). The default is 0 V.
n ⋅ ( J1 – J2 ) = Q j
–n ⋅ J = Jn
n ⋅ J = n ⋅ J0
The normal current density is positive when the current flows inward in the domain.
Add a contribution as a Harmonic Perturbation by right-clicking the parent node or
clicking Harmonic Perturbation on the Physics toolbar. For more information see
Harmonic Perturbation — Exclusive and Contributing Nodes in the COMSOL
Multiphysics Reference Manual.
• For Inward current density enter a value or expression for the Normal current density
Jn (SI unit: A/m2). Use a positive value for an inward current flow or a negative
value for an outward current flow. The default is 0 A/m2.
• For Current density enter values or expressions for the components of the Current
density J0 (SI unit: A/m2). The defaults are 0 A/m2.
Distributed Impedance
The Distributed Impedance node adds a distributed impedance boundary condition to
a model.
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or on the Physics toolbar, click the
Attributes menu and select Harmonic Perturbation. For more information see Harmonic
Perturbation — Exclusive and Contributing Nodes in the COMSOL Multiphysics
Reference Manual.
Use this boundary condition to model a thin sheet of a resistive material connected to
a reference potential Vref.
The layer impedance can be specified either with the bulk material conductivity σs, the
relative permittivity εr and layer thickness ds, or directly with the surface resistance ρs
and capacitance Cs. Assuming DC currents, the equation is:
σs
n ⋅ ( J 1 – J 2 ) = ----- ( V – V ref )
ds
1
n ⋅ ( J 1 – J 2 ) = ----- ( V – V ref )
ρs
For the Frequency Domain and Time Dependent study types, this boundary condition
is slightly more sophisticated and accounts also for capacitive coupling. The equations
are:
∂
n ⋅ ( J 1 – J 2 ) = ------- σ ( V – V ref ) + ε 0 ε r ( V – V ref )
1
dL ∂t
DISTRIBUTED IMPEDANCE
Enter the reference potential Vref (SI unit: V). The default is 0 V.
The Terminal node also enables the computation of the lumped parameter of the
system, such as the impedance.
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or on the Physics toolbar, click the
Attributes menu and select Harmonic Perturbation. The Harmonic Perturbation subnode
only has effect if Current or Voltage are selected as terminal type. For more information
see Harmonic Perturbation — Exclusive and Contributing Nodes in the COMSOL
Multiphysics Reference Manual.
TERMINAL
Specify the terminal’s properties. The Terminal name identifies the terminal and the
variables it creates. The name should be numeric for sweeps to work properly.
The Circuit and Power types cannot be used together with a terminal
sweep.
Select:
• Current to specify a total current flowing from the terminal. Enter the current I0 (SI
unit: A). The default is 0 A for an open circuit.
- Enter a Terminal power P0 (SI unit: W) to specify the input power at this terminal.
The default is 1 W.
Be aware that, in general, the values of the voltage and the current might
not uniquely be determined by this constraint and may depend on the
initial value provided for the current and voltage variables. Particularly in
the frequency domain, the absolute phase of the quantities can be
arbitrary. This indeterminacy can have an impact on the solution process.
See Compile Equations in the COMSOL Multiphysics Reference
Manual for a possible solution to this problem.
Initial Values
The initial values subsection is used to provide initial values for the internal state
variables in the Terminal. These settings can be useful to provide an initial guess to the
nonlinear solver or the initial value in a time-dependent study when the initial values
The Initial value for voltage is available for all Terminal types except Voltage. Enter an
Initial value for voltage Vinit (SI unit: V). The default is 0 V.
The Initial value for current is available for Power only. In this case, it provides an initial
guess for the nonlinear solver, and should not be set to zero. The initial value provided
may affect the solution. Enter an Initial value for current Iinit (SI unit: A). The default
is 1 mA.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Lumped Parameters
Electric Shielding
The Electric Shielding node provides an electric shielding boundary condition. Use this
node to model a thin layer of a highly conductive medium that shields the electric field.
The sheet has the electrical conductivity σs and the surface thickness d.
The condition is represented by the following equation for interior boundaries and
(setting J2 = 0) exterior boundaries assuming DC currents
n ⋅ ( J 1 – J 2 ) = – ∇ t ⋅ d ( σ s ∇tV )
For the frequency domain and Time Dependent study types, also displacement
currents are accounted for via the bulk relative permittivity of the sheet; εrs and the
conservation laws change to:
n ⋅ ( J 1 – J 2 ) = – ∇ t ⋅ d ( ( σ s + jωε 0 ε rs ) ∇tV )
∂
n ⋅ ( J 1 – J 2 ) = – ∇ t ⋅ d σ s ∇tV + ε 0 ε rs ∇tV
∂t
Use this boundary condition when approximating a thin domain with a boundary to
reduce the number of mesh elements. The Electric Field and Conduction Current sections
are the same as for Current Conservation.
Contact Impedance
Use the Contact Impedance node on interior boundaries to model a thin layer of
resistive material. It can also be added as a pair using a Pair Contact Impedance node.
The feature allows specifying the contact impedance either by entering the properties
of the material together with the layer thickness, or by entering the impedance
properties of the thin layer directly.
The feature, in stationary studies, applies the following conditions that relate the
normal electric current density with the jump in the electric potential:
σ
n ⋅ J 1 = ------ ( V 1 – V 2 )
ds
σ
n ⋅ J 2 = ------ ( V 2 – V 1 )
ds
1
n ⋅ J 1 = ----- ( V 1 – V 2 )
ρs
1
n ⋅ J 2 = ----- ( V 2 – V 1 )
ρs
The first two equations refer to a layer impedance specified using the bulk material
conductivity σs and the layer thickness ds, while the last two equations refer to the case
in which the surface resistance ρs is specified. The indices 1 and 2 refer to the two sides
of the boundary. These parameters work the same as with Distributed Impedance.
For the frequency domain and Time Dependent study types, this boundary condition
is slightly more sophisticated and accounts also for capacitive coupling. The bulk
material relative permittivity εr or the layer capacitance Cs can be specified.
n ⋅ J 1 = ----- + jωC s ( V 1 – V 2 )
1
ρs
n ⋅ J 2 = ----- + jωC s ( V 2 – V 1 )
1
ρs
∂
n ⋅ J 1 = ----- σ ( V 1 – V 2 ) + ε 0 ε r ( V 1 – V 2 )
1
ds ∂t
n ⋅ J 2 = ----- σ ( V 2 – V 1 ) + ε 0 ε r ∂ ( V 2 – V 1 )
1
ds ∂t
∂
n ⋅ J 1 = ----- ( V 1 – V 2 ) + C s ( V 1 – V 2 )
1
ρs ∂t
∂
n ⋅ J 2 = ----- ( V 2 – V 1 ) + C s ( V 2 – V 1 )
1
ρs ∂t
CONTACT IMPEDANCE
Select a potentially complex-valued Layer specification—Thin layer (the default) or
Surface impedance.
The conductance h is involved in the current flow across the surfaces in contact
according to:
n ⋅ J1 = –hc ( V1 – V2 )
n ⋅ J2 = –hc ( V2 – V1 )
ELECTRICAL CONTACT
Select the Constriction conductance — Cooper-Mikic-Yovanovich correlation (the default),
Mikic elastic correlation, or User defined. For User defined enter a value or expression for
hc (SI unit: S/m2). The default is 0 S/m2.
• Surface roughness, asperities average height σ asp (SI unit: m). The default is 1 μm.
• Surface roughness, asperities average slope masp (dimensionless). The default is 0.4.
• Contact pressure p (SI unit: Pa). The default is 0 Pa.
• For Microhardness enter a value for Hc (SI unit: Pa). The default is 3 GPa.
• For Vickers hardness enter a value for the Vickers correlation coefficient c1 (SI unit:
Pa) and Vickers size index c2 (dimensionless). The defaults are 5 GPa and −0.1,
respectively.
• For Brinell hardness enter a value for HB (SI unit: Pa). The default is 3 GPa. HB
should be between 1.30 and 7.60 GPa.
• For Weighted harmonic mean, enter values or expressions for the Young’s modulus,
upside, Eu (SI unit Pa), the Young’s modulus, downside, Ed (SI unit Pa), the Poisson’s
ratio, upside, νu (dimensionless), and the Poisson’s ratio, downside, νd
(dimensionless). If this node is selected from the Pairs menu, enter instead values or
expressions for the Young’s modulus, source, Esrc (SI unit Pa), the Young’s modulus,
destination, Edst (SI unit Pa), the Poisson’s ratio, source, νsrc (dimensionless), and
the Poisson’s ratio, destination, νdst (dimensionless). The defaults are 0.
• For User defined enter another value or expression for the Contact interface Young’s
modulus Econtact (SI unit: Pa). The default is 1 GPa.
Sector Symmetry
Select Sector Symmetry at interfaces between rotating objects where sector symmetry
is used. It is only available for pairs. A default subnode is added. Right-click to select
additional features from the Fallback Features submenu. In 2D, this feature assumes
rotation around the origin.
PAIR SELECTION
When using nonconforming meshes on the source and destination of a pair, for
numerical stability, a finer mesh should be applied on the destination side for any pair
with a condition that imposes a coupling or a constraint across the pair. The sector
symmetry feature falls into this category.
SECTOR SETTINGS
Enter the Number of sectors (<50) nsect. The default is 2.
Based on space dimension, enter values or expressions in the table for the Axis of
rotation arot.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
POINT SELECTION
POINT SELECTION
DIPOLE SPECIFICATION
Select a Dipole specification—Magnitude and direction (the default) or Dipole moment.
DIPOLE PARAMETERS
Under Dipole Specification:
• For Magnitude and direction enter coordinates for the Electric current dipole moment
direction np and the Electric current dipole moment, magnitude p (SI unit: A·m).
• For Dipole moment enter the components of the Electric current dipole moment p (SI
unit: A·m).
Piezoresistive Material
The Piezoresistive Material is normally used together with a Piezoresistive Effect, Domain
Currents multiphysics coupling node. The node is added by default to the Electric
Currents interface when adding a Piezoresistivity, Domain Currents predefined
multiphysics coupling interface. It is also available from the context menu (right-click
the Electric Currents interface parent node) or from the Physics toolbar.
• Specify a Electrical conductivity, zero stress (SI unit: S/m). This typically comes from
the material added under the Materials node.
• For Piezoresistance form, select a Piezoresistance coupling matrix Πl
(SI unit: m4/(s⋅A2); note that this is equivalent to Ω⋅m/Pa).
• For a Elastoresistance form, select an Elastoresistance coupling matrix Ml
(SI unit: Ω⋅m).
In earlier version this interface was called Electric Currents, Shell. It has
been replaced by The Electric Currents, Layered Shell Interface, that is
found under the AC/DC branch ( ).
The physics interface solves a current conservation equation based on Ohm's law using
the scalar electric potential as the dependent variable.
Current Conservation is the main node, which adds the equation for the electric
potential and provides a Settings window for defining the electrical conductivity as well
as the constitutive relation for the electric displacement field and its associated material
properties such as the relative permittivity.
When this physics interface is added, these default nodes are also added to the Model
Builder— Current Conservation, Electric Insulation (the default edge or point condition),
and Initial Values. Then, from the Physics toolbar, add other nodes that implement, for
Except where described here, the majority of the Settings windows are the
same as for the Electrostatics and Electric Currents interfaces as
referenced. The only difference are:
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 ecs.
OUT-OF-PLANE THICKNESS
When modeling with shells, the Shell Thickness section defines a parameter
that enters the equations everywhere. It is available in all dimensions.
When modeling in 2D, this means that on the physics interface Settings
window, there are two sections—the Shell Thickness and Out-of-Plane
Thickness—available, as well as two feature nodes—Change Shell Thickness
and Change Thickness (Out-of-Plane).
Select the Activate manual terminal sweep check box to switch on the sweep and invoke
a parametric sweep over the terminals. Enter a Sweep parameter name to assign a
specific name to the variable that controls the terminal number solved for during the
sweep. The default is PortName. The generated lumped parameters are in the form of
capacitance matrix elements. The terminal settings must consistently be of either fixed
voltage or fixed charge type.
The lumped parameters are subject to Touchstone file export. Enter a file path or Browse
for a file. Select a Parameter format (value pairs) for the Touchstone export—Magnitude
and angle (MA) (the default), Magnitude in dB and angle (DB), or Real and imaginary
parts (RI). Select an option from the If file exists list—Overwrite or Create new. Select a
Parameter to export—Z (the default), Y, or S.
DEPENDENT VARIABLES
The dependent variable (field variable) is for the Electric potential V. The name can be
changed but the names of fields and dependent variables must be unique within a
model.
• Lumped Parameters
• Boundary, Edge, Point, and Pair Nodes for the Electric Currents,
Single Layer Shell Interface
Boundary, Edge, Point, and Pair Nodes for the Electric Currents,
Single Layer Shell Interface
The Electric Currents, Single Layer Shell Interface has these boundary, edge, point,
and pair nodes 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).
n2 ⋅ ( J1 – J2 ) = 0
AVAILABLE NODES
These nodes and subnodes are described in this section:
Current Conservation
The Current Conservation node adds the continuity equation for the electrical potential
and provides an interface for defining the electric conductivity as well as the
constitutive relation and the relative permittivity for the displacement current.
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation and other variables
evaluated in a material reference configuration (material frame). Select Non-solid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
ELECTRIC FIELD
See Electric Field as described for the Charge Conservation node for the Electrostatics
interface.
User Defined
For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic depending on the
characteristics of the electrical conductivity, and then enter values or expressions in the
field or matrix.
Linearized Resistivity
Select Linearized resistivity for a temperature-dependent conductivity (this occurs in,
for example, Joule heating, and is also called resistive heating). The equation
describing the conductivity:
1
σ = -----------------------------------------------
ρ0 ( 1 + α ( T – T0 ) )
where ρ0 is the resistivity at the reference temperature T0, and α is the temperature
coefficient of resistance, which describes how the resistivity varies with temperature.
The default Reference temperature Tref (SI unit: K), Resistivity temperature coefficient
α (SI unit: 1/K), and Reference resistivity ρ0 (SI unit: Ω⋅m) are taken From material,
which means that the values are taken from the domain material.
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, an Archie’s Law subnode is available.
Initial Values
Initial Values adds an initial value for the electric potential that can serve as an initial
condition for a transient simulation or as an initial guess for a nonlinear solver. Add
more Initial Values nodes from the Physics toolbar.
INITIAL VALUES
Enter a value or expression for the initial value of the Electric potential V (SI unit: V).
The default is 0 V.
Current Source
The Current Source node adds a distributed current source Qj in the equation that the
physics interface defines. Use this node with caution as it might violate the current
conservation law that is inherent in Maxwell-Ampère’s law.
Electric Shielding
The Electric Shielding node, available for 3D components, can be used to model a
geometrically thin section of shell made of a highly conductive medium. The layer
thickness parameter dl specifies the thickness in the direction tangential to the shell,
while the thickness in the direction normal to the shell is taken from the physics
interface or, if present, from a Change Shell Thickness node.
ELECTRIC SHIELDING
The default Relative permittivity εr (dimensionless) and Electrical conductivity σ
(SI unit: S/m) take values From material. For User defined enter different values or
expressions.
Enter a value or expression for the Layer thickness dl (SI unit: m). The default is 1 cm.
The physics interface solves a current conservation equation based on Ohm’s law using
the scalar electric potential as the dependent variable. The electric potential is added
either by the Layered Shell feature or the Homogenized Shell feature. Both features
have advantages and disadvantages, particularly considering shell thickness and
electrical connections between the shells.
For (moderately) thick shells, the layered shell implementation is superior as it provides
a full 3D representation (perpendicular potential gradients and -currents are included
by means of an extra dimension). The homogenized shell implementation instead only
solves for in-plane potential gradients. It requires less degrees of freedom, and is
numerically more stable. It is especially useful for very thin layers, for which the layered
shell implementation may fail. Moreover, it naturally preserves electrical continuity
across edges interior to the shell. Both the Layered Shell feature and the Homogenized
Shell feature are equipped with a number of subfeatures, allowing you to take
advantage of the strength of each, in different parts of the model.
When the physics interface is added, the following default nodes are added to the Model
Builder— Layered Shell and Electric Insulation (the default edge condition). Then, from
the Physics toolbar, other nodes can be added that implement, for example, boundary
• This interface supports the Layered Shell feature, which means each
boundary may represent a number of thin layers with different
electrical properties. Potential variations in the normal direction (in
between and inside layers) are supported. For a general description of
layer and interface selections, see The Layer Selection Section. See also
the section on how to set up of a Layered Material node in the
COMSOL Multiphysics Reference Manual.
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 ecls.
DEPENDENT VARIABLES
The dependent variable (field variable) is for the Electric potential V. The name can be
changed but the names of fields and dependent variables must be unique within a
model. This variable is explicitly used only by Homogenized Shell features. The Layered
Shell features have their own variable name (for the degrees of freedom in the extra
dimension).
DISCRETIZATION
The setting for the discretization order is shared by the Homogenized Shell degrees of
freedom and the Layered Shell degrees of freedom.
The physics features can be of two, fundamentally different, types: Those which act on
the layers themselves, and those which act on the interfaces between layers. The
settings for these two types of features are slightly different.
The default selections in the Layer Selection section differ between different physics
features. This reflects the fact that some physics features (such as the Layered Shell
feature) are more likely to be applied to all layers, whereas other features (such as the
• All layered materials. This option is not available for all features.
• Any Layered Material Link defined under Materials in the current component.
• Any Layered Material Stack defined under Materials in the current component.
• Any Single Layer Material defined under Materials in the current component.
If you select All layered materials, you cannot control individual layers; the contribution
is given to all layers. All information is taken from the definitions made in the layered
material features (Layered Material Link, Layered Material Stack, or Single Layer Material)
under Materials in the current component. This means that a single physics node can
accommodate several different stacking sequences. As geometrical selection, you can
use any object selected in any of the layered material features.
If you select an individual Layered Material Link, Layered Material Stack, or Single Layer
Material, then you can only select geometrical objects which are part of the selections
of that feature. In most cases, you will get access to a list of check boxes where you can
further limit the contributions to individually selected layers or interfaces to which this
contribution is to be added.
LAYER SELECTIONS
When a physics feature acts on layers, the interpretation of layer selections are:
Figure 4-3: Boundary selection with interfaces ‘layer 2 - layer 3’ and ‘layer 5 up’ selected.
• The Homogenized Shell feature imposes electrical continuity among all adjacent
homogenized shells, regardless whether they belong to the same Layered Material.
• The Layered Shell feature exclusively imposes electrical continuity within layered
shells sharing the same Layered Material.
• The Layered Shell feature and the Homogenized Shell feature are electrically insulated
by default (even when sharing the same Layered Material). To override this insulation
you can use the Homogenized Shell Connection feature.
For more information; these cases are demonstrated in the Solid to Multilayered Shell
Comparison model available in the Application Library.
The extra dimension coordinate has a name like x_llmat1_xdim. The middle part of
the coordinate name is derived from the tag of the layered material definition where it
is created; in this example a Layered Material Link.
Boundary and Edge nodes for the Electric Currents, Layered Shell
Interface
The The Electric Currents, Layered Shell Interface has these boundary and edge nodes
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).
AVAILABLE NODES
These nodes and subnodes are described in this section:
Layered Shell
The Layered Shell node adds the continuity equation for the electric potential and
allows for defining the electric conductivity of the layers. This interface closely relates
to Current Conservation under The Electric Currents Interface. Even though the
interface works on the boundary level, the resulting equations correspond directly to
a volumetric model. The Layered Shell includes an extra dimension to describe
perpendicular gradients.
Notice that the Layered Shell requires a Layered Material. This has consequences for the
continuity equation. For more information about continuity handling, see Electrical
Continuity in the Electric Currents, Layered Shell Interface.
LAYER SELECTION
Select the applicable layers (the default setting is All layered materials). Additionally,
there is a shorthand available for creating a Single Layer Material (the plus, next to the
Layer Selection setting).
User Defined
For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic depending on the
characteristics of the electrical conductivity, and then enter values or expressions in the
field or matrix.
Interlayer Connection
The Interlayer Connection is available both for boundaries, and edges. In both cases it
removes for the selected layers, the potential gradient in the perpendicular direction
(in the direction of the extra dimension). The layers will become electrically
connected, and the top and bottom of all involved layers will have the same potential.
Multiple Interlayer Connection features can be applied to create groups of connected
layers.
An Interlayer Connection on a single layer can be used to locally convert the full 3D
formulation into an in-plane conservation problem. For thin highly conductive layers,
this is numerically more stable.
LAYER SELECTION
Select the applicable layers. The layered material is given by the Layered Shell parent
feature.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
The Homogenized Shell Connection removes for the selected layers, the potential
gradient in the direction normal to the shell (in the direction of the extra dimension),
by setting the potential equal to the one in the base selection. The base selection
potential is the one used by the Homogenized Shell feature. Setting the potential equal
to the base potential creates an electrical connection between the selected layers, and
Homogenized Shell boundaries adjacent to the selected edge.
LAYER SELECTION
Select the applicable layers. The layered material is given by the Layered Shell parent
feature.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Homogenized Shell
The Homogenized Shell feature takes the material properties from one or more selected
layers, and determines their effective material properties. It integrates the conductivity
in the direction normal to the shell (in the direction of the extra dimension), taking
the local layer thickness into account:
1
σ eff = ---
d σ dx ,
XD
where d refers to the total thickness of all layers involved. For a discrete set of layers
having their own uniform material properties, this integral effectively translates to a
sum of conductivities, multiplied by their corresponding local layer thickness:
This represents a parallel connection of resistors. The Varying Thickness subfeature can
be used to locally affect the degree to which a particular layer contributes to the sum.
The resulting material properties are then used in a current conservation problem, that
only considers the potential gradients tangential to the shell (in-plane conservation
law):
∇ T ⋅ ( dσ eff E T ) = ∇ T ⋅ J s = Qs dx ,
XD
where the product dσeffET is effectively the surface current density Js in A/m. ET is
the tangential electric field, and Qs represents a volumetric source of current. Without
sources, the right-hand side becomes zero, and in-plane current should be conserved.
When applied to a single layer, the Homogenized Shell is virtually identical to the
traditional, non-layered Current Conservation feature in the Electric Currents, Single
Layer Shell Interface.
LAYER SELECTION
Select the layers to be included in the integral for the effective material properties.
Varying Thickness
The Varying Thickness feature is a subfeature to the Homogenized Shell feature. It allows
to set a spatially or temporally varying thickness. This is typically used for modeling
material deposition and such. For more information on varying thickness and the
homogenized shell, see the Homogenized Shell feature.
LAYER SELECTION
Select the layers for which the thickness should be varied.
VARYING THICKNESS
Enter a value or expression for the Selected layers thickness dl (SI unit: m). May depend
on location, time, temperature, or other quantities. The default value is 0.1 mm.
In case of an edge, it represents an influx of current Qj,s from the edge into the
neighboring layers (SI unit: A/m2). In the 3D representation of the layer stack, you
could consider it a “boundary current source” applied to one of the boundaries on the
outer perimeter of the layer. For more information on the 3D representation of the
layer stack, see the The Layer Selection Section.
Use this node with caution as it might violate the current conservation law
that is inherent in Maxwell-Ampère’s law.
In the 3D representation of the layer stack, you could consider it a “boundary current
source” applied to either the top or bottom of the layer stack, or on an interface
in-between layers. For more information on the 3D representation of the layer stack,
see the The Layer Selection Section.
Use this node with caution as it might violate the current conservation law
that is inherent in Maxwell-Ampère’s law.
INTERFACE SELECTION
Select the interfaces for which the source contribution should be added.
If you do not know whether to use the Electric Currents or the Electrostatics interface,
which both solve for the scalar electric potential V, consider using an explicit charge
transport model. See Charge Relaxation Theory.
In this section:
• Electrostatics Equations
• The Electrostatics Interface in Time Dependent or Frequency Domain Studies.
• Effective Material Properties in Effective Media and Mixtures
• Effective Conductivity in Effective Media and Mixtures
• Effective Relative Permittivity in Effective Media and Mixtures
• Effective Relative Permeability in Effective Media and Mixtures
• Archie’s Law Theory
For a discussion about the boundary element method used with The
Electrostatics, Boundary Elements Interface, see Theory for the
Boundary Elements PDE in the COMSOL Multiphysics Reference
Manual.
J c = σE
∂-----
ρ-
+ ∇ ⋅ Jc = 0
∂t
∇ ⋅ ( εE ) = ρ
By combining these, one can deduce the following differential equation for the space
charge density in a homogeneous medium
∂-----
ρ- σ
+ --- ρ = 0
∂t ε
–t ⁄ τ
ρ ( t ) = ρ0 e
where
ε
τ = ---
σ
is called the charge relaxation time. For a good conductor like copper, τ is of the order
of 10−19 s, whereas for a good insulator like silica glass, it is of the order of 103 s. For
a pure insulator, it becomes infinite.
When modeling real-world devices, there is not only the intrinsic time scale of the
charge relaxation time but also an external time scale t at which a device is energized
or the observation time. It is the relation between the external time scale and the
charge relaxation time that determines what physics interface and study type to use.
The results are summarized in Table 4-3 below,
TABLE 4-3: SUITABLE PHYSICS INTERFACE AND STUDY TYPE FOR DIFFERENT TIME-SCALE REGIMES.
By combining the definition of the potential with Gauss’ law, you can derive the
classical Poisson’s equation. Under static conditions, the electric potential V is defined
by the equivalence E = −∇V. Using this together with the constitutive relation D = ε0E
+ P between D and E, you can rewrite Gauss’ law as a variant of Poisson’s equation
– ∇ ⋅ ( ε 0 ∇V – P ) = ρ
This equation is used in the Electrostatics interface. It is worth noting that Gauss’ law
does not require the charge distribution to be static. Thus, provided dynamics are slow
enough that induced electric fields can be neglected and hence a scalar electric
potential is justified, the formulation can be used also in the Time Dependent study
type. That typically involves either prescribing the charge dynamics or coupling a
separate formulation for this.
J c = σE + J e
where Je is an externally generated current density. The static form of the equation of
continuity then reads
∇ ⋅ J c = – ∇ ⋅ ( σ ∇V – J e ) = 0
– ∇ ⋅ ( σ ∇V – J e ) = Q j
This equation is used in the static study type for the Electric Currents interface.
∇ ⋅ Jc = ∇ ⋅ ( σE + J e ) = –jωρ
with the equation (∇⋅ D = ρ) yields the following equation for the frequency domain
study type:
For the Time Dependent study type, use the transient equation of continuity
∂ρ
∇ ⋅ J c = ∇ ⋅ ( σE + J e ) = –
∂t
– ∇ ⋅ ∂ ( ε 0 ∇V + P ) – ∇ ⋅ ( σ∇V – J e ) = 0
∂t
These dynamic formulations are valid as long as induced electric fields can be ignored
and hence the electric field is essentially curl free. This condition is fulfilled provided
that skin effect and wave propagation effects can be ignored. The skin depth must be
much larger than the geometrical dimensions of the modeled device and so must the
wavelength. Note also that these formulations can be used to model dielectric regions
of capacitive/resistive devices even though the interior of electrodes sometimes does
not meet the large skin depth condition. In that case, the electrodes must only be
represented as boundary conditions (fixed or floating potential). The interior, metallic
domains are not included in the analysis. Obviously, this is only a valid approach for
devices where metallic electrodes do not entirely bypass (short circuit) the
capacitive/resistive layers. If metallic electrodes short circuit the capacitive/resistive
layers, the time evolution of the current is determined by inductive and resistive effects
Electrostatics Equations
Under static conditions, the electric potential, V, is defined by the relationship:
E = – ∇V
Combining this equation with the constitutive relationship D = ε0E + P between the
electric displacement D and the electric field E, it is possible to represent Gauss’ law
as the following equation:
– ∇ ⋅ ( ε 0 ∇V – P ) = ρ
In this equation, the physical constant, ε0 (SI unit: F/m) is the permittivity of vacuum,
P (SI unit: C/m2) is the electric polarization vector, and ρ (SI unit: C/m3) is a space
charge density. This equation describes the electrostatic field in dielectric materials.
For in-plane 2D modeling, the Electrostatics interface assumes a symmetry where the
electric potential varies only in the x and y directions and is constant in the z direction.
This implies that the electric field, E, is tangential to the xy-plane. With this symmetry,
the same equation is solved as in the 3D case. The physics interface solves the following
equation where d is the thickness in the z direction:
– ∇ ⋅ d ( ε 0 ∇V – P ) = ρ
The axisymmetric version of the physics interface considers the situation where the
fields and geometry are axially symmetric. In this case, the electric potential is constant
in the ϕ direction, which implies that the electric field is tangential to the rz-plane.
Users that have a license supporting the piezoelectric material model can
also include conduction losses in the frequency domain.
The support for dynamic studies simplifies the coupling of the Electrostatics interface
with other physics interfaces. Using the physics interface in a dynamic study is a valid
approximation only if the time-scale (or the frequency) of the study is so slow that
transient electromagnetic effects can be neglected; for example, in acoustic or
structural problems.
Most of the features in the Electrostatics interface have the same behavior in all study
types, with the important exception of the Terminal node.
The Electrostatics interface also supports the small-signal analysis study sequence,
which can be used when a time-harmonic perturbation is superposed on a static bias
charge or voltage.
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
1 θi θ1 θ2
--- =
σ ----σ-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
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.
In this section:
J = σE + J e
where σ is the electrical conductivity (SI unit: S/m), and Je is an externally generated
current density (SI unit: A/m2). The static form of the equation of continuity then
states:
∇ ⋅ J = – ∇ ⋅ ( σ ∇V – J e ) = 0
– ∇ ⋅ ( σ ∇V – J e ) = Q j
In planar 2D the Electric Currents interface assumes that the model has a symmetry
where the electric potential varies only in the x and y directions and is constant in the
z direction. This implies that the electric field, E, is tangential to the xy-plane. The
– ∇ ⋅ d ( σ ∇V – J e ) = dQ j (4-1)
In 2D axisymmetry, the Electric Currents interface considers the situation where the
fields and geometry are axially symmetric. In this case, the electric potential is constant
in the ϕ direction, which implies that the electric field is tangential to the rz-plane.
∇⋅J = ∇ ⋅ ( σE + Je ) = –jωρ
∇⋅D = ρ
∂ρ
∇ ⋅ J = ∇ ⋅ ( σE + J e ) = –
∂t
– ∇ ⋅ ∂ ( ε 0 ∇V – P ) – ∇ ⋅ ( σ∇V – J e ) = Q j (4-2)
∂t
In planar 2D the dynamic formulations also involves the thickness d in the z direction:
∂
– ∇ ⋅ d (( ε 0 ∇V – P ) – ∇ ⋅ d ( σ∇V – J e )) = dQ j
∂t
The current densities at the upside and downside boundaries depend on the potential
difference according to the relations:
n ⋅ J1 = –hc ( V1 – V2 )
n ⋅ J2 = –hc ( V2 – V1 )
masp,u asp,u
Y
SURFACE ASPERITIES
The microscopic surface asperities are characterized by the average height σu, asp and
σd, asp and the average slope mu, asp and md, asp. The RMS values σasp and masp are
(4.16 in Ref. 1):
σ asp = σ u2, asp + σ d2, asp m asp = m u2, asp + m d2, asp
CONSTRICTION CONDUCTANCE
Here, Hc is the microhardness of the softer material, p is the contact pressure, and
σcontact is the harmonic mean of the contacting surface conductivities:
2σ u σ d
σ contact = -------------------
σu + σd
-----------------------------------
1
p p ( 1 + 0.071c2 )
------- = ------------------------------------------------------
Hc σ asp c2
c 1 1.62 ----------- - m asp
σ 0
The coefficients c1 and c2 are the Vickers correlation coefficient and size index,
respectively, and σ0 is equal to 1 µm. For materials with Brinell hardness between 1.30
and 7.60 GPa, c1 and c2 are given by the correlation below (4.16.1 in Ref. 1):
c1 HB HB 2 HB 3
------- = 4.0 – 5.77 -------- + 4.0 -------- – 0.61 --------
H0 H0 H0 H0
HB
c 2 = – 0.37 + 0.442 --------
c1
m asp 0.94
h c = 1.54σ contact ------------- ------------------------
2p
σ asp mE contact
1 - 1 – ν u2 1 – ν d2
------------------ = --------------- + ---------------
E contact Eu Ed
where Eu and Ed are the Young’s moduli of the two contacting surfaces and νu and νd
are the Poisson’s ratios.
Electrical Contact
REFERENCE
1. M.M. Yovanovich and E.E. Marotta, “Thermal Spreading and Contact Resistance,”
Heat Transfer Handbook, A. Bejan and A.D. Kraus, John Wiley & Sons, 2003.
In this section:
– ∇ t ⋅ d ( σ∇ t V – J e ) = dQ j
– ∇ t ⋅ d ( ( σ + jωε 0 )∇ t V – ( J e + jωP ) ) = dQ j
∂
– ∇ t ⋅ d (( ε 0 ∇ t V – P ) – ∇ t ⋅ d ( σ∇ t V – J e )) = dQ j
∂t
THEORY FOR THE ELECTRIC CURRENTS, SINGLE LAYER SHELL INTERFACE | 211
212 | CHAPTER 4: ELECTRIC FIELD INTERFACES
5
This chapter summarizes the functionality of the magnetic field interfaces found
under the AC/DC branch ( ).
In this chapter:
213
The Magnetic Fields Interface
The Magnetic Fields (mf) interface ( ), found under the AC/DC branch ( ) when
adding a physics interface, is used to compute magnetic field and induced current
distributions in and around coils, conductors, and magnets. Depending on the licensed
products, stationary, frequency-domain, small-signal analysis, and time-domain
modeling are supported in 2D and 3D. Note that the frequency and time domain
formulations become ill-posed when approaching the static limit. One may extend the
useful frequency range downward by adding a low conductivity.
The physics interface solves Maxwell's equations, which are formulated using the
magnetic vector potential and, optionally for coils, the scalar electric potential as the
dependent variables.
The main node is Ampère’s Law, which adds the equation for the magnetic vector
potential and provides an interface for defining the constitutive relations and its
associated properties, such as the relative permeability.
When this physics interface is added, these default nodes are also added to the Model
Builder — Magnetic Fields, Ampère’s Law, Magnetic Insulation (the default boundary
condition), and Initial Values. Then, from the Physics toolbar, add other nodes that
implement boundary conditions and external currents. You can also right-click
Magnetic Fields 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 mf.
BACKGROUND FIELD
This section allows the specification of a background magnetic vector potential (that
generates a background magnetic flux density).
For Magnetic vector potential, enter values or expressions for the Background magnetic
vector potential Ab (SI unit: Wb/m). The defaults are 0 Wb/m. For Uniform magnetic
flux density, in 2D or 3D components, enter the values or expressions of the
components of the Uniform magnetic flux density Bb (SI unit: T); in 2D axisymmetric
components, enter the value or expression of the Uniform axial magnetic flux density
Bb,z (SI unit: T). The defaults are 0 T. The specified background magnetic flux
densities must be uniform (space-independent), but they can be a function of time.
When a background field is active, this solves for the relative (perturbated) field only.
This setting, together with the External Magnetic Vector Potential feature, is useful to
introduce an external field generated by systems not included in the model.
COMPONENTS
This section is only available in 2D and 2D axially symmetric components.
The current vector has the same direction as the magnetic vector potential. This setting
also controls the direction in which applied and induced currents can flow in the
model. The default option is to solve for the out-of-plane component only for 2D and
2D axisymmetric components.
From the practical viewpoint this choice is equivalent to deciding in what directions
the electric current is allowed to flow (out-of-plane currents, in-plane currents, or
currents flowing in all three coordinate directions) and affects other settings in the
model, for example, the Port Properties>Type of port section for the Lumped Port
node.
Use the Change Thickness (Out-of-Plane) node (described for the Electrostatics
interface) to define specific geometric entities (for example, domains) instead of a
global setting for the thickness.
ERROR CHECK
To display this section, click the Show button ( ) and select Advanced Physics Options.
When the Check applicability of features in study check box is selected, any features that
are incompatible with the study will generate an error message when trying to solve or
show the default solver. No solver will be generated. Deselect it and you will be able
to run the model, possibly with runtime errors instead. It is available to allow the
advanced user to tweak any feature and use it outside of its intended study scope.
PHYSICS-CONTROLLED MESH
Select the Enable check box to allow the physics interface to control the meshing
process. Information from the physics, such as the presence of an infinite elements
DEPENDENT VARIABLES
The dependent variable is the Magnetic vector potential A. You can change both its field
name and the individual component variable names. If the new field name coincides
with the name of another magnetic vector potential field in the model, the physics
interfaces share degrees of freedom and component names. The new field name must
not coincide with the name of a field of another type, or with a component name
belonging to some other field. Component names must be unique within a model,
except for fields of the same type sharing a common field name.
DISCRETIZATION
• Domain, Boundary, Point, and Pair Nodes for the Magnetic Fields
Interface
• Lumped Parameters
• Theory for the Magnetic Fields Interface
Domain, Boundary, Point, and Pair Nodes for the Magnetic Fields
Interface
The Magnetic Fields interface has these domain, boundary, point, and pair nodes
available, which are listed in alphabetical order.
n2 × ( A1 – A2 ) = 0
n2 × ( H1 – H2 ) = 0
Table 5-1 lists the interior and exterior boundary conditions available with this physics
interface.
TABLE 5-1: INTERIOR AND EXTERIOR BOUNDARY CONDITIONS FOR THE MAGNETIC FIELDS INTERFACE
Available Nodes
These nodes, listed in alphabetical order, are available from the Physics ribbon toolbar
(Windows users), Physics context menu (Mac or Linux users), or by right-clicking to
access the context menu (all users).
In the COMSOL Multiphysics Reference Manual, see Table 2-3 for links
to common sections and Table 2-4 for common feature nodes. You can
also search for information: press F1 to open the Help window or Ctrl+F1
to open the Documentation window.
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation, and other variables
evaluated in a material reference configuration (material frame). Select Non-solid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
CONDUCTION CURRENT
This section is described for the Current Conservation feature.
The options Effective medium and Archie’s law require additional subnodes. If Effective
medium is selected, a Effective Medium subnode is available from the context menu
(right-click the parent node) as well as from the Physics toolbar, Attributes menu. If
Archie’s law is selected, add an Archie’s Law subnode in the same way. These subnodes
contain additional settings to specify how the material properties are computed.
Effective medium models a mixture of materials whose properties are computed by
averaging the properties of the components. Archie’s law models a conductive liquid in
a nonconductive matrix.
ELECTRIC FIELD
The default Relative permittivity εr (dimensionless) for the media is used From material
and defined on the shell domain. For User defined, select Isotropic, Diagonal, Symmetric,
or Anisotropic based on the characteristics of the permittivity and then enter values or
expressions in the field or matrix. If Effective medium is selected, a Effective Medium
subnode is available from the context menu (right-click the parent node) as well as
from the Physics toolbar, Attributes menu, which can specify the relative permittivity of
the mixture.
The equation for the selected constitutive relation displays under the list. For all
options, the default uses values From material, or select User defined to enter a different
value or expression.
Select a Constitutive relation — Relative permeability (the default), H-B curve, Magnetic
losses, Remanent flux density, Magnetization, Effective H-B curve, or Hysteresis
Jiles-Atherton model.
Relative Permeability
Select Relative permeability μr (dimensionless) to use the constitutive relation
B = μ0μrH. For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic and
enter values or expressions in the field or matrix. If Effective medium is selected, a
Effective Medium subnode is available from the context menu (right-click the parent
node) as well as from the Physics toolbar, Attributes menu, which can specify the
relative permeability of the mixture.
H-B Curve
Select H-B curve |H| (SI unit: A/m) to use a curve that relates magnetic flux density B
and the magnetic field H as |H| = f(|B|).
The Magnetic field norm setting can take the values From material, External, or User
defined.
When External is selected, specify the External material to use (from the Materials node
under Global Definitions). This setting allows using material models or constitutive
relations defined in an external library. See Working with External Materials for more
information.
When User defined is selected, specify a user-defined expression for the magnetic field
norm. The direction of the magnetic field is taken to be the same as the direction of
the magnetic flux density at each point.
• The default relative permeability μr (dimensionless) uses values From material. For
User defined, select Isotropic, Diagonal, Symmetric, or Anisotropic based on the
characteristics of the relative permeability and enter another value or expression in
the field or matrix.
• Enter x and y components for the Remanent flux density Br. For 3D components,
enter x, y, and z components.
Magnetization
Select Magnetization M (SI unit: A/m) to use the constitutive relation B = μ0H + μ0M.
Enter x and y components. For 3D components, enter x, y, and z components.
The Direction of magnetization is the only input that normally should be entered in the
physics.
For a more through description of the model, its parameters and literature
references, see The Jiles-Atherton Hysteresis Model in the theory section.
CONDUCTION CURRENT
This section is described for the Current Conservation feature.
ELECTRIC FIELD
The default Relative permittivity εr (dimensionless) for the media is used From material
and defined on the shell domain. For User defined, select Isotropic, Diagonal, Symmetric,
or Anisotropic based on the characteristics of the permittivity and then enter values or
expressions in the field or matrix.
MAGNETIC FIELD
Specify the constitutive relation that describes the macroscopic properties of the
medium (relating the magnetic flux density B and the magnetic field H) for a
magnetostrictive material:
B = μ 0 [ H + M ( H, S mech ) + M r ]
• Magnetostriction
• Modeling Magnetostrictive Materials
• Magnetostrictive Material
Initial Values
The Initial Values node adds an initial value for the magnetic vector potential A that can
serve as an initial value for a transient simulation or as an initial guess for a nonlinear
solver.
INITIAL VALUES
Enter values or expressions for the Magnetic vector potential A (SI unit: Wb/m). The
defaults are 0 Wb/m.
The external current density does not contribute to the losses (due to Joule heating),
since there is no electric field associated with it. To include the contribution to the
losses from the external current density, select the Add contribution of the external
current density to the losses check box. Then select an option from the External losses
list—From domain conductivity (the default) or User defined. If From domain
conductivity is selected, the heat source is computed using the conductivity specified in
The Velocity (Lorentz term) node adds velocity v. The external current is equal to
σv × B. It also only valid when solving for both the electric potential and the magnetic
vector potential using The Magnetic and Electric Fields Interface.
An operational definition of when it can be used is that the moving domain should only
contain an induced magnetic source (magnetization plus eddy currents) that has to be
stationary with respect to the motion. Thus, it cannot be used for modeling projectiles
of finite length or projectiles containing magnets. It can be used to model conductive,
homogeneous spinning disks (magnetic brakes); magnets over a moving infinite
homogeneous plane (maglev trains); and flow of homogeneous conducting fluid past
a magnet (liquid metal pumps or Hall generators/thrusters, for example).
If you are not sure how to proceed, contact the COMSOL Support
Center: http://www.comsol.com/support.
Magnetic insulation is a special case of the magnetic potential boundary condition that
sets the tangential component of the magnetic potential to zero.
This node is used for the modeling of a lossless metallic surface, for example, a ground
plane or as a symmetry type boundary condition. The node imposes symmetry for
magnetic fields and “magnetic currents.” In the transient and time harmonic
formulations, it also imposes antisymmetry for electric fields and electric currents. The
node 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 are automatically balanced by induced surface currents.
The Magnetic Insulation node can also be applied on interior boundaries. The
boundary will then support two surface current densities on the two sides, denoted bu
Jsu (upside) and Jsd (downside).
Js
J
I'
I
Js
Figure 5-1: The magnetic insulation boundary condition is used on exterior and interior
boundaries representing the surface of a lossless metallic conductor or (on exterior
boundaries) 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 vector
potential (and electric field) vanishes at the boundary.
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
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or on the Physics toolbar by clicking the
Attributes menu and selecting Harmonic Perturbation. For more information see
Harmonic Perturbation — Exclusive and Contributing Nodes in the COMSOL
Multiphysics Reference Manual.
MAGNETIC FIELD
Enter the value or expression for the Magnetic Field H0 (SI unit: A/m) vector
coordinates. The defaults are 0 A/m.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
–n × H = Js
n × ( H1 – H2 ) = Js
SURFACE CURRENT
Enter values or expressions for the Surface current density Js0 (SI unit: A/m)
coordinates. The defaults are 0 A/m.
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.
Magnetic Potential
The Magnetic Potential node adds a boundary condition for the magnetic vector
potential:
n × A = n × A0
MAGNETIC POTENTIAL
Enter a value or expression for the Magnetic vector potential A0 (SI unit: Wb/m)
coordinates.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Js=0
I'
I
J=0
Figure 5-2: 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.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
DIPOLE SPECIFICATION
Select a Dipole specification — Magnitude and direction (the default) or Dipole moment.
DIPOLE PARAMETERS
Under Dipole specification:
• For Magnitude and direction enter coordinates for the Electric current dipole moment
direction np (the defaults are 0) and Electric current dipole moment, magnitude p
(SI unit: A⋅m) (the default is 0 A⋅m).
• For Dipole moment enter coordinates for the Electric current dipole moment p
(SI unit: A⋅m). The defaults are 0 A⋅m.
DOMAIN SELECTION
From the Selection list, choose the domains to define the gauge-fixing potential ψ. In
most cases, the feature should be applied to all domains where the magnetic vector
potential A is solved for. By default, the selection is set to All domains, ensuring that
the gauge fixing is applied to all the valid domains in the model.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
This section allows a more fine control on the boundary conditions for ψ applied by
the Gauge Fixing feature. The domain equation for ψ only imposes a condition on the
gradient, so it is important to constrain the absolute value of ψ to ensure a non-singular
model. The Constant value on insulation boundaries check box (selected by default)
imposes a constant value on the conductive boundaries in the model, such as Magnetic
Insulation. Select the Method to enforce this condition: Constrain value (the default) or
Constrain tangential gradient.
If you apply Gauge Fixing across pair boundaries, the internal logic for
constraining ψ in at least one point does not work. The problem is most
likely to appear in gauge fixing for rotating machinery. Assume you have
one gauge fixed Ampère's Law domain in the stator and one in the rotor,
then one must use separate Gauge Fixing features for these as the geometry
analysis for setting up point constraints on ψ does not work across pairs.
It can also be handled by adding manual constraints on ψ.
Coil
The Coil node can be used to model coils, cables and other conductors subject to a
lumped excitation, such as an externally applied current or voltage. The Coil feature
transforms this lumped excitation into local quantities (electric field and electric
current density), and computes lumped parameters of interest such as impedance, and
inductance.
• Single conductor, which models a conductive body such as a wire, busbar, or other
metallic conductor in which the current flows freely due to the material’s
conductivity. This model can be used when the current flow has a well-defined
beginning and end (for example, connections to an external source) or is closed in
a loop.
The Domain Selection has to be complete in the sense that selecting only
part of a contiguous conductor will lead to unphysical results.
The Coil feature is available both for domain and for boundary selections. In the latter
case, it represents a flat coil or a conductor with a thickness negligible compared with
the other dimensions. Different subnodes can be added to the Coil node in different
cases.
The global Harmonic Perturbation subnode is available from the context menu
(right-click the parent node and select it from the Global menu) or from the Physics
toolbar, Attributes menu. The subnode can be used to apply a harmonic perturbation
to the coil excitation.
In 2D and 2D axisymmetric components, the Coil feature supports the Coil group
functionality, that can be activated by selecting the corresponding check box.
The Coil group option assumes that the selected domains represent cross
sections of the same conductor going in and out of the modeling plane.
These are expected to have the same areas. The same total current will be
imposed in each domain, even if the domain areas are not equal but the
computed concatenated flux, coil voltage and inductance will be
incorrect. For cases with varying cross section areas, it is recommended to
use separate coil features that are coupled using The Electrical Circuit
Interface.
Refer to the Coil Features section in the modeling guide for more information about
this node.
See Coil Features in this guide to learn more about using this feature.
COIL
Coil Name
Enter a Coil name. This name is appended to the global variables (current, voltage)
defined by this coil, and it can be used to identify the coil in a Coil Geometry Analysis
study step.
Conductor Model
Select the Conductor model for the coil. The choices correspond to rather different
physical model, although the set up is similar. The Single conductor model (the default)
is appropriate for solid, massive current-carrying conductors. The Homogenized
multi-turn model represent a bundle of tiny wires that are not geometrically resolved
but taken into account in their average effect. The choice of Conductor model affects
the controls that are visible in the GUI and the available subnodes for the Coil feature.
Select a Coil Type—Linear (the default), Circular, Numeric, or User defined. The different
alternatives are described in the following sections. Also see Using Coils in 3D Models
for more information.
To respect the current conservation law, the applied currents cannot originate from
interior boundaries. A Linear coil should therefore be terminated on exterior
boundaries.
Coil Excitation
Select a Coil excitation—Current (the default), Voltage, Circuit (voltage), Circuit
(current), or Power (2D and 2D axisymmetric components only).
• Current forces a total current flowing in the coil wire. Enter a Coil current Icoil (SI
unit: A). The default is 1 A. See the box below for study limitations on this setting.
• Voltage applies a total voltage across the coil enter a Coil voltage Vcoil (SI unit: V).
The default is 1 V.
• Circuit (current) works similarly to the Current excitation, but in this case the inputs
are provided by a circuit connection.
• Circuit (voltage) works similarly to the Voltage excitation, but in this case the inputs
are provided by a circuit connection.
• Power (only available for 2D and 2D axisymmetric components) forces the coil input
power (cycle-average in frequency studies) to the specified value. Choosing this
option makes the problem nonlinear. For Power enter a Coil power Pcoil (SI unit: W).
The default value is 1 W.
When using the Current, Circuit (current) and Power options, the coil
feature sets up a control problem for the coil voltage and current. Due to
its complexity, the following limitations apply:
CONDUCTION CURRENT
This section is available only when Single conductor is selected as the Coil model. In this
case, the coil represents a solid, massive conductor and the conductivity of the material
is required to compute the current density flowing in it.
Number of Turns
Enter the Number of turns N. The default is 10. This is the number of tiny wires
constituting the coil. The coil resistance is affected by this number and so is the current
density in the coil as it together with the Current setting defines the number of
Ampère-turns in the coil.
• For User defined, enter the value of the cross section area acoil (SI unit: m2). The
default is 10−6 m2.
• For Standard wire gauge, enter the SWG size. Sizes between 7/0 and 50 are
available. The default size is 0.
• For American wire gauge (Brown & Sharpe), enter the AWG size. Sizes between 0000
and 40 are available. Sizes such as 0000 can be also written as 4/0. The default size
is 0.
• For From round wire diameter, enter the diameter of the individual wire dcoil (SI
unit: m). The cross-section area of the round wire will be computed from it. The
default value of dcoil is 1 mm.
STABILIZATION
To display this section, click the Show button ( ) and select Stabilization. This section
is available only in 3D components when using Homogenized multi-turn as the
Conductor model and it contains advanced settings relative to the accuracy and
stabilization of the solution.
The Accurate coil voltage calculation check box enables a current filtering functionality
that improves the accuracy of the computed electric field and the induced coil voltage,
at the cost of a slightly increased number of degrees of freedom. This functionality is
only applicable for frequency domain studies, and is active by default.
For the purpose of stabilizing the solution, the coil feature can apply a small electric
conductivity to the coil domain. Use the Stabilization combo box to specify the value
of the conductivity. Choose Automatic (the default) to use a conductivity automatically
computed by the coil. In frequency domain studies, the conductivity is chosen so that
the skin depth in the coil is 1000 times larger than the coil length (see the sections Coil
Geometry and User Defined Coil Geometry below). It is deduced from the formula
2
δ = -----------
ωμσ
by setting the coil length, equal to the skin depth δ. In other study types the
conductivity is set to 0.6 S/m.
If None is chosen, no conductivity is used in the coil domain. Choose User defined to
specify the Electrical conductivity in the coil domain σΩ (SI unit: S/m). The default
value is 1 S/m. The purpose of this electrical conductivity is only to stabilize the
solution. According to the Homogenized multi-turn model, the domain should not
be conductive and all the currents should flow in the direction of the wires only.
• Coil Excitation
• Power Excitation
• Using Coils in 3D Models
COIL
A harmonic perturbation can be specified only if the Coil excitation of the parent feature
is Coil current or Coil voltage. In the first case, enter a Coil current Icoil (SI unit: A). The
default is 1 A. In the second case, enter a Coil voltage Vcoil (SI unit: V). The default is
1 V.
Coil Geometry
The Coil Geometry subnode of the Coil is available if the Conductor model of the parent
feature is Homogenized multi-turn and the Coil type is Linear or Circular. This feature was
previously called Reference Edge. This subnode is used to specify the direction of the
tiny wires constituting the multi-turn coil.
• In Linear coils, the wires are all parallel and straight lines. Select an edge or a group
of co-linear edges. The direction of the wires and the coil length is taken to be the
direction and the length of the edges. To avoid unphysical currents, a Linear coil
should be terminated on external boundaries.
• In Circular coils, the wires are wound in circles around the same axis. Select a group
of edges forming a circle or a part of a circle around the coil’s axis. From the selected
edges, the coil axis is computed, and the direction of the wires is taken to be the
azimuthal direction around the axis. The coil length used is simply the length of the
edges; the best approximation is obtained when the radius of the edges is close to
the average radius of the coil.
Select the Use robust geometry analysis method check box to use an alternative
algorithm for the computation of the coil axis and direction. This method works
even in models without geometry (for example with an imported mesh), but
requires that the selected edges form a complete circle.
COIL GEOMETRY
• Enter the Coil current flow ecoil, a vector field representing the local direction of the
wires, which can be an analytic expression, or the solution of another physics (for
example, Curvilinear Coordinates). The vector field will be normalized.
• Enter the Domain length Ldomain (for domain features) or the Boundary length
Lboundary (for boundary features). The SI unit is meters. This value should
correspond to the average length of the wires in the coil, before any symmetry
correction.
• If the model represents only a part of a geometry obtained from symmetry cuts (for
example, a quarter of a circle), use the settings under Symmetry specification to
specify appropriate correction factors. Enter the Coil length multiplication factor FL
and Coil area multiplication factor FA (dimensionless integer numbers). The actual
length of the coil, used to compute the coil voltage and resistance, is then computed
as the product FL·Ldomain(or F·Lboundary for boundary features). The cross-section
area of the coil is computed as FA·Adomain.
This node should be used in combination with the Output (for User Defined Coil
Geometry) node.
This feature is available if the Accurate coil voltage calculation functionality is enabled in
the Stabilization section of the parent feature. This functionality performs a current
filtering to improve the accuracy of the computed electric field and coil voltage in
frequency domain. For this filtering to work properly, special conditions must be
defined on the inlet/outlet boundaries, where the wire direction has a nonzero normal
component. For Linear and Circular coils, the geometry analysis algorithm detects
automatically these boundaries, but the Override Coil Inlet/Outlet feature can be used in
case the detection fails. For Numeric coils, or when the Conductor model is Single
conductor, the subnodes of the Geometry Analysis feature are used to identify the
inlets/outlets. For User defined coils, the User Defined Coil Geometry can be used to
select the inlets/outlets instead.
Geometry Analysis
In 3D components, when Conductor model of the parent feature is Homogenized
multi-turn and Numeric is selected as Coil type, or if Conductor model of the parent
feature is Single conductor in the Coil node, the Geometry Analysis subnode is added by
default. To this, an Input (for Geometry Analysis) subnode is also added by default.
This subnode is needed to set up the automatic analysis of the coil geometry that can
determine the direction of the wires (and of the current flow). The inlet and outlet
boundaries for the wires are specified using the Input (for Geometry Analysis) and
Output (for Geometry Analysis) subnodes.
To complete the setup for the computation of the coil direction, a Coil Geometry
Analysis study step must be added to the study before the main study step.
When applied on one or more boundaries, a red arrow in the Graphics window will
show the direction of the wire. Select the Reverse direction check box to reverse the
direction of the wires.
For an open coil, the Input subnode should never be added to boundaries
that are interior to the coil’s domain selection. Note also that an open coil
should have its start (Input) and end (Output) on boundaries that are
exterior to the domain selection of the parent physics (e.g. Magnetic Fields)
and of a type that supports closure of the coil by means of an induced
surface current density. The latter condition means Magnetic Insulation
and/or, in the frequency domain, Impedance Boundary Condition.
For an open coil, the Input subnode should never be added to boundaries
that are interior to the coil’s domain selection. Note also that an open coil
should have its start (Input) and end (Output) on boundaries that are
exterior to the domain selection of the parent physics (e.g. Magnetic Fields)
and of a type that supports closure of the coil by means of an induced
surface current density. The latter condition means Magnetic Insulation
and/or, in the frequency domain, Impedance Boundary Condition.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Figure 5-3: A coil geometry with turns that are geometrically unified, yet electrically
separated. The electric insulation boundary condition can be used to override the electrical
continuity that applies to interior boundaries by default.
• When the Conductor model of the Coil feature is set to Single conductor,
the Electric Insulation affects the externally applied currents only, not
the currents resulting from electromagnetic induction in the model
itself. For this reason, the single conductor with electric insulation
should be used in the (quasi) static regime only.
• When the Conductor model of the Coil feature is set to Homogenized
multi-turn, the Electric Insulation works for higher frequencies as well,
as both external and induced currents are forced along the same path.
• For low frequencies, the Homogenized multi-turn option can be
considered as an approximation of the Single conductor one. For higher
frequencies, you should consider having a gap between the turns, or
using the Contact Impedance boundary feature available in The
Magnetic and Electric Fields Interface.
Connected Boundaries
The Connected Boundaries subnode is available from the context menu (right-click the
Geometry Analysis parent node) or from the Physics toolbar, Attributes menu. It can
The vector field computed in the Coil Geometry Analysis step (which corresponds to
the current flow or the wire direction, according to the Conductor model selected) will
have zero total flux into the selected boundaries. For example, if two non-connected
boundaries are selected, the feature will enforce the condition that the current entering
one boundary must be equal to the current leaving the other boundary.
The feature models an electrical connection between the selected boundaries, and can
therefore be used for coils constituted of multiple series-connected domains.
Single-Turn Coil
The Single-Turn Coil node is available both as a domain and as a boundary feature. It
models a conductive (for example, metallic) domain or boundary subject to a lumped
excitation, such as a total voltage or current. The excitation specified is translated into
a conduction current flowing in the domain (as a current density) or the boundary (as
a surface current density). The Single-Turn Coil domain node is a contributing feature,
it is applied on top of an Ampère’s Law feature that provides the material model
(electrical conductivity and relative permittivity). The boundary node is instead
exclusive, since the material model is provided in its Thin Layer section.
The Single-Turn Coil feature is obsolete in the Magnetic Field and Rotating
Machinery, Magnetic features and will be removed in future versions. The
recommended alternative is to use a Coil feature with the Conductor model
set to Single conductor.
See Coil Features in the modeling section to learn more about using this
node.
This node is not available in 3D for The Magnetic and Electric Fields Interface and
The Induction Heating Interface. In 2D and 2D axisymmetric components, the
direction of the applied electric field is assumed to be out-of-plane. The settings specify
how to compute the electric field and the relative current density.
For 3D components add boundary conditions for the coil potential. The
Gap Feed (at the domain level only), Boundary Feed, Ground, and
Floating Potential subnodes are available from the context menu
(right-click the parent node) or from the Physics toolbar, Attributes menu.
Use these to define the geometry of the coil.
Select a Coil conductivity σcoil (SI unit: S/m) — From material (the default) or User
defined. For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic based on
the characteristics of the coil conductivity, and then enter values or expressions in the
field or matrix. The default is 1 S/m.
Select a Coil relative permittivity εr,coil (dimensionless) — From material (the default)
or User defined. For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic
based on the characteristics of the coil relative permittivity, and then enter values or
expressions in the field or matrix. The default is 1.
Enter the Thickness ds (SI unit: m) of the thin conductive layer. The default is 1 mm.
SINGLE-TURN COIL
3D Components
For 3D components, the parameters Coil conductivity and Coil relative permittivity
define the material model to be used with the current continuity equation. Normally,
these parameters should be the same used in the Ampère’s Law feature active in the
domain.
Select a Coil conductivity σcoil (SI unit: S/m) — From material or User defined. For User
defined select Isotropic, Diagonal, Symmetric, or Anisotropic based on the characteristics
of the coil conductivity, and then enter values or expressions in the field or matrix.
Select a Coil relative permittivity εr,coil (dimensionless) — From material or User defined.
For User defined select Isotropic, Diagonal, Symmetric, or Anisotropic based on the
characteristics of the coil relative permittivity, and then enter values or expressions in
the field or matrix.
When defining this node on a boundary, also enter a Thickness ds (SI unit: m). The
default is 1 mm.
The settings specifying coil name and excitation are the similar to the ones for the Coil
node. See the section Coil for a description of these settings.
Gap Feed
The Gap Feed subnode is available from the context menu (right-click the Single-Turn
Coil parent node) or from the Physics toolbar, Attributes menu for 3D components (at
the domain level). Gap Feed is used to excite the coil from an internal boundary. It
applies a discontinuity in the coil potential across a boundary to enforce a potential
difference or a prescribed current. This feature must be applied to an internal
boundary in the coil domain, and is typically used to excite a closed loop.
SINGLE-TURN COIL
See Single-Turn Coil for all settings.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Boundary Feed
The Boundary Feed subnode is available from the context menu (right-click the
Single-Turn Coil parent node) or from the Physics toolbar, Attributes menu for 3D
components. The Boundary Feed node excites the coil acting on the coil potential
variable. The feature must be applied to an external boundary.
SINGLE-TURN COIL
See Single-Turn Coil for all settings.
Ground
The Ground subnode is available from the context menu (right-click the Single-Turn
Coil parent node) or from the Physics toolbar, Attributes menu for 3D components.
The Ground subnode enforces the condition Vcoil = 0 on a boundary.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Floating Potential
The Floating Potential subnode is available from the context menu (right-click the
Single-Turn Coil parent node) or from the Physics toolbar, Attributes menu. It can be
used to model a boundary at a constant, unknown, coil potential. If applied on
multiple non connected boundaries, it models an electrical connection between them.
The feature can be used for coils constituted of multiple series-connected domains.
Unlike the corresponding feature in other AC/DC physics, this feature only allows
zero total current flowing in the boundary. This is to ensure that all currents in the
coils are excited by the Boundary Feed or Gap Feeds and are properly accounted for in
the postprocessing variables.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Coil Groups
Use this subnode to mark domains or boundaries in which the coil current flows in a
reversed direction, that is antiparallel to the out-of-plane direction (z or azimuthal
direction). This subnode can be used in any combination with the Domain
Group/Boundary Group subnode.
Coil Groups
Lumped Port
Use the Lumped Port condition to apply a uniform electric field between two metallic
boundaries. The excitation at the port can be expressed as a voltage or as a current, or
via the connection to a circuit interface. The use of this feature is justified in
frequency-domain studies when the distance between the metallic boundaries is much
smaller than the wavelength of the electromagnetic radiation. Subject to the same
See S-Parameters and Ports and Lumped Ports with Voltage Input for
more information.
The geometry of the port is specified by the Type of port. A Uniform lumped port
applies a constant electric field between the metallic electrodes. A Coaxial lumped port
applies a radial electric field between two concentric circular metallic boundaries. For
these two cases, the dimension of the port is computed automatically by analyzing the
geometry. Generally, select User defined to manually specify the direction between the
lumped port terminals ah and the dimensions of the port: hport (Height of lumped port)
and wport (Width of lumped port), both with SI unit: m.
Type of Port
The geometry of the port is specified by the Type of port. A Uniform lumped port
applies a constant electric field between the metallic electrodes. A Coaxial lumped port
applies a radial electric field between two concentric circular metallic boundaries. For
these two cases, the dimension of the port is computed automatically by analyzing the
geometry.
Select User defined for non uniform ports, for example, a curved port and enter values
or expressions in the fields for these additional settings:
Select On or Off from the Wave excitation at this port list to set whether it
is an inport or a listener port (connected to a passive transmission line).
For On enter a Voltage V0 (SI unit: V), and Port phase θ (SI unit: rad).
The section is deactivated if a Port sweep is active, since in this case the
sweep controls the ports.
For the purpose of computing S-parameters, only one port at the time
should be excited, for example by performing a Port sweep. In other cases,
for example, when studying microwave heating, more than one inport
might be desired, but in this case the S-parameter variables cannot be
correctly computed. If several ports are excited at the same time, the
S-parameter variables are not available for postprocessing.
SETTINGS
• For Cable enter the Characteristic impedance Zref (SI unit: Ω).
• For Current enter a current I0 (SI unit: A).
Lumped Element
Use a Lumped Element node to mimic the insertion of a resistor, capacitor, inductor or
general impedance between two metallic boundaries. It is similar to the Lumped Port
but does not allow for an excitation and does not generate S-parameter output.
• For User defined enter a Lumped element impedance Zelement (SI unit: Ω; default
50 Ω.)
• For Inductor enter a Lumped element inductance Lelement (SI unit: H; default 1 nH).
• For Capacitor enter a Lumped element capacitance Celement (SI unit: F; default 1 pF).
Edge Current
Use the Edge Current node to specify a line current along one or more edges. A positive
edge current flows in the direction marked by the red arrow.
EDGE CURRENT
Specify the Edge current I0 (SI unit: A). The default is 0 A.
The External Magnetic Vector Potential boundary condition forces the reduced
magnetic vector potential to be zero on the boundary, or, equivalently, forces the total
field to be equal to the background field. Apply this boundary condition on external
boundaries that are at a distance far enough from the system so that its effect on the
background field is negligible.
The boundary condition approximates this penetration to avoid the need to include
another domain in the model. The material properties that appear in the equation are
those for the conductive material excluded from the model.
The skin depth (that is, the distance where the electromagnetic field has decreased by
a factor e−1) is for a good conductor
2
δ = -----------
ωμσ
The impedance boundary condition is a valid approximation if the skin depth is small
compared to the size of the conductor. The source electric field Es can be used to
specify a source surface current on the boundary.
Js
J
I'
I
Js
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or on the Physics toolbar, click the
Attributes menu and select Harmonic Perturbation. For more information see Harmonic
Perturbation — Exclusive and Contributing Nodes in the COMSOL Multiphysics
Reference Manual.
Based on space dimension, enter coordinate values or expressions for the Source electric
field Es (SI unit: V/m).
( 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 defaults use the values From material, taking the properties from the material
specified for the boundary. For User defined, enter different values or expressions.
d
n × ( H 1 – H 2 ) = ∇ t × ------------ ∇ t × A
μ0 μr
Enter a Surface thickness ds (SI unit: m) for the gap. The default is 1 m.
DIPOLE PARAMETERS
If Magnitude and direction is selected under Dipole Specification, enter coordinates for
the Magnetic dipole moment direction nm and the Magnetic dipole moment, magnitude m
(SI unit: m2⋅A).
If Dipole moment is selected under Dipole Specification, enter coordinates for the
Magnetic dipole moment m (SI unit: m2⋅A).
Magnetic Shielding
The Magnetic Shielding node adds a boundary condition for magnetic shielding. It
describes a thin layer of a highly permeable medium that shields the magnetic field.
Because of the presence of this thin layer, the tangential components of the magnetic
field and vector potential have discontinuities given by:
n × ( H1 – H2 ) = n × Ht ( Bt )
A1 – A2
B t = n × --------------------
ds
In this equation, ds is the surface layer thickness and the relation between the
tangential magnetic field Ht and flux density Bt inside the layer can be linear and given
by a relative permeability or nonlinear and given by the H-B curve of the layer material.
Relative Permeability
Select Relative permeability μr (dimensionless) to use the constitutive relation
B = μ0μrH. The default uses values From material. For User defined enter another value
or expression.
H-B Curve
Select H-B curve |H| (SI unit: A/m) to use a curve that relates magnetic flux density B
and the magnetic field H as |H| = f(|B|). The Magnetic field norm setting can take the
values From material or User defined.
The physics interface solves Maxwell’s equations formulated using the magnetic field
as the dependent variable.
The main node is Faraday’s Law, which adds the equation for the magnetic field and
provides an interface for defining the constitutive relations and its associated properties
such as the relative permeability.
When this physics interface is added, these default nodes are also added to the Model
Builder — Magnetic Field Formulation, Faraday’s Law, Magnetic Insulation (the default
boundary condition), Initial Values, and Magnetic Gauss’ Law. Then, from the Physics
toolbar, add other nodes that implement, for example, boundary conditions and
external currents. You can also right-click Magnetic Field Formulation 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 mfh.
The electric currents and the electric field vectors are orthogonal to the
magnetic field. This setting also controls the direction in which applied
and induced currents can flow in the model. The default option is to solve
for the in-plane components only, allowing an out-of-plane current flow.
Select Components — Out-of-plane magnetic field, In-plane magnetic field (the default),
or Three-component magnetic field for the magnetic field. This choice affect also the
directions in which the electric current is allowed to flow: the magnetic field and
electric currents are orthogonal.
DEPENDENT VARIABLES
The dependent variable is the Magnetic field H. You can change both its field name and
the individual component variable names. If the new field name coincides with the
name of another magnetic vector potential field in the model, the interfaces share
degrees of freedom and component names. The new field name must not coincide with
the name of a field of another type or with a component name belonging to some other
field. Component names must be unique within a model, except for fields of the same
type sharing a common field name.
DIVERGENCE CONSTRAINT
To display this section, click the Show button ( ) and select Advanced Physics Options.
The Activate divergence constraint check box is selected by default, which enables the
default Magnetic Gauss’ Law node. Click to clear the check box to disable the
Magnetic Gauss’s Law node and remove it from the physics interface.
ERROR CHECK
To display this section, click the Show button ( ) and select Advanced Physics Options.
When the Check applicability of features in study check box is selected, any features that
are incompatible with the study will generate an error message when trying to solve or
show the default solver. No solver will be generated. Deselect it and you will be able
to run the model, possibly with runtime errors instead. It is available to allow the
advanced user to tweak any feature and use it outside of its intended study scope.
• Domain, Boundary, Point, and Pair Nodes for the Magnetic Field
Formulation Interface
• Theory for the Magnetic Field Formulation Interface
Domain, Boundary, Point, and Pair Nodes for the Magnetic Field
Formulation Interface
The Magnetic Field Formulation Interface has these domain, boundary, point, and
pair nodes 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).
n2 × ( H1 – H2 ) = 0
need to be fulfilled. Because H is being solved for, the tangential component of the
magnetic field is always continuous, and thus the condition is automatically fulfilled.
Table 5-2 lists the interior and exterior boundary conditions available with this physics
interface.
TABLE 5-2: INTERIOR AND EXTERIOR BOUNDARY CONDITIONS FOR THE MAGNETIC FIELD FORMULATION
INTERFACE
Electric Field x x
Impedance Boundary Condition x
Magnetic Field x x
Magnetic Insulation x x
These nodes are described for the Magnetic Fields or Electrostatics interface (listed in
alphabetical order):
In the COMSOL Multiphysics Reference Manual see Table 2-3 for links
to common sections and Table 2-4 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.
Faraday’s Law
The Faraday’s Law node adds Faraday’s law for the electric field and provides an
interface for defining the constitutive relation and its associated properties such as the
relative permeability as well as electric properties.
In the Magnetic Field section, the B-H curve and Effective B-H curve
constitutive relations replace respectively the H-B curve and the Effective
H-B curve constitutive relations available in Ampère’s Law.
Initial Values
The Initial Values node adds an initial value for the magnetic field H that can serve as
an initial value for a transient simulation or as an initial guess for a nonlinear solver.
INITIAL VALUES
Enter values or expressions for the Magnetic field H (SI unit: A/m) components. The
defaults are 0 A/m.
∇⋅B = 0
on the magnetic flux density by introducing an auxiliary dependent variable ψ with the
corresponding equation.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
This section allows a more fine control on the boundary conditions for ψ applied by
the Magnetic Gauss’ Law feature. The domain equation for ψ only imposes a condition
on the gradient, so it is important to constrain the absolute value of ψ to ensure a
nonsingular model. Select the Constrain variable in at least one point to ensure that
there is always a constraint set on the value of the divergence condition variable ψ. The
feature also sets up automatically the appropriate boundary conditions for ψ. Use the
Method combo box to select the approach used to enforce the boundary condition:
Constrain value (the default) or Constrain tangential gradient.
Magnetic Insulation
The Magnetic Insulation node is the default boundary condition for the Magnetic Field
Formulation interface. In this formulation, this condition corresponds to setting the
tangential components of the electric field to zero at the boundary n × E = 0.
Electric Field
The Electric Field boundary condition
n × E = n × E0
specifies the tangential component of the electric field. This feature can be used to
specify a source electric field on the boundary of the conductor.
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or on the Physics toolbar, click the
Attributes menu and select Harmonic Perturbation. For more information see Harmonic
Perturbation — Exclusive and Contributing Nodes in the COMSOL Multiphysics
Reference Manual.
– n × E = J ms
n × ( E 1 – E 2 ) = J ms
The physics interface solves Gauss’ Law for the magnetic field using the scalar magnetic
potential as the dependent variable.
The main node is the Magnetic Flux Conservation feature, which adds the equation
for the magnetic potential and provides an interface for defining the material
properties and the constitutive relation for the magnetic flux density. It is used when
there are no currents and all the magnetic fields are originated by permanent magnets
or external systems not included in the model.
When this physics interface is added, these default nodes are also added to the Model
Builder — Magnetic Fields, No Currents, Magnetic Flux Conservation, Magnetic Insulation
(the default boundary condition), and Initial Values. Then, from the Physics toolbar,
add other nodes that implement, for example, additional boundary conditions and
point conditions. You can also right-click Magnetic Fields, No Currents to select physics
features from the context menu.
Except where described below, the Settings windows are described for the
Magnetic Fields and Electrostatics interfaces.
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 mfnc.
When a background field is active, the physics interface solves for the reduced
(perturbation) field only. This setting, together with the External Magnetic Flux
Density feature, is useful to introduce an external field generated by systems not
included in the model.
THICKNESS
ERROR CHECK
To display this section, click the Show button ( ) and select Advanced Physics Options.
When the Check applicability of features in study check box is selected, any features that
are incompatible with the study will generate an error message when trying to solve or
show the default solver. No solver will be generated. Deselect it and you will be able
to run the model, possibly with runtime errors instead. It is available to allow the
advanced user to tweak any feature and use it outside of its intended study scope.
DEPENDENT VARIABLES
The dependent variable (field variable) is for the Magnetic scalar potential Vm. The
name can be changed but the names of fields and dependent variables must be unique
within a model.
• Domain, Boundary, Point, and Pair Nodes for the Magnetic Fields, No
Currents Interface
• Theory for the Magnetic Fields, No Currents Interface
Domain, Boundary, Point, and Pair Nodes for the Magnetic Fields,
No Currents Interface
The Magnetic Fields, No Currents Interface has these domain, boundary, point, and
pair nodes 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).
n2 ⋅ ( B1 – B2 ) = 0
This condition provides continuity of the normal component of the magnetic flux
density and is automatically satisfied by the natural boundary condition for interior
boundaries, which is
AVAILABLE NODES
These nodes are available for this physics interface and listed in alphabetical order. Also
see Table 5-3 for a list of interior and exterior boundary conditions.
Table 5-3 lists the interior and exterior boundary conditions available with this physics
interface.
TABLE 5-3: INTERIOR AND EXTERIOR BOUNDARY CONDITIONS FOR THE MAGNETIC FIELDS, NO CURRENTS
INTERFACE
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation and other variables
evaluated in a material reference configuration (material frame). Select Nonsolid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
MAGNETIC FIELD
See the settings for Magnetic Field under Ampère’s Law, except for these
differences:
Initial Values
The Initial Values node adds an initial value for the magnetic scalar potential that can
serve as an initial guess for a nonlinear solver.
INITIAL VALUES
Enter a value or expression for the initial value of the Magnetic scalar potential Vm. The
default value is 0 A.
V m = V m0
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
n ⋅ B = n ⋅ B0 (5-1)
Using this boundary condition specify the normal component of the magnetic flux
density at the boundary.
Alternatively, specify an inward (or outward) flux density using the following equation:
–n ⋅ B = Bn (5-2)
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
The External Magnetic Flux Density boundary condition forces the reduced magnetic
flux density to be zero on the boundary, or, equivalently, forces the total field to be
equal to the background field. Apply this boundary condition on external boundaries
that are at a distance far enough from the system so that its effect on the background
field is negligible.
Magnetic Shielding
The Magnetic Shielding node adds a boundary condition for magnetic shielding. It
describes a thin layer of a highly permeable medium that shields the magnetic field.
n ⋅ ( B1 – B2 ) = ∇t ⋅ ( ds Bt ( Ht ) )
Ht = –∇t Vm
MAGNETIC SHIELDING
Enter a value or expression for the Surface thickness ds (SI unit: m).
Relative Permeability
Select Relative permeability μr (dimensionless) to use the constitutive relation
B = μ0μrH. The default uses values From material. For User defined enter another value
or expression.
B-H Curve
Select B-H curve |B| (SI unit: T) to use a curve that relates magnetic field H and the
magnetic flux density B as |B| = f(|H|). The Magnetic flux density norm setting can take
the values From material or User defined.
u d
V m – V m = V m, d
where Vmu is the magnetic scalar potential on the upside of the boundary selection and
d u
V m – V m
n ⋅ ( B 1 – B 2 ) = μ 0 μ r -----------------------
ds
where Vmu is the magnetic scalar potential on the upside of the boundary selection and
Archie’s Law
See the Archie’s Law section in the Electric Field Interfaces chapter.
– ∇ ⋅ ( μ 0 μ r ∇V m ) = 0 .
For a discussion about the boundary element method, see Theory for the
Boundary Elements PDE in the COMSOL Multiphysics Reference
Manual.
When this physics interface is added, these default nodes are also added to the Model
Builder — Magnetic Flux Conservation, Magnetic Insulation (on exterior boundaries) and
Initial Values. Then, from the Physics toolbar, add other nodes that implement, for
example, magnetic scalar potential and magnetic flux density conditions. You can also
right-click Magnetic Fields, No Currents, Boundary Elements to select physics features
from the context menu.
DOMAIN SELECTION
From the Selection list, select any of the options — Manual, All domains, All voids, or All
domains and voids (the default). The geometric entity list displays the selected domain
entity numbers. Edit the list of selected domain entity numbers using the selection
toolbar buttons to the right of the list or by selecting the geometric entities in the
Graphics window. Entity numbers for voids can be entered by clicking the Paste ( )
button in the selection toolbar and supplying the entity numbers in the in the dialog
box. The entity number for the infinite void is 0, and finite voids have negative entity
numbers.
Selections can also be entered using the Selection List window, available from the
Windows menu on the Home toolbar.
PHYSICS SYMBOLS
Select the Enable physics symbols check box to display symmetry planes (in 3D) and
lines (in 2D) in the Graphics window, as specified in the Symmetry settings.
SYMMETRY
For more information about the Far Field Approximation settings, see
Far-Field Approximation Settings in the COMSOL Multiphysics
Reference Manual.
QUADRATURE
To display this section, click the Show button ( ) and select Advanced Physics Options.
INFINITY CONDITION
DISCRETIZATION
From the Magnetic scalar potential/Magnetic flux density list, choose from predefined
options for the boundary element discretization order for the magnetic scalar potential
variable and the magnetic flux density variable, respectively. The predefined options
represent the suitable combinations of element orders such as Quadratic/Linear (the
default).
The settings under Value types when using splitting of complex variables are important
for sensitivity and optimization computations. See the description of the built-in
operators fsens and fsensimag.
AVAILABLE NODES
These nodes and subnodes are described in this section:
See the settings for Magnetic Field under Ampère’s Law. Notice that the
Magnetic Fields, No Currents, Boundary Elements interface supports the
following constitutive relations only: Relative permeability, Remanent flux
density, and Magnetization.
Magnetic Insulation
The Magnetic Insulation node for the Magnetic Fields, No Currents, Boundary
Elements interface provides the default boundary condition. The magnetic insulation
boundary condition sets the normal component of the magnetic flux density to zero:
n⋅B = 0
Initial Values
Initial Values adds an initial value for the magnetic scalar potential that can serve as an
initial condition for a transient simulation or as an initial guess for a nonlinear solver.
Add more Initial Values nodes from the Physics toolbar, Global menu.
INITIAL VALUES
Enter a value or expression for the initial value of the Magnetic scalar potential Vm. The
default value is 0 A.
V m = V m0
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
See the settings for Magnetic Flux Density in the Magnetic Fields, No
Currents documentation. Notice that the Magnetic Fields, No Currents,
Boundary Elements interface supports the Inward flux density option only.
The Magnetic Fields, No Currents Interface can be set to solve for the
Reduced field using the Background Magnetic Field setting. The Magnetic
Fields, No Currents, Boundary Elements Interface does not have such a
setting but it can be made to represent an open boundary for the reduced
field by adding, in The Magnetic Fields, No Currents Interface, an
External Magnetic Flux Density feature on the boundary selection for the
Magnetic Scalar-Scalar Potential Coupling.
BOUNDARY SELECTION
Select Manual or All boundaries from the Selection list. Make additional edits to the list
of boundary entity numbers using the Selection toolbar buttons. When All boundaries
is selected from the Selection list, the boundaries exterior to the Magnetic Fields, No
Currents interface that intersect the exterior boundaries to the Magnetic Fields, No
Currents, Boundary Elements interface are available in the boundary entity number
list.
COUPLED INTERFACES
Select Magnetic Fields, No Currents as Primary interface and Magnetic Fields, No Currents,
Boundary Elements as Secondary interface.
The Magnetic Fields Interface can be set to solve for the Reduced field
using the Background Magnetic Field setting. The Magnetic Fields, No
Currents, Boundary Elements Interface does not have such a setting but
it can be made to represent an open boundary for the reduced field by
adding, in The Magnetic Fields Interface, a Magnetic Field feature with
the background magnetic field on the boundary selection for the Magnetic
Vector-Scalar Potential Coupling.
BOUNDARY SELECTION
Select Manual or All boundaries from the Selection list. Make additional edits to the list
of boundary entity numbers using the Selection toolbar buttons. When All boundaries
is selected from the Selection list, the boundaries exterior to the Magnetic Fields
interface that intersect the exterior boundaries to the Magnetic Fields, No Currents,
Boundary Elements interface are available in the boundary entity number list.
When this physics interface is added, these default nodes are also added to the Model
Builder — Electric Field Transformation, Ampère’s Law, Mixed Formulation Boundary,
Magnetic Insulation (the default boundary condition), and Initial Values. Then, from the
Physics toolbar, add other nodes that implement, for example, boundary conditions
and point conditions. You can also right-click Rotating Machinery, Magnetic 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.
BACKGROUND FIELD
Select an option from the Solve for list — Full field (the default) or Reduced field. For
Reduced field specify a Background magnetic vector potential Ab (SI unit: Wb/m). The
entered expressions must be differentiable.
The total field used in the physics and equations are given by the sum of the reduced
and background fields.
COMPONENTS
For 2D components, select the Components — Out-of-plane vector potential (the
default), In-plane vector potential, or Three-component vector potential.
THICKNESS
For 2D components, enter a value or expression for the global Out-of-plane thickness d
(SI unit: m). The default value of 1 m is typically not representative for a thin domain.
Instead it describes a unit thickness that makes the 2D equation identical to the
equation used for 3D components.
PHYSICS-CONTROLLED MESH
Select the Enable check box to allow the physics interface to control the meshing
process. Information from the physics, such as the presence of infinite elements
domain or periodic condition, will be used to set up automatically an appropriate mesh
sequence.
ERROR CHECK
To display this section, click the Show button ( ) and select Advanced Physics Options.
When the Check applicability of features in study check box is selected, any features that
are incompatible with the study will generate an error message when trying to solve or
show the default solver. No solver will be generated. Deselect it and you will be able
to run the model, possibly with runtime errors instead. It is available to allow the
advanced user to tweak any feature and use it outside of its intended study scope.
DEPENDENT VARIABLES
The dependent variables (field variables) are the Magnetic vector potential A and the
Magnetic scalar potential Vm. The name can be changed but the names of fields and
dependent variables must be unique within a model.
DISCRETIZATION
Domain, Boundary, Edge, Point, and Pair Nodes for the Rotating
Machinery, Magnetic Interface
Domain, Boundary, Edge, Point, and Pair Nodes for the Rotating
Machinery, Magnetic Interface
The Rotating Machinery, Magnetic Interface has these boundary, edge, pair, and point
nodes 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).
In the COMSOL Multiphysics Reference Manual see Table 2-3 for links
to common sections and Table 2-4 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.
These boundary, edge, pair, and point nodes are described for the main physics
interfaces that solve for the magnetic vector potential (Magnetic Fields) and the
Initial Values
The Initial Values node adds an initial value for the magnetic vector potential.
INITIAL VALUES
Enter a value or expression for the initial value of the Magnetic vector potential A
(SI unit: Wb/m).
PRESCRIBED ROTATION
Enter the Rotation angle α rot (SI unit: radians).
• For 2D components, enter the X and Y coordinates for the Rotation axis base point
rax (SI unit: m) identifying the fixed point of the rotation. The axis of rotation is
the z (out-of-plane) axis.
• For 3D components, specify the rotation axis by entering the coordinates of the
Rotation axis base point rax (SI unit: m) and the components of the Rotation axis
urot (dimensionless).
• For 2D components, enter the X and Y coordinates for the Rotation axis base point
rax (SI unit: m) identifying the fixed point of the rotation. The axis of rotation is
the z (out-of-plane) axis.
• For 3D components, specify the rotation axis by entering the coordinates of the
Rotation axis base point rax (SI unit: m) and the components of the Rotation axis
urot (dimensionless).
Magnetic Insulation
The Magnetic Insulation node is the default external boundary condition for the
Rotating Machinery, Magnetic interface. It adds a boundary condition that, on
boundaries with magnetic vector potential (Ampère’s Law domains) sets the tangential
For the magnetic vector potential, Magnetic insulation is a special case of the magnetic
potential boundary condition that sets the tangential component of the magnetic
potential to zero and can be used for the modeling of a lossless metallic surface, for
example, a ground plane or as a symmetry type boundary condition. The Magnetic
Insulation node can also be applied on interior boundaries in vector potential domains.
For more details on this, see also the corresponding documentation section in the
Magnetic Fields interface on Magnetic Insulation.
The physics interface solves Maxwell’s equations formulated using the magnetic vector
potential and the scalar electric potential as the dependent variables.
Use the Magnetic and Electric Fields interface when the current continuity equation
is needed everywhere in the simulation domain.
The main node is the Ampère’s Law and Current Conservation feature, which adds the
equation for the electric potential and magnetic vector potential and provides an
interface for defining the constitutive relations and their associated properties such as
the relative permeability, relative permittivity, and electrical conductivity.
The Ampère’s Law and Current Conservation feature can be overridden with an
Ampère’s Law feature or a Coil feature, that are identical to the ones in the Magnetic
Fields interface, removing the electric potential and the current conservation equation
from the selected domains.
When this physics interface is added, these default nodes are also added to the Model
Builder— Ampère’s Law and Current Conservation, Magnetic Insulation (the default
boundary condition for the magnetic vector potential), and Initial Values. Then, from
the Physics toolbar, add other nodes that implement, for example, boundary
Except where described below, most Settings windows are the same as for
the Magnetic Fields, Electrostatics, and Electric Currents interfaces. Also
see About the Magnetic and Electric Field Interface Boundary
Conditions for more information.
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 mef.
BACKGROUND FIELD
The only option available from the Solve for list is Full field.
COMPONENTS
THICKNESS
Select the Activate input sweep check box to switch on the sweep and invoke a
parametric sweep over the Lumped ports or the Terminal nodes.
Select an option from the Sweep on list — Terminals or Ports. This setting controls
which features are activated during the sweep.
Enter a Sweep parameter name to indicate the name of the model parameter that
controls the terminal or port activated in each step of the sweep. The default is
PortName. The name given must match the model parameter, defined under Global
Definitions, that is the object of a Parametric Sweep node in the current Study.
When Terminals is selected under Sweep on, also select a Parameter to export — Z (the
default), Y, or S. When the sweep is carried on Ports, the exported parameter is always S.
ERROR CHECK
To display this section, click the Show button ( ) and select Advanced Physics Options.
When the Check applicability of features in study check box is selected, any features that
are incompatible with the study will generate an error message when trying to solve or
show the default solver. No solver will be generated. Deselect it and you will be able
to run the model, possibly with runtime errors instead. It is available to allow the
advanced user to tweak any feature and use it outside of its intended study scope.
DEPENDENT VARIABLES
The dependent variables (field variables) are for the Electric potential V and Magnetic
vector potential A. The name can be changed but the names of fields and dependent
variables must be unique within a model.
The basic steps for this are to add the magnetic boundary conditions from
the Physics toolbar. Then right-click these magnetic boundary condition
nodes to add compatible electric boundary conditions as subnodes.
n2 × ( A1 – A2 ) = 0
n2 × ( H1 – H2 ) = 0
need to be fulfilled. Because the physics interface solves for A, the tangential
component of the magnetic potential is always continuous, and thus the first condition
is automatically fulfilled. The second condition is equivalent to the natural boundary
condition and is hence also fulfilled unless surface currents are explicitly introduced.
n2 ⋅ ( J1 – J2 ) = 0
A default Electric Insulation node is also added under the default Magnetic
Insulation node. This is to ensure that default boundary conditions are
always specified both for the magnetic vector potential and the electric
potential. Additional subnodes are available from the context menu
(right-click the parent node) or from the Physics toolbar, Attributes menu.
Table 5-4 lists the interior and exterior boundary conditions available with this physics
interface.
TABLE 5-4: INTERIOR AND EXTERIOR ELECTRIC AND MAGNETIC BOUNDARY CONDITIONS FOR THE
MAGNETIC AND ELECTRIC FIELDS INTERFACE
MAGNETIC BOUNDARIES
Contact Impedance x
Electric Insulation x x
Electric Potential x x
Electric Shielding x
Floating Potential x x
Ground x x
Normal Current Density x
Periodic Condition x
Domain, Boundary, Edge, Point, and Pair Nodes for the Magnetic
and Electric Fields Interface
The Magnetic and Electric Fields Interface has these domain, boundary, edge, point,
and pair nodes 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).
Also see About the Magnetic and Electric Field Interface Boundary Conditions and
Table 5-1 for a list of interior and exterior boundary conditions.
• Sector Symmetry
• Infinite Elements, Perfectly Matched Layers, and Absorbing Layers in
the COMSOL Multiphysics Reference Manual
In the COMSOL Multiphysics Reference Manual see Table 2-3 for links
to common sections and Table 2-4 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.
MATERIAL TYPE
The Material type setting decides how materials behave and how material properties are
interpreted when the mesh is deformed. Select Solid for materials whose properties
change as functions of material strain, material orientation and other variables
evaluated in a material reference configuration (material frame). Select Nonsolid for
materials whose properties are defined only as functions of the current local state at
each point in the spatial frame, and for which no unique material reference
configuration can be defined. Select From material to pick up the corresponding setting
from the domain material on each domain.
The Conduction Current and Magnetic Field settings are the same as
Ampère’s Law for the Magnetic Fields interface. The Electric Field settings
are the same as for Current Conservation for the Electric Currents
interface.
Initial Values
The Initial Values node adds an initial value for the magnetic vector potential and
electric potential that can serve as an initial value for a transient simulation or as an
initial guess for a nonlinear solver.
INITIAL VALUES
Enter a value or expression for the initial value of the Magnetic vector potential A
(SI unit: Wb/m) and initial value of the Electric potential V (SI unit: V). The default
values are 0 Wb/m and 0 V, respectively.
Surface Current
The Surface Current node adds a boundary condition for a surface current density Js.
An externally applied surface current density can be specified directly in the Settings
window for the node. Alternatively, the Electric Shielding subnode is available from the
context menu (right-click the parent node) or from the Physics toolbar, Attributes
menu. Use the subnode to model a thin conductive layer (a layer whose electrical
conductivity s is higher than in the surrounding domains). In this case, the electric
current flows preferably along the conductive boundary, and the resulting surface
current density is automatically incorporated in the parent Surface Current node to act
as a source for the magnetic field.
SURFACE CURRENT
Enter values or expressions for the Surface current density Js0 (SI unit: A/m)
coordinates.
Magnetic Shielding
The Magnetic Shielding node adds a boundary condition for magnetic shielding. By
itself, this feature represents a geometrically thin layer with a relative permeability
substantially higher than the surrounding material, that shields the magnetic field.
Subnodes are available from the context menu (right-click the parent node) or from
the Physics toolbar, Attributes menu. Use the subnodes to provide an electrical model
• Magnetic Shielding with no subnode: this represents a highly magnetic layer with no
special electrical properties.
• Magnetic Shielding and Electric Shielding: this combination represents a layer in which
both the relative permeability μr and the electrical conductivity σ are much higher
than in the adjacent domain. This layer shields the magnetic field and the current
preferentially flows along it. The resulting surface current density is incorporated in
the model as a source for the magnetic field.
• Magnetic Shielding and Electric Insulation: this combination represents a
nonconductive magnetic layer. The layer shields the magnetic field and there is no
current flowing through it. The electric potential is discontinuous across the
boundary.
• Magnetic Shielding and Contact Impedance: this combination represents a resistive
magnetic layer. The layer shields the magnetic field and the current flows preferably
orthogonally to the layer. The electric potential is discontinuous across the
boundary.
MAGNETIC SHIELDING
The default Relative permeability μr (dimensionless) uses values From material. For User
defined select Isotropic, Diagonal, Symmetric, or Anisotropic and enter other values or
expressions. For anisotropic material, the relative permeability is a tensor.
Enter a value or expression for the Surface thickness ds (SI unit: m).
Magnetic Continuity
The Magnetic Continuity node does not apply any condition or constraint to the model.
Its purpose is to allow the application of electrical boundary conditions by means of
the subnodes Electric Insulation and Contact Impedance. The subnodes are available
from the context menu (right-click the parent node) or from the Physics toolbar,
Attributes menu. Use this feature to model thin layer with particular electrical
properties but with the same magnetic properties of the surrounding medium.
Apply this feature to a group of domains representing the cross sections of the coil
turns on the modeling plane. Since the electric potential is assumed to be constant in
the cross section of each turn, the V variable must be removed in the selected domains
by applying an Ampère’s Law node. The feature applies an external current density
flowing orthogonally to the plane, and also imposes a voltage constraint on the
boundary of each cross section, computed from the coil and the excitation properties.
The feature also enforces a balance of the current flowing out-of-plane and the current
leaking in the plane between the coil turns.
DOMAINS LATTICE
In order to apply the correct voltage at each coil turn, it is necessary to precisely define
the order in which the domains are connected. The RLC Coil Group feature assumes
that all domains are of the same shape and are disposed on a two-dimensional Bravais
lattice; that is, that there exist two vectors a and b (called primitive vectors) such that
given two lattice points r1 and r2 the following relation holds:
r 1 – r 2 = na + mb
The RLC Coil Group node can automatically determine two lattice vectors (or one, in
the case of a linear lattice) from the geometry and position of the selected domains.
The choices Row-wise and Column-wise in the Ordering settings use this functionality.
The two choices control how the vectors a and b are chosen.
Since the primitive vectors of a Bravais lattice are not unique, the lattice recognition
algorithm could find two vectors different from the desired one. In this case, it is
possible to manually specify User defined lattice vectors.
DOMAIN ORDERING
By default, the RLC Coil Group feature assigns the first turn index (closest to the
reference potential) to the turn with the smallest n and m. When using the automatic
recognition algorithm, this usually corresponds to the coil turn with smallest x and y
(or r and z in axisymmetry). The other coil turns are numbered in order of increasing
m first, and then increasing n. Refer to the diagram in the node’s Settings window for
a visual explanation.
d
Vi = Vi – 1 – V1
where V0 is the value specified as the Ground voltage in the Settings window, that drives
an out-of-plane external current density computed as in the Single-Turn Coil case:
d
Vi
J e = σ ------- e z in 2D
d
d
Vi
J e = σ ---------- e φ in 2D axisym.
2πr
The current flowing in-plane is balanced with the difference of the currents between
two adjacent turns. If Ii is the out-of-plane current flowing in the i-th turn,
Ii = (J ⋅ e)
Ωi
Ii – ( J ⋅ n ) – Ii + 1 = 0 i = 1, …, n – 1
∂Ω i
In – ( J ⋅ n ) – I Coil = 0
∂Ω n
The coil current ICoil is computed from the coil excitation in the same way as for the
Single-Turn Coil feature.
GEOMETRY
This section specifies the information needed to identify the Bravais lattice on which
the coil turns are placed The feature automatically detects the two primitive vectors of
the lattice if Ordering is set as Row-wise (the default) or Column-wise. For User defined
use the diagram as a guide to manually enter the Primitive vectors a and b (SI unit: m)
of the lattice. The diagram shows how the domains are ordered for each choice.
Select a Winding style — Regular (the default) or Alternated. Look at the diagram for an
illustration of these options.
The Magnetic Fields Interfaceis used for 3D, 2D in-plane, and 2D axisymmetric
components. Unless you have a license for the AC/DC Module, only 2D components
involving out-of-plane currents and axisymmetric components involving azimuthal
currents are supported.
In this section:
• Magnetostatics Equation
• Frequency Domain Equation
• Transient Equation
• Maxwell’s Equations
• Magnetic and Electric Potentials
• Gauge Transformations
• Selecting a Particular Gauge
• The Gauge and Equation of Continuity for Dynamic Fields
• Explicit Gauge Fixing/Divergence Constraint
• Ungauged Formulations and Current Conservation
• Time-Harmonic Magnetic Fields
Magnetostatics Equation
To derive the magnetostatic equation, start with Ampère’s law for static cases
∇ × H = J. The current is
J = σv × B + J e
–1
∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) = J e
The term involving the velocity only applies in the 2D and 2D axisymmetric
formulations.
∂D ∂D
∇ × H = J + ------- = σE + σv × B + J e + -------
∂t ∂t
Now assume time-harmonic fields and use the definitions of the fields,
B = ∇×A
E = – j ωA
and combine them with the constitutive relationships B = μ0(H + M) and D = ε0E to
rewrite Ampère’s law as
2 –1
( jωσ – ω ε 0 )A + ∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) = J e
The term involving the velocity only applies in the 2D and 2D axisymmetric
formulations.
Transient Equation
The transient equation this physics interface solves is Ampère’s law, illustrated here
with the constitutive relation B = μ0(H + M).
∂A –1
σ + ∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) = J e
∂t
Maxwell’s Equations
Quasi-static analysis of magnetic and electric fields is valid under the assumption that
∂D/∂t = 0.
This implies that it is possible to rewrite Maxwell’s equations in the following manner:
∇ × H = J = σ ( E + v × B ) + Je
∂B
∇×E = –
∂t
∇⋅B = 0
∇⋅D = ρ
∇⋅J = 0
B = ∇×A
E = – ∇V – ∂ A
∂t
and the constitutive relation B = μ0(H + M), Ampère’s law can be rewritten as
∂A –1
σ + ∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) + σ ∇V = J e (5-3)
∂t
The equation of continuity, which is obtained by taking the divergence of the above
equation, adds the following equation:
∂A
∇ ⋅ – σ ------- + σv × ( ∇ × A ) – σ ∇V + J e = 0 (5-4)
∂t
Gauge Transformations
The electric and magnetic potentials are not uniquely defined from the electric and
magnetic fields through
∂A
E = – ------- – ∇V
∂t
B = ∇×A
˜
A = A + ∇Ψ
˜ ∂Ψ
V = V – -------
∂t
˜ ˜
∂A ∂( A – ∇Ψ ) ˜ ∂Ψ ∂A
E = – ------- – ∇V = – --------------------------- – ∇ V + ------- = – ------- – ∇V
˜
∂t ∂t ∂t ∂t
˜ ˜
B = ∇ × A = ∇ × ( A – ∇Ψ ) = ∇ × A
∂A
∇ ⋅ – σ ------- + J e = 0
∂t
It is clear that unless the electrical conductivity is uniform, the particular gauge used
to eliminate V cannot be the Coulomb gauge as that would violate the equation of
continuity and would thereby also violate Ampère’s law.
–1
∇ × ( μ0 ∇ × A – M ) = Je
The equation for ψ is used to impose the Coulomb gauge: ∇⋅ A = 0. However, to get
a closed set of equations, ψ must be able to affect the first equation and this is obtained
by modifying the first equation to:
–1
∇ × ( μ 0 ∇ × A – M ) = J e + ∇ψ
The additional term on the right-hand side can be seen as a Lagrange multiplier that
not only imposes the Coulomb gauge but also eliminates any divergence in the
externally generated current density, Je and makes it comply with the current
continuity inherent in Ampère’s law.
The gauge fixing feature similarly imposes the Coulomb gauge also for the dynamic
(frequency domain) study type in the Magnetic and Electric Fields interface.
∇⋅J = 0
∇ ⋅ ( σA ) = 0
The main benefit of using this kind of divergence constraint is improved numerical
stability, especially when approaching the static limit when the inherent gauge
deteriorates.
When using the Magnetic and Electric Fields interface the electric potential is used to
state current conservation so unless nonphysical current sources are specified inside the
computational domain current conservation is fulfilled.
When using the Magnetic Fields interface, current conservation is usually imposed
either by the solver (for magnetostatics) or in the transient or time harmonic case by
the induced current density. The explicit gauge or divergence constraint can also help
imposing current conservation as described in Explicit Gauge Fixing/Divergence
Constraint.
∇ × H = J = σ ( E + v × B ) + ( jωD + J e )
In the transient case, the inclusion of this term leads to a second-order equation in
time, but in the harmonic case there are no such complications. Using the definition
of the electric and magnetic potentials, the system of equations becomes:
The constitutive relation D = ε0E + P has been used for the electric field.
To obtain a particular gauge that reduces the system of equation, choose Ψ = −jV/ω in
the gauge transformation. This gives:
˜ j ˜
A = A – ---- ∇V V = 0
ω
˜
When V vanishes from the equations, only the second one is needed,
2 ˜ –1 ˜ ˜
( jωσ – ω ε 0 )A + ∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) = J e + jωP
˜
Working with A is often the best option when it is possible to specify all source
currents as external currents Je or as surface currents on boundaries.
Equation System
To derive the equation system in this physics interface, start with Faraday’s law:
∂B
∇ × E = – -------
∂t
The electric field is computed from the current density, using an appropriate
constitutive relation:
E = Em ( J – Je ) + v × B
∇×H = J
Combining these equations give the first equation in the system applied by the
Magnetic Field Formulation interface:
∂B
∇ × ( E m ( J – J e ) + v × B ) + ------- = 0
∂t
∇ × ( Em ( J – Je ) + v × B ) = 0
∇ × ( E m ( J – J e ) + v × B ) + jωB = 0
∇⋅B = 0
• Make the problem numerically stable when using curl elements. The feature adding
this equation can only be applied when the magnetic field is discretized using curl
elements, that is, in 3D components, and 2D or 2D axisymmetric components when
the in-plane components are solved for.
• In magnetostatics, the magnetic flux density B does not appear in the first equation.
To allow modeling of magnetic materials and permanent magnets, it is necessary to
include B (with its constitutive relation) through the Magnetic Gauss’ Law node.
Em = ρ J
–1
Em = σ J
–1 –1
Em = ( ( ρ ) + j ωε 0 ε r ) J
–1
E m = ( σ + j ωε 0 ε r ) J
If the E-J characteristic is nonlinear, such in the case of a superconductor, the function
can be expressed explicitly or as a material property.
– ∇ ⋅ ( μ 0 ∇V m – μ 0 M ) = 0 (5-5)
The Magnetic Fields, No Currents Interface uses this equation for modeling of
magnetostatics in the absence of electric currents.
In planar 2D the dynamic formulations also involves the thickness d in the z direction:
– ∇ ⋅ d ( μ 0 ∇V m – μ 0 M ) = 0
For a discussion about the boundary element method used with The
Magnetic Fields, No Currents, Boundary Elements Interface, see Theory
for the Boundary Elements PDE in the COMSOL Multiphysics Reference
Manual.
The Magnetic and Electric Currents interface only supports the stationary
and frequency domain study types—that is, there is no transient
formulation available.
In this section:
• Magnetostatics Equations
• Frequency Domain Equations
Magnetostatics Equations
To derive the magnetostatics equations, start with Ampère’s law for static cases:
∇ × H = J. Define the potentials,
B = ∇×A
E = – ∇V
and the current
J = σv × B – σ ∇V + J e
–1
∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) + σ ∇V = J e
– ∇ ⋅ ( – σv × ( ∇ × A ) + σ∇V – J e ) = 0
–1
∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) + σ ∇V = J e
∇ × H = J = σ ( E + v × B ) + ( jωD + J e )
B = ∇×A
E = –∇V –j ωA
and combine them with the constitutive relationships B = μ0(H + M) and D = ε0E + P
to rewrite Ampère’s law as
2 –1
( jωσ – ω ε 0 )A + ∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) + ( σ + jωε 0 )∇V – jωP = J e
The equation of continuity is again obtained by taking the divergence of Ampère’s law.
It is the equation solved for the electric potential. Thus the following equations for V
and A apply:
2
– ∇ ⋅ ( ( jωσ – ω ε 0 )A – σv × ( ∇ × A ) + ( σ + jωε 0 )∇V – ( J e + jωP ) ) = 0
2 –1
( jωσ – ω ε 0 )A + ∇ × ( μ 0 ∇ × A – M ) – σv × ( ∇ × A ) + ( σ + jωε 0 )∇V – jωP = J e
The Electrical Circuit interface, found under the AC/DC branch ( ) when adding
interfaces, has the equations for modeling electrical circuits with or without
connections to a distributed fields model, solving for the voltages, currents and
charges associated with the circuit elements.
In this chapter:
321
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.
DEVICE NAMES
Each circuit component has an associated Device name, which is constructed from a
prefix identifying the type of the device and a string. The string can be specified in the
feature’s Settings window. The Device name is used to identify variables defined by the
component, and for the SPICE Import and Export functionality.
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 Node
The Ground Node ( ) feature adds a ground node with the default node number zero
to the electrical circuit. This is the default node in the Electrical Circuit interface. More
ground nodes can be added but those must have unique node numbers and are by
default given higher node numbers.
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.
DEVICE NAME
Enter a Device name for the ammeter. The name is also used in the probe.
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.
DEVICE NAME
Enter a Device name for the resistor. The prefix is R.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
DEVICE PARAMETERS
Enter the Resistance of 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.
DEVICE NAME
Enter a Device name for the inductor. The prefix is L.
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.
DEVICE NAME
Enter a Device name for the voltage source. The prefix is V.
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.
DEVICE NAME
Enter a Device name for the current source. The prefix is I.
DEVICE PARAMETERS
Enter the Source type that should be adapted to the selected study type. It can be
DC-source, AC-source, or a time-dependent Sine source or Pulse source. Depending on
the choice of source, also specify the following parameters:
• 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.
DEVICE PARAMETERS
Enter the voltage Gain. The resulting voltage is this number multiplied by the control
voltage.
DEVICE NAME
Enter a Device name for the voltage-controlled current source. The prefix is G.
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
Enter the source Gain (SI units: S). The resulting current is this number multiplied by
the control voltage. It represents the transconductance of the source.
DEVICE NAME
Enter a Device name for the current-controlled voltage source. The prefix is H.
DEVICE PARAMETERS
Enter the voltage Gain and select the Device whose current is taken as the control
current. The resulting voltage is this number multiplied by the control current through
the named Device (any two-pin device). Thus it formally has the unit of resistance.
DEVICE NAME
Enter a Device name for the current-controlled current source. The prefix is F.
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
Enter the current Gain and select the Device whose current is taken as the control
current. The resulting current is this number multiplied by the control current
through the Device.
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
Subcircuit Instance
The Subcircuit Instance ( ) feature represents an instance of a subcircuits defined by
a Subcircuit Definition feature.
DEVICE NAME
Enter a Device name for the subcircuit instance. The prefix is X.
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.
DEVICE NAME
Enter a Device name for the BJT. The prefix is Q.
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.
DEVICE NAME
Enter a Device name for the MOSFET. The prefix is M.
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.
M = k L1 L2
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
Chose 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.
DEVICE NAME
Enter a Device name for the diode. The prefix is D.
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.
DEVICE NAME
Enter a Device name for the External I vs. U node.
NODE CONNECTIONS
Specify the two Node names for the connecting nodes for the voltage source. The first
node represents the positive reference terminal.
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.
DEVICE NAME
Enter a Device name for the External U vs. I node.
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 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.
DEVICE NAME
Enter a Device name for the External I-terminal.
NODE CONNECTIONS
Set the Node name for the connecting node for the voltage assignment.
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.
• 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 6-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 bc v be I KF A I KR A
2 1 – ----------- – -----------
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
BF
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 –M Jx
1 – ---------
-
V Jx v bx < F C V Jx
C jbx
= AC Jx ×
v bx v bx ≥ F C V Jx
( 1 – F ) – 1 – MJx 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
–MJ
bd
--------
-
PB v bx < F C P B
C jbd = C BD ×
v bx v bx ≥ F C P B
( 1 – F ) – 1 – MJ 1 – F ( 1 + M ) + M -------
-
C C J J PB
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
VJ
k B T NOM
V T = ------------------------
q
The AC/DC Module license includes a physics interface found under the
Heat Transfer>Electromagnetic Heating branch ( ). This physics interface
combines magnetic fields with heat transfer for modeling of electromagnetic
heating such as induction heating. The Induction Heating Interface is discussed in
this chapter.
349
T he I nduc ti o n H eat i n g In t erfac e
The Induction Heating interface ( ) is used to model induction heating and eddy
current heating. This multiphysics interface adds a Magnetic Fields interface and a
Heat Transfer in Solids interface. The multiphysics couplings add the electromagnetic
power dissipation as a heat source, and the electromagnetic material properties can
depend on the temperature.
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
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 Magnetic Fields and Heat Transfer in Solids interfaces are added,
COMSOL Multiphysics adds an empty Multiphysics Couplings 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 Couplings node) or from the Physics toolbar, Multiphysics
menu.
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 Magnetic Fields Interface are listed in the
section Domain, Boundary, Point, and Pair Nodes for the Magnetic Fields
Interface.
• See The Heat Transfer Interfaces in the COMSOL Multiphysics Reference Manual
for information about the available physics features for heat transfer.
Electromagnetic Heating
The Electromagnetic Heating multiphysics coupling 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 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 in the model) is emh1.
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.
You can also select None from either list to uncouple the Electromagnetic Heating node
from a physics interface. If the physics interface is removed from the Model Builder, for
example Heat Transfer in Solids is deleted, then the Heat transfer list defaults to None as
there is nothing to couple to.
If a physics interface is deleted and then added to the model again, and in
order to reestablish 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 Approaches in the
COMSOL Multiphysics Reference Manual.
Glossary
355
Glossary of Terms
anisotropy Variation of material properties with direction.
constitutive relation The relation between the D and E fields and between the B and
H fields. These relations depend on the material properties.
electric dipole See electric currents point dipole, electrostatics point dipole, and
magnetic fields electric dipole.
electric currents point dipole Two equal and opposite point current sources +q and −q
separated a short distance d. The electric current dipole moment is given by p = qd,
where d is a vector going from −q to +q.
electrostatics point 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.
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.
magnetic fields electric dipole The limiting case when the length d of a current
filament carrying uniform current I approaches zero while maintaining the product
between I and d.
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| 359
change cross-section (node) 135 electric currents, shell 179
change shell thickness (node) 179 electrical circuits 327
change thickness (out-of-plane) (node) current sources, theory 199
135 current-controlled current source
charge conservation (node) 119, 144 (node) 330
charge conservation, piezoelectric current-controlled voltage source
(node) 135 (node) 329
charge relaxation theory 197 cylindrical coordinates 22, 65
circuit import, SPICE 337
D device models, electrical circuits 340
circular coil 86
dielectric shielding (node) 126
coenergy 54
dielectrics and perfect conductors 40
coil (node) 233
diode (node) 334
coil geometry (node) 240
diode transistor model 346
coil wires inlet/outlet (node) 243
dispersive materials 39
coil, single-turn 247
distributed capacitance (node) 131
collector node 331
distributed impedance (node) 160
common settings 22
documentation 23
conductive media 211
domain group (node) 252
conductor model (coil) 233
domain nodes
connected boundaries (node)
electric currents 151
coil 246–247
electrostatics 116
constitutive relations theory 33–34
magnetic and electric fields 300
contact impedance (node) 165
magnetic field formulation 264
continuum mechanics forces 45
magnetic fields 217
Cooper-Mikic-Yovanovich (CMY) corre-
magnetic fields, no currents 271
lation 208
drain node 332
Coulomb gauge 311
E edge current (node) 256
coupling, to the electrical circuits inter-
edge nodes
face 107
electric currents 151
curl elements 316
electric currents, shell 176
curl shape functions 231
electrostatics 116
current conservation (node)
magnetic and electric fields 300
electric currents 153
effective medium (node) 138
electric currents, shell 177
E-J characteristic 315
current conservation, ungauged formula-
elastic material stress tensor
tions 313
electromagnetics 51
current source (node)
moving bodies 54
electric currents 156
elastic pure conductor 51, 53
360 | I N D E X
electric coenergy 54 electromagnetic quantities 57
electric currents interface 149 electromagnetic sources, applying 67
theory 206 electromagnetic stress tensors 49
electric currents solid to shell via bound- electrostatic point dipole (node) 136
ary 195 electrostatics interface 114
electric currents, layered shell interface theory 196
181 electrostatics, boundary elements inter-
electric currents, shell interface 173 face 139
theory 211 emailing COMSOL 25
electric displacement field (node) 124 emitter node 331
electric field (node) 267 equation of continuity 197
electric field transformation (node) 291 error message, electrical circuits 107
electric fields theory 196 exporting
electric forces and torques 56 SPICE netlists 111
Electric Insulation 245 external current density (node)
electric insulation (node) electric currents 155
electric currents 157 magnetic fields 225
electric point dipole (node) external I vs. U (node) 334
electric currents 171 external I-terminal (node) 336
magnetic fields 231 external magnetic flux density (node) 275
electric point dipole (on axis) (node) 171 external magnetic vector potential
electric potential (node) 122, 144 (node) 256
electric scalar-scalar potential coupling external surface charge accumulation
(node) 147 (node) 123
electric shielding (node) external U vs. I (node) 335
electric currents 164
F Faraday’s law 315
electric currents, shell 180
Faraday’s law (node) 265
electrical circuit interface 322
field variables in 2D 69
theory 339
file, Touchstone 96
electrical circuits
fixed current, ports 94
modeling techniques 106
floating potential (node)
electrical conductivity, porous media 204
electric currents 157
electrical contact (node) 167
electrostatics 129, 145
electrical size 63
single-turn coil 251
electromagnetic energy, theory 36
fluid saturation 137
electromagnetic forces and torques 72
force calculation (node) 120
electromagnetic forces, calculating 45
forced voltage, port 94
electromagnetic heating (node) 353
forces
INDEX| 361
calculating 72 SPICE netlists 110, 337
continuum mechanics and 45 induction heating interface 350
elastic solids, and 47 inductor (node) 326
electromagnetic, calculating 45 inhomogeneous materials 38
in moving objects 52 initial values (node)
stationary fields 48 electric currents 155
torque 48 electric currents, shell 179
formation factor 205 electrostatics 120
frequency domain study 207 magnetic and electric fields 302
force calculation and 121 magnetic field formulation 266
theory, electric currents 199 magnetic fields 225
magnetic fields, no currents 273
G Galilei invariants and transformations 52
magnetic fields, no currents, boundary
gap feed (node) 250
elements 282
gate node 332
rotating machinery, magnetic 291
gauge fixing 312
input (node), coil 242
gauge fixing for A-field (node) 232
input (node), coil domains 244
gauge transformation 311
interface current source 194
Gauss’ law and charge relaxation theory
internet resources 23
197
Gauss’ law equation 200 J Jiles-Atherton model 42
general elastic materials Jiles-Atherton, hysteresis 42
electromagnetics 51
K Kirchhoff’s circuit laws 339
moving bodies 54
knowledge base, COMSOL 26
geometry, simplifying 64
L layer current source 193
ground (node) 122, 145
layered shell 189
ground (node), single-turn coil domain
line charge (node) 131
251
line charge (on axis) (node) 132
ground node (node) 324
line charge (out-of-plane) (node) 133
H harmonic perturbation (node), global
line current (on axis) (node) 231
240
line current (out-of-plane) (node) 230
Helmholtz’s theorem 311
line current source (node) 169
homogenized multi-turn (coil) 234
line current source (on axis) (node) 169
homogenized shell 191
linear coil 86
homogenized shell connection 191
Lorentz force equation 38
I impedance boundary condition (node) Lorentz forces, calculating 72
256 Lorentz term 37
importing lumped element (node) 255
362 | I N D E X
lumped parameters magnetic potential (node) 229
converting 97 magnetic scalar potential (node) 274
energy method 95 magnetic fields, no currents, boundary
fixed current 94 elements 283
Ohm’s law and 93 magnetic scalar potential discontinuity
studying 96 (node)
lumped port (node) 252 magnetic fields, no currents 276
lumped ports 101 magnetic scalar-scalar potential coupling
(node)
M magnetic and electric fields interface 294
magnetic fields, no currents, boundary
theory 319
elements 283
magnetic coenergy 54
magnetic shielding (node)
magnetic continuity (node) 304
magnetic and electric fields 303
magnetic field (node) 228
magnetic currents 260
magnetic field formulation interface 262
magnetic fields, no currents 275
theory 315
magnetic vector-scalar potential coupling
magnetic fields interface 214
(node)
theory 308
magnetic fields, no currents, boundary
magnetic fields, no currents interface 269
elements 285
theory 318
magnetostatics, no currents, boundary
magnetic flux conservation (node) 273
elements interface (node)
magnetic fields, no currents, boundary
magnetic fields, no currents, boundary
elements 281
elements 278
magnetic flux density (node) 274
material properties
magnetic fields, no currents, boundary
electromagnetics 38
elements 283
Maxwell capacitance 97
magnetic forces and torques 55
Maxwell stress tensor, calculating 72
magnetic Gauss’ law 316
Maxwell’s equations
magnetic Gauss’ law (node) 266
dielectrics 40
magnetic insulation (node)
electrical circuits and 339
magnetic and electric fields 303
quasi-static approximation 310
magnetic currents 227
theory 32
magnetic field formulation 267
mesh resolution 70
magnetic fields, no currents 274
method of virtual work 54
magnetic fields, no currents, boundary
Mikic elastic correlation 209
elements 282
mixed formulation boundary (node) 293
magnetic point dipole (node) 259
MPH-files 25
magnetic point dipole (on axis) (node)
multiphysics couplings
260
INDEX| 363
induction heating 350 physics interfaces, common settings 22
multi-turn, homogenized (coil) 234 PMC. see perfect magnetic conductor.
mutual capacitance matrix 98 PNP BJT (node) 331
mutual inductance (node) 333 point charge (node) 133
point charge (on axis) (node) 134
N n-Channel MOS transistor 332, 343
point current source (node) 170
n-Channel MOSFET (node) 332
point current source (on axis) (node)
netlists, SPICE 110, 337
170
nodes, common settings 22
point dipoles, magnetic 259
nonlinear materials 39
point nodes
normal current density (node)
electric currents 151
electric currents 159
electric currents, shell 176
electric currents, shell 180
electrostatics 116
NPN bipolar junction transistor 340
magnetic and electric fields 300
NPN BJT (node) 331
magnetic field formulation 264
numeric coil 87
magnetic fields 217
O Ohm’s law and charge relaxation theory
magnetic fields, no currents 271
197
ports, lumped 101
output (node), coil 242
potentials, scalar and magnetic 35
output (node), coil domains 245
power law, porous media
P pair contact impedance (node) 165 conductivity 202
pair electrical contact (node) 167 permeability 204
pair nodes permittivity 203
electric currents 151 Poynting’s theorem 37
electric currents, shell 176 predefined couplings, electrical circuits
electrostatics 116 107
magnetic field formulation 264 prescribed rotation (node) 292
magnetic fields 217 prescribed rotational velocity (node) 292
magnetic fields, no currents 271 principle of virtual displacement 54
pair thin low permeability gap (node) 277 pure conductor stress tensor 51, 53
p-Channel MOS transistor 332
Q quasi-static approximation 37, 310
p-Channel MOSFET (node) 332
R reciprocal permeability, volume average
PDE formulations 35
204
perfect conductors and dielectrics 40
reciprocal permittivity, volume average
perfect magnetic conductor (node) 230
203
periodic condition (node) 124
resistance matrix 98
phasors 40
resistor (node) 325
physics interface guide 17
RLC coil group (node) 305
364 | I N D E X
robust geometry analysis method (coil) physics interface availability 21
86 selecting 68
rotating machinery, magnetic interface subcircuit definition (node) 330
287 subcircuit instance (node) 331
surface charge density (node) 123, 145
S saturation coefficient 205
surface current (node)
saturation exponent 137
magnetic and electric fields 303
scattering parameters. see S-parameters
magnetic currents 228
sector symmetry (node) 168
surface magnetic current density (node)
selecting
229, 268
solvers 71
sweeps and lumped parameters 96
space dimensions 64
symmetry specification (coil) 86
study types 68
semiconductor device models 340 T technical support, COMSOL 25
shape functions 231 terminal (node)
single conductor (coil) 233 electric currents 162
single-turn coil (node) 247 electrostatics 127
skin effects 70 electrostatics, boundary elements 146
solver settings 71 theory
solving models 64 constitutive relations 33–34
source node 332 electric currents interface 206
space charge density (node) 120 electric currents, shell 211
space dimension, solving models and 64 electric fields 196
space dimensions, selecting 64 electrical circuit interface 339
S-parameter calculations electromagnetics 32
electric field, and 104 electrostatics interface 196
SPICE lumped ports 101
exporting 111 magnetic and electric fields 310, 319
SPICE capacitance matrix 97 magnetic field formulation 315
SPICE netlists 110, 337 magnetic fields interface 308
standard settings 22 magnetic fields, no currents 318
stationary fields, forces 48 thin low permeability gap (node)
stationary source sweep (node) 99 magnetic fields 259
stress tensors 53 magnetic fields, no currents 277
electromagnetics 51 thin low permittivity gap (node) 125
moving bodies 54 thin shells, conductive media 211
study types time dependent study 68, 199, 207
electric currents 207 torque forces 48
electromagnetic theory 35 torques, calculating 72
INDEX| 365
Touchstone and port sweeps 96
transformations, Galilei 52
transformer (node) 333
transient study 68
transition boundary condition (node)
258
W wavelength 71
websites, COMSOL 26
366 | I N D E X