MEMSModuleUsersGuide
MEMSModuleUsersGuide
User’s Guide
MEMS Module User’s Guide
© 1998–2023 COMSOL
Protected by patents listed on www.comsol.com/patents, or see Help>About COMSOL Multiphysics on the
File menu in the COMSOL Desktop for less detailed lists of U.S. Patents that may apply. Patents pending.
This Documentation and the Programs described herein are furnished under the COMSOL Software License
Agreement (www.comsol.com/sla) and may be used or copied only under the terms of the license
agreement.
COMSOL, the COMSOL logo, COMSOL Multiphysics, COMSOL Desktop, COMSOL Compiler,
COMSOL Server, and LiveLink are either registered trademarks or trademarks of COMSOL AB. All other
trademarks are the property of their respective owners, and COMSOL AB and its subsidiaries and products
are not affiliated with, endorsed by, sponsored by, or supported by those trademark owners. For a list of such
trademark owners, see www.comsol.com/trademarks.
Version: COMSOL 6.2
Contact Information
Visit the Contact COMSOL page at www.comsol.com/contact to submit general inquiries
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 on the COMSOL Access
page at www.comsol.com/support/case. Other useful links include:
Chapter 1: Introduction
Structural Mechanics 45
Damping and Losses . . . . . . . . . . . . . . . . . . . . . 46
Electrostatics 47
Electrostatics . . . . . . . . . . . . . . . . . . . . . . . . 47
Electric Currents in Conductive Media . . . . . . . . . . . . . . 48
Lumped Parameters. . . . . . . . . . . . . . . . . . . . . . 49
CONTENTS |3
Piezoresistivity 54
Microfluidics 60
Fluid-Structure Interaction . . . . . . . . . . . . . . . . . . . 60
Thin-Film Damping . . . . . . . . . . . . . . . . . . . . . . 60
Reference for Thin-Film Damping . . . . . . . . . . . . . . . . 62
4 | CONTENTS
Force Calculation. . . . . . . . . . . . . . . . . . . . . . . 80
Zero Charge . . . . . . . . . . . . . . . . . . . . . . . . 81
Ground . . . . . . . . . . . . . . . . . . . . . . . . . . 82
Electric Potential . . . . . . . . . . . . . . . . . . . . . . . 82
Surface Charge Density . . . . . . . . . . . . . . . . . . . . 83
External Surface Charge Accumulation . . . . . . . . . . . . . . 83
Symmetry Plane (for Electric Field) . . . . . . . . . . . . . . . . 84
Electric Displacement Field . . . . . . . . . . . . . . . . . . . 84
Periodic Condition . . . . . . . . . . . . . . . . . . . . . . 85
Thin Low Permittivity Gap . . . . . . . . . . . . . . . . . . . 86
Dielectric Shielding . . . . . . . . . . . . . . . . . . . . . . 87
Terminal . . . . . . . . . . . . . . . . . . . . . . . . . . 87
Floating Potential . . . . . . . . . . . . . . . . . . . . . . . 89
Distributed Capacitance . . . . . . . . . . . . . . . . . . . . 91
Line Charge . . . . . . . . . . . . . . . . . . . . . . . . . 91
Line Charge (on Axis) . . . . . . . . . . . . . . . . . . . . . 92
Line Charge (Out-of-Plane) . . . . . . . . . . . . . . . . . . . 92
Point Charge . . . . . . . . . . . . . . . . . . . . . . . . 93
Point Charge (on Axis) . . . . . . . . . . . . . . . . . . . . 94
Change Cross Section . . . . . . . . . . . . . . . . . . . . . 94
Change Thickness (Out-of-Plane). . . . . . . . . . . . . . . . . 95
Charge Conservation, Ferroelectric. . . . . . . . . . . . . . . . 95
Charge Conservation, Piezoelectric . . . . . . . . . . . . . . . . 95
Dispersion . . . . . . . . . . . . . . . . . . . . . . . . . 96
Electrostatic Point Dipole . . . . . . . . . . . . . . . . . . . 99
CONTENTS |5
The Electric Currents Interface 110
Domain, Boundary, Edge, Point, and Pair Nodes for the Electric
Currents Interface . . . . . . . . . . . . . . . . . . . . 113
Current Conservation . . . . . . . . . . . . . . . . . . . . 115
Dispersion . . . . . . . . . . . . . . . . . . . . . . . . 118
Initial Values . . . . . . . . . . . . . . . . . . . . . . . 120
External Current Density . . . . . . . . . . . . . . . . . . 121
Current Source . . . . . . . . . . . . . . . . . . . . . . 121
Electric Insulation . . . . . . . . . . . . . . . . . . . . . 122
Floating Potential . . . . . . . . . . . . . . . . . . . . . . 122
Boundary Current Source . . . . . . . . . . . . . . . . . . 124
Normal Current Density . . . . . . . . . . . . . . . . . . . 124
Distributed Impedance. . . . . . . . . . . . . . . . . . . . 125
Terminal . . . . . . . . . . . . . . . . . . . . . . . . . 127
Electric Shielding . . . . . . . . . . . . . . . . . . . . . . 130
Contact Impedance . . . . . . . . . . . . . . . . . . . . . 130
Sector Symmetry . . . . . . . . . . . . . . . . . . . . . . 132
Line Current Source . . . . . . . . . . . . . . . . . . . . 133
Line Current Source (on Axis). . . . . . . . . . . . . . . . . 134
Point Current Source . . . . . . . . . . . . . . . . . . . . 134
Point Current Source (on Axis) . . . . . . . . . . . . . . . . 135
Electric Point Dipole . . . . . . . . . . . . . . . . . . . . 135
Electric Point Dipole (on Axis) . . . . . . . . . . . . . . . . 136
Current Conservation, Piezoresistive . . . . . . . . . . . . . . 136
6 | CONTENTS
Switch . . . . . . . . . . . . . . . . . . . . . . . . . . 145
Subcircuit Definition . . . . . . . . . . . . . . . . . . . . 146
Subcircuit Instance . . . . . . . . . . . . . . . . . . . . . 147
NPN BJT and PNP BJT. . . . . . . . . . . . . . . . . . . . 147
n-Channel MOSFET and p-Channel MOSFET . . . . . . . . . . . 148
Mutual Inductance . . . . . . . . . . . . . . . . . . . . . 148
Transformer . . . . . . . . . . . . . . . . . . . . . . . 149
Diode . . . . . . . . . . . . . . . . . . . . . . . . . . 149
External I vs. U . . . . . . . . . . . . . . . . . . . . . . 150
External U vs. I . . . . . . . . . . . . . . . . . . . . . . 151
External I-Terminal . . . . . . . . . . . . . . . . . . . . . 152
SPICE Circuit Import . . . . . . . . . . . . . . . . . . . . 153
SPICE Circuit Export . . . . . . . . . . . . . . . . . . . . 153
CONTENTS |7
Chapter 4: Fluid Flow Interfaces
8 | CONTENTS
The Piezoresistivity, Shell Interface 226
Piezoresistivity, Shell . . . . . . . . . . . . . . . . . . . . 228
Chapter 7: Glossary
Index 249
CONTENTS |9
10 | C O N T E N T S
1
Introduction
This user’s guide describes the MEMS Module. The MEMS Module User’s Guide
introduces the modeling stages in COMSOL Multiphysics® and this module and
serves as a reference for more advanced modeling techniques and details about the
physics interfaces.
In this chapter:
11
About the MEMS Module
In this section:
• About MEMS
• About the MEMS Module
• The MEMS Module Physics Interface Guide
• Common Physics Interface and Feature Settings and Nodes
• Coupling to Other Physics Interfaces
• Where Do I Access the Documentation and Application Libraries?
About MEMS
The field of MEMS evolved as engineers and scientists explored new avenues to make
use of the fabrication technologies developed by the microelectronics industry. These
technologies enabled complex micron and submicron structures to be integrated with
electronic systems and batch fabricated at low cost. Mechanical devices fabricated using
these technologies have become known as microelectromechanical systems (MEMS)
or alternatively as microsystems. A proper description of these devices usually requires
multiple physical effects to be incorporated.
12 | CHAPTER 1: INTRODUCTION
common with the widespread adoption of FBAR (thin–film bulk acoustic resonator)
technology in mobile communication applications. Although there are few
commercially important devices that employ thermal actuation (largely due to the high
power consumption of such devices), thermal phenomena, such as packaging-induced
thermal stresses, are often important in the design of MEMS devices. Finally, for
resonant MEMS devices such as accelerometers, damping phenomena are important
and there is a need for detailed modeling of phenomena such as thin-film damping.
Very high quality factor resonators are often used in MEMS timing applications, and
for these devices the quality factor can be limited by thermoelasticity: the intrinsic
coupling of the thermal and structural domains that exists as a result of the laws of
thermodynamics.
Microfluidic systems are often included within the field of MEMS and represent a
growing market for MEMS devices. For modeling microfluidics devices, the dedicated
Microfluidics Module contains tailored physics interfaces for the simulation of
microfluidic devices, including creeping flow, multiphase flow, electrokinetic effects,
and slip flow.
AC/DC
14 | CHAPTER 1: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
Fluid Flow
Single-Phase Flow
16 | CHAPTER 1: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
Structural Mechanics
Thermal–Structure Interaction
18 | CHAPTER 1: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
Piezoresistivity
20 | CHAPTER 1: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
The first column of the table lists phenomena, couplings, and devices that are often
associated with the word electromechanical in the narrow and literal meaning of
MEMS. The devices in this category are usually various kinds of actuators and sensors.
The microfluidic devices, although using some of the same manufacturing and
miniaturizing techniques, form a totally different application area. The
Fluid-Structure Interaction column lists various techniques and phenomena that are
useful for both electromechanical and microfluidic applications, where movement and
deformation of solids are of concern.
This table, however, shows only the tip of the iceberg — our view of the most
important applications where you can use the MEMS Module. In your hands, the
multiphysics combinations and applications are unlimited.
TABLE 1-1: EXAMPLES OF COUPLED PHENOMENA AND DEVICES YOU CAN MODEL USING THE MEMS MODULE.
ELECTROMECHANICAL FLUID-STRUCTURE
INTERACTION
22 | CHAPTER 1: INTRODUCTION
THE MEMS MODULE STUDY CAPABILITIES BY PHYSICS INTERFACE
This section lists the physics interfaces, the physical quantities they solve for, and the
standard abbreviation each one uses. The physical quantities in these physics interfaces
are:
Table 1-2 lists the physics interfaces specific to MEMS modeling and this module.
TABLE 1-2: MEMS MODULE INTERFACE DEPENDENT VARIABLES AND PRESET STUDY OPTIONS.
SMALL-SIGNAL ANALYSIS,
PRESTRESSED ANALYSIS,
PRESTRESSED ANALYSIS,
FREQUENCY. DOMAIN
FREQUENCY-DOMAIN
FREQUENCY DOMAIN
FREQUENCY DOMAIN
EIGENFREQUENCY
EIGENFREQUENCY
TIME DEPENDENT
STATIONARY
MODAL
AC/DC
Electrostatics es V
Electric Currents ec V
Electrical Circuit cir none
FLUID FLOW
TIME-DEPENDENT MODAL
SMALL-SIGNAL ANALYSIS,
PRESTRESSED ANALYSIS,
PRESTRESSED ANALYSIS,
FREQUENCY. DOMAIN
FREQUENCY-DOMAIN
FREQUENCY DOMAIN
FREQUENCY DOMAIN
EIGENFREQUENCY
EIGENFREQUENCY
TIME DEPENDENT
STATIONARY
MODAL
Thin-Film Flow tff pf
Thin-Film Flow, tff pf
Domain
Fluid-Structure — u, v, w
Interaction p
Fluid-Structure — u, v, w
Interaction, Fixed p
Geometry
STRUCTURAL MECHANICS
24 | CHAPTER 1: INTRODUCTION
TABLE 1-2: MEMS MODULE INTERFACE DEPENDENT VARIABLES AND PRESET STUDY OPTIONS.
TIME-DEPENDENT MODAL
SMALL-SIGNAL ANALYSIS,
PRESTRESSED ANALYSIS,
PRESTRESSED ANALYSIS,
FREQUENCY. DOMAIN
FREQUENCY-DOMAIN
FREQUENCY DOMAIN
FREQUENCY DOMAIN
EIGENFREQUENCY
EIGENFREQUENCY
TIME DEPENDENT
STATIONARY
MODAL
Piezoresistivity, pzrb u, v, w
Boundary V
Currents
Piezoresistivity, pzrs u, v, w
Shell V, ar
• In the Model Builder or Physics Builder, click a node or window and then
press F1.
• In the main toolbar, click the Help ( ) button.
• From the main menu, select Help>Help.
• Press Ctrl+F1.
• In the main toolbar, click the Documentation ( ) button.
• From the main menu, select Help>Documentation.
26 | CHAPTER 1: INTRODUCTION
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.
MEMS MODELING
Creating and Analyzing MEMS Models familiarizes you with modeling procedures
useful when working with this module and MEMS applications such as microactuators
and microsensors. Because the module is fully integrated with COMSOL Multiphysics,
the design and analysis process is similar to the one you use in the base environment.
The MEMS-specific modeling aspects should prove valuable during the modeling
phase. There is also a short section (with model examples) about the materials
databases available with this module.
28 | CHAPTER 1: INTRODUCTION
The Thin-Film Flow Interfaces (in the CFD Module User’s Guide) can be added
either singularly or in combination with other physics interfaces such as mass and
energy transfer. The Thin-Film Flow interface is for 2D, 2D axisymmetric and 3D
models, and the Thin-Film Flow, Domain interface is for 2D models.
The Fluid–Solid Interaction Interface (in the Structural Mechanics Module User’s
Guide) is an interface found under the Fluid Flow branch. It has the equations and
physics for fluid-structure interaction, solving for the displacements, fluid velocity, and
fluid pressure. You have access to all functionality from the Solid Mechanics interface
and Single-Phase Flow interfaces for modeling the solid and fluid domains. From the
Laminar Flow (the default) and Solid Mechanics submenus on the solid and boundary
level you can access materials, sources, loads, and boundary conditions for the
individual physics interfaces.
The Piezoresistivity, Domain Currents Interface is appropriate for situations when the
thickness of the conducting and piezoresistive layers are both resolved by the mesh.
The Piezoresistivity, Boundary Currents Interface should be used when the thicknesses
of the conducting and piezoresistive layers are much smaller than those of the
structural layers (this is often the case in practice).
The Piezoresistivity, Shell Interface is used when the structural layer is thin enough to
be treated by the shell interface, but the conducting and piezoresistive layers are still
much thinner than the structural layers.
30 | CHAPTER 1: INTRODUCTION
2
This chapter gives an overview of the physics interfaces and features available for
modeling MEMS devices in COMSOL Multiphysics®.
In this chapter:
| 31
Physics and Scaling for MEMS Devices
Microelectromechanical systems (MEMS) is a term that covers a broad range of
manufactured mechanical components characterized by a small size. Many of these
devices are fabricated with micron-scale dimension and represent mixed
mechanical/electrical systems, although the term now encompasses systems that do
not necessarily contain electrical parts, such as microfluidic devices. Thousands of
MEMS devices can be batch fabricated on a single wafer, using technologies derived
from the semiconductor industry. As a result of the smaller size and reduced
manufacturing costs, MEMS-based equivalents of macroscopic devices such as
accelerometers have been able to penetrate new markets. Additionally, MEMS
technology has enabled entirely new devices to be created, such as the micromirror
arrays that enable Texas Instruments’ DLP projectors.
The smaller size of MEMS devices has important effects on the physical processes that
govern their operation (Ref. 1). In general, as the length scale (L) of the device is
reduced, the scaling of a physical effect with respect to L determines its relative
importance (see Table 2-1). The inertial force required to produce a fixed acceleration
of a solid body scales volumetrically as L3. The scaling of other forces in comparison
to this inertial force has important consequences for MEMS devices. For example, the
effective spring constant for a body scales as L1. The spring stiffness therefore
decreases much more slowly than the system mass as the size of the system is reduced,
resulting in higher resonant frequencies for smaller devices (resonant frequency scales
as L1). This means that micromechanical systems typically have higher operating
frequencies and faster response times than macroscopic systems.
Electrostatic forces scale favorably as the device dimensions are reduced (for example,
the force between parallel plates with a fixed applied voltage scales as L0) (see
Table 2-1). Additionally, electrostatic actuators consume no DC power and can be
Electrostatically actuated and sensed MEMS devices dominate the market for inertial
sensors (usually employed in the automobile and consumer electronics industries).
Magnetic MEMS devices can also be modeled in COMSOL Multiphysics, but since
Piezoelectric forces also scale well as the device dimension is reduced (the force
produced by a constant applied voltage scales as L1). Furthermore, piezoelectric
sensors and actuators are predominantly linear and do not consume DC power in
operation. Piezoelectrics are more difficult to integrate with standard semiconductor
processes, but significant progress has been made with commercial successes in the
market (for example, Avago Technologies (formerly Agilent/HP) FBAR filters). High
frequency FBAR resonators fabricated from aluminum nitride thin films on silicon
wafers are now widely used as filters used in consumer devices such as mobile phones.
Thermal forces scale as L2, assuming that the forces are generated by a fixed
temperature change. This scaling is still favorable in comparison to inertial forces, and
the thermal time scale also scales well (as L2), making thermal actuators faster on the
microscale (although thermal actuators are typically slower than capacitive or
piezoelectric actuators). Thermal actuators are also easy to integrate with
semiconductor processes although they usually consume large amounts of power and
thus have had a limited commercial applicability. Thermal effects play an important
role in the manufacture of many commercial MEMS technologies with thermal stresses
in deposited thin films being critical for many applications.
The interaction of MEMS devices with a fluid leads to a range of microfluidic devices.
COMSOL also has the Microfluidics Module, which is geared to model microfluidic
devices, but some functionality to model fluids is included with the MEMS Module.
In particular, the interaction of moving components with a fluid can be modeled using
the Fluid-Structure Interaction Interface. In addition to modeling microfluidic
actuators, this multiphysics interface can be used to model the damping of moving
Microfluidics
This discussion has focused on actuators; however, many of the actuators discussed can
also be used as sensors. Many MEMS devices (for example, gyroscopes and resonant
sensors) require both actuators and sensors to be integrated into the same device.
Electrostatic, magnetic, piezoelectric, and even thermal devices have all been employed
as MEMS sensors. Another commonly used sensing mechanism found in commercial
MEMS devices is that of piezoresistivity. The piezoresistive effect refers to the change
in a material’s conductivity that occurs in response to an applied stress. Piezoresistors
are relatively easy to manufacture using common semiconductor processes and are
inherently shielded from electromagnetic interference, and it is usually straightforward
to implement electronics to interface with them. Pressure sensors based on the
piezoresistive effect were some of the first MEMS devices to be mass produced.
Piezoresistive devices do, however, usually consume more power and generate more
electrical noise than capacitive sensors, which are displacing piezoresistive devices in
some applications.
Piezoresistivity
Electrical Circuits
Many MEMS devices include accompanying electrical circuits, either integrated on the
same chip or, increasingly, through connections to an adjacent chip that is packaged
with the device. COMSOL Multiphysics includes tools to connect physical finite
ꞏꞏ + Duꞏ + Ku = F
Mu (2-1)
where u is the displacement vector (size: n-by-1), K is the stiffness matrix (size:
n-by-n), D is the damping matrix (size: n-by-n), and M is the mass matrix (size:
n-by-n). In the frequency domain the problem takes the form
2
– M u 0 + iDu 0 + Ku 0 = F
where u = u0eit.
2
Kû i = i Mû i (2-2)
These eigenvectors are orthogonal with respect to M, as can be seen from the
following construction, based on Equation 2-2:
T 2 T
û j Kû i – i û j Mû i = 0
T 2 T
û i Kû j – j û i Mû j = 0
T 2 T T T 2 T
û i Kû j – j û i Mû j = û j Kû i – j û j Mû i = 0
using the results MT M, and KT K, which apply for physical reasons. Subtracting
the first equation from the third gives:
2 2 T
– j – i û j Mû i = 0
T
û j Mû i = 0 i j , i j (2-3)
T
û j Kû i = 0 i j , i j (2-4)
Next the following n-by-n matrix is constructed, with columns taken from the n
eigenvectors:
U = û 1 , û 2 û n
From Equation 2-3 it is clear that this is a diagonal matrix. Similarly from Equation 2-4
it is clear that UTKU is also diagonal.
From the properties of the eigenvectors it is possible to expand any function in terms
of the eigenvectors. Thus the displacement u can be written as:
u = ai ûi
i=1
u = Ua (2-5)
Now consider the original equation: Equation 2-1. First substitute for u using
Equation 2-5. Then transform the equation to the modal coordinate system by
premultiplying by UT. This gives:
T ꞏꞏ + U T DUaꞏ + U T KUa = U T F
U MUa (2-6)
It has already been established that the matrices UTMU and UTKU are diagonal, and
frequently a damping model is chosen that results in a diagonal damping matrix. For
example, in Rayleigh damping D M + K, where and are constants.
Alternatively, a damping ratio, i, is assigned to each mode.
so that
T 2
U KU = diag i
where diag(i2) is the diagonal matrix with diagonal elements i2. Similarly, if
damping ratios for each mode are defined, the damping matrix can be expressed in the
form
T
U DU = diag 2 i i
Thus if mass matrix scaling is used Equation 2-6 takes the form
ꞏꞏ + diag c
diag m eff , i a ꞏ T
eff , i a + diag k eff , i a = U F (2-8)
Here meff,i is the effective mass of the ith mode, ceff,i = 2meff,iii is the effective
damping parameter for the mode, and keff,i is the effective spring constant. Each
element of the vector UTF gives the force component that acts on each of the
respective modes.
For electromechanical systems it is common to express the force term, UTF, in terms
of a driving voltage, V, applied to a terminal in the system. Given this:
T
U F = V
ꞏꞏ + diag c
diag m eff , i a ꞏ
eff , i a + diag k eff , i a = V (2-9)
ꞏ
i mot , i = i a i
For each mode, the terms on the left of Equation 2-10 can be identified with a resistor
capacitor and inductor, such that the mechanical system behaves like a set of parallel
series LCR circuits. The equivalent modal resistances, capacitances, and inductances
are given by:
m eff , i
L i = --------------
-
2
i
c eff , i
R i = ------------
-
2
i
2
i
C i = ------------
-
k eff , i
For a given mode, these quantities are frequently referred to as the motional
inductance, resistance, and capacitance.
REDUCED-ORDER MODELS
The preceding discussion did not consider how to reduce the number of degrees of
freedom in the system. For systems in which the vector UTF has only a few significant
components (for example, components i = 1, …, m where m « n) the following
approximation can be made:
u ai ûi
i=1
where U' is now an m-by-n and a' is a vector with m components. The equation
system in modal coordinates now takes the form
T ꞏꞏ + U' T DU'a'
U' MU'a' ꞏ + U' T KU'a' = U' T F (2-11)
The matrices U'TMU', U'TDU', and U'TKU' now have dimensions m-by-m. Similarly
the vector U'TF has m components. This results in a significant reduction in the
system complexity.
2
– i M û r ,i + i i Dû r ,i + Kû r ,i = F
u = U' r a'
where U'r is the n-by-m matrix containing the right eigenvectors chosen for the modal
analysis. Once again a' is a vector with m components. The system in modal
coordinates takes the form
where U'l is the n-by-m matrix containing the left eigenvectors chosen for the modal
analysis.
The matrices U'lTMU'r, U'lTDU'r, and U'lTKU'r are no longer necessarily diagonal.
The modal solver accepts any linearly independent set of vectors to project the solution
vector and equations onto and constructs the reduced-order system accordingly.
After the model has solved, right-click the Results>Derived Values node and select
System Matrices. In the output section choose the Matrix to display in the list. The mass
matrix corresponds to the matrix U'lTMU'r the stiffness matrix corresponds to
U'lTKU'r, and the damping matrix corresponds to U'lTDU'r. The vector U'lTF is
available as the load vector. These matrices are given in a format that respects the
normalization of the preceding Eigenvalue Solver. To change this select the Eigenvalue
Solver node, and change the Scaling of Eigenvectors setting under the Output section.
Use the Max setting if an equivalent Mass-Spring-Damper system is required, in which
case the modal amplitude corresponds to the maximum displacement of the mode.
To convert this into an equivalent circuit representation, see the discussion in the
section Eigenvalue scaling and Equivalent Electrical and Mechanical Systems.
Use the ECAD Import Module to import mask data into COMSOL
Multiphysics. See the ECAD Import Module User’s Guide or go to
https://www.comsol.com/ecad-import-module for more information
about this and other COMSOL products.
Many important aspects of the creation of structural mechanics models are covered in
the chapter Structural Mechanics Modeling in the Structural Mechanics Module
User’s Guide. The following sections are of special importance for MEMS modeling;
STRUCTURAL MECHANICS | 45
Damping and Losses
For frequency domain (eigenfrequency or frequency domain studies) and transient
simulations it is often necessary to include damping in the model. This is particularly
important for resonant MEMS devices.
• Modeling losses through anchor structures in the time domain using a spring
damper system.
In this section:
• Electrostatics
• Electric Currents in Conductive Media
• Lumped Parameters
Electrostatics
The Electrostatics Interface solves for the electric potential V using the equation
– V =
For most MEMS devices the details of the charge distribution, , inside material are
unimportant and the charge usually resides only on surfaces held at fixed potentials.
The Electric Potential, Ground, and Terminal boundary conditions are usually used to
fix the potential on the surfaces of a MEMS device. At symmetry boundaries the Zero
Charge boundary condition applies. To represent infinite regions adjacent to the
model a common approach is to use infinite elements. Note that it is possible to model
electric circuits associated with a terminal, as described in the Connecting to Electrical
Circuits section.
ELECTROSTATIC FORCES
The Maxwell stress tensor is a general method for calculating electromagnetic forces
and torques. Force and torque calculations using Maxwell’s stress tensor are available
in the Electrostatics interface. In electrostatics the force is calculated by integrating
1 T
n 1 T 2 = – --- n 1 E D + n 1 E D (2-12)
2
ELECTROSTATICS | 47
on the surface of the object that the force acts on where E is the electric field, D the
electric displacement, and n1 the outward normal from the object.
In practice, electric fields are modeled not only in the air but also inside the solid
materials. By entering the electromagnetic force variable name into the domain
settings, the force calculations are activated, and the physics interface automatically
generates new domain and boundary variables.
CAPACITANCE CALCULATIONS
The capacitance of a system can be computed using the functionality built into the
Terminal boundary condition. This is discussed in more detail in the section Lumped
Parameters.
e
– V – J = Q j
Conductors in MEMS designs are often thin sheets. For cases in which large aspect
ratio differences occur it is often convenient to represent thin layers using a planar
approximation.
Lumped Parameters
Lumped parameters are matrices describing electromagnetic properties such as
resistance, capacitance, and inductance. In the time-harmonic case the lumped
parameter matrix is either an impedance matrix or an admittance matrix depending on
how the model is excited (current or voltage). In a static calculation only the resistive,
capacitive, or inductive part of the lumped parameter matrix is obtained.
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 2-1.
V1 V3
Ground
V2 V4
Figure 2-1: A five-electrode system with 4 terminals and one ground electrode.
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
ELECTROSTATICS | 49
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:
V1 Z 11 Z 12 Z 13 Z 14 I 1
V2 Z 21 Z 22 Z 23 Z 24 I 2
=
V3 Z 31 Z 32 Z 33 Z 34 I 3
V4 Z 41 Z 42 Z 43 Z 44 I 4
–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 ji
Vi
C ii = ------2- W e d Vj =
Vi j = i
0 k i j
1 Vi Vj
C ij = ------------ W e d – --- ------ C + ------ C jj
1
Vi Vj 2 V j ii V i
Vk = Vi
k = i
Vj k = j
S-PARAMETERS
Scattering parameters (or S-parameters) are complex-valued, frequency-dependent
matrices describing the transmission and reflection of electromagnetic waves at
different ports of devices like filters, antennas, waveguide transitions, and transmission
lines. S-parameters originate from transmission-line theory and are defined in terms of
transmitted and reflected voltage waves. All ports are assumed to be connected to
matched loads/feeds, that is, there is no reflection directly at a port.
S 11 S 12 . . S 1n
S 21 S 22 . . .
S = . . . . .
. . . . .
S n1 . . . S nn
ELECTROSTATICS | 51
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.
S-Parameter Calculations
The MEMS interfaces have built-in support for S-parameter calculations. In the
Electric Currents and Electrostatics interfaces, use the terminal boundary feature with
the terminated setting to approximate a connecting transmission line or a voltage
source with a known internal impedance. For a terminal the voltage measurement is
always with respect to ground so at least one ground feature is also required in the
model.
S-Parameter Variables
This module automatically generates variables for the S-parameters. The port names
(use numbers for sweeps to work correctly) determine the variable names. If, for
example, there are two ports with the numbers 1 and 2 and Port 1 is the inport, the
software generates the variables S11 and S21. S11 is the S-parameter for the reflected
wave and S21 is the S-parameter for the transmitted wave. For convenience, two
variables for the S-parameters on a dB scale, S11dB and S21dB, are also defined using
the following relation:
S 11dB = 20 log 10 S 11
The model and physics interface names also appear in front of the variable names so
they can vary. The S-parameter variables are added to the predefined quantities in
appropriate plot lists.
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
You can compute conversions between the impedance matrix, Z, the admittance
matrix, Y, and the S-parameter matrix S in a results table using the Settings window
for the Global Matrix Evaluation node, which you can add under Results>Derived Values.
ELECTROSTATICS | 53
P i e z o r e s i s ti v i t y
The piezoresistive effect describes the change in a material’s conductivity when a stress
is applied to the material. Unlike the piezoelectric effect, piezoresistivity is not
reversible, so an applied current does not induce a stress (unless other secondary effects
are present, such as heating). Piezoresistance is usually associated with semiconductor
materials. In semiconductors, piezoresistance results from the strain-induced alteration
of the material’s band structure and the associated changes in carrier mobility and
number density.
E = J + J
where is the resistivity and is the induced change in the resistivity. In the general
case, and are both rank-2 tensors (matrices). The change in resistance is related
to the stress, , (for the piezoresistance form of the equations) or the strain, , (for the
elastoresistance form of the equations) by the constitutive relationship:
= (piezoresistance form)
(2-13)
= M (elastoresistance form)
where is the piezoresistance tensor (SI unit: Pa1m) and M is the elastoresistance
tensor (SI unit: m). Note that both of these quantities are material properties. and
M are in this case rank-4 tensors; however, they can be represented as matrices if the
resistivity, stress, and strain are converted to vectors within a reduced subscript
notation.
In practice, piezoresistors are fabricated by doping a thin layer of a lightly doped wafer,
either using the process of ion implantation or by diffusing dopant into the wafer from
its surface at high temperatures. Both processes usually produce thin piezoresistors,
typically 100s of nm thick in structures that are microns or tens of microns thick. These
aspect ratios are challenging to resolve with a practical mesh, so the MEMS Module
includes physics interfaces that make it possible to model the current conduction in a
thin layer, represented by a boundary in the model geometry. The Piezoresistivity,
Boundary Currents interface is appropriate for situations in which the thickness of the
piezoresistive layer is much less than that of the structural layer. If the structural layer
is also sufficiently thin that it can be represented using the Shell interface, then the
Piezoresistivity, Shell interface can be used (requires the Structural Mechanics Module
PIEZORESISTIVITY | 55
Thermal Effects in MEMS Devices
In this section:
T- + – k T = Q
C p ------ (2-14)
t
Here is the density of the solid, Cp is its heat capacity, k is Boltzmann’s constant, T
is the temperature, and Q is a term that includes heat sources and sinks. The Heat
Transfer in Solids interface solves Equation 2-14, without additional coupling
assumed.
Equation 2-14 is solved in the Joule Heating interface, along with the equations for
electric current flow. In this case the volumetric heat source Q takes the value:
Q = JE (2-15)
where J is the current density, and E is the electric field. The electric field component
in the direction of the current flow expresses the voltage drop per unit length of
material in the direction of flow. Since the current density gives the current per unit
area through a surface with normal parallel to the direction of flow, the quantity Q is
The Modeling Thermal Stresses section explains how to use this feature
to include residual stresses in the film in addition to the thermal stresses
generated by Joule heating.
Thermal stresses in a single added layer can be modeled in multiphysics interfaces using
a coupled multiphysics interface that include the Solid Mechanics interface by adding
a The Fluid–Solid Interaction Interface subnode under the Multiphysics node (in the
Solid Mechanics interface you can also add it as subnode to the material model). A
reference temperature, Tref, for zero thermal strain must be specified, and the
operating temperature, T, is set by selecting a heat transfer interface in a multiphysics
model, which sets up this temperature coupling automatically, or by selecting a
temperature as a model input (that is, temperatures generated by other physics
interfaces can be selected or a user-defined value or expression for the temperature can
be supplied), if the thermal expansion is defied in the Solid Mechanics interface. If the
device is operated at a constant temperature, then The Solid Mechanics Interface
should be used with T specified as a user-defined constant value in the model input. If
Thermoelastic Damping
When the product of the coefficient of thermal expansion and the temperature of the
structure at rest is large, mechanical vibrations cause significant oscillations of the
temperature field. The resulting energy transfer from the mechanical domain to the
thermal domain is called thermoelastic damping. The Thermoelasticity Interface
provides a predefined multiphysics interface to model thermoelastic damping.
In the MEMS Module, the microfluidics physics interfaces are focused on modeling
the interaction of fluids with microstructures.
In this section:
• Fluid-Structure Interaction
• Thin-Film Damping
Fluid-Structure Interaction
The Fluid-Structure Interaction interface can be used to compute the flow of a liquid
or gas around a moving microstructure. A typical application is the simulation of
microvalves. The physics interface can be used to model both turbulent and laminar
flows at low Mach numbers.
Thin-Film Damping
In many MEMS devices a released structure moves relative to the substrate,
compressing or shearing a thin layer of gas between the two structures. The fluid film
It is often convenient to model the film using the Reynolds equation — enabling the
3D gas layer to be modeled on a 2D surface. The Reynolds equation applies for
isothermal flows of a thin film in which the film height, h, is much less than a typical
lateral dimension of the device, L. It must also be possible to neglect inertial flow
effects in comparison to the viscous effects (for MEMS devices, this assumption is
reasonable below MHz frequencies). Additionally, the curvature of the reference
surface is neglected.
For isothermal, ideal gases the Reynolds equation can be simplified, leading to a
modified Reynolds equation, formulated in the pressure only (the full Reynolds
equation is formulated using both the pressure and the density). COMSOL
Multiphysics includes both forms of the Reynolds equations in the Thin-Film Flow
interfaces, and it is possible to switch between them using the Equation Type setting.
Both forms of the equations include squeeze-film damping and slide-film damping.
The Thin-Film Flow interfaces are available as standalone physics interfaces that can be
manually coupled to, for example, the physics interfaces in the Structural Mechanics
branch. Thin-Film Damping has also been added as a boundary condition to the
structural mechanics physics interfaces — these include Solid Mechanics,
Piezoelectricity, Electromechanics, and Joule Heating and Thermal Expansion. The
boundary condition automatically couples the damping forces to the model.
For gas damping, especially in vacuum packaged devices, the Knudsen number for the
thin film gap is an important parameter to consider. This is the ratio of the gas mean
free path to the gap height. Below Knudsen numbers of 0.01 it is necessary to use a
slip boundary condition for the fluid by specifying a characteristic slip length. As the
Knudsen number decreases further additional rarefaction effects need to be accounted
MICROFLUIDICS | 61
for. These can be included empirically in the Reynolds equation using models that
require the mean free path and the tangential momentum accommodation coefficient
(TMAC) to be specified. The TMAC is a measure of the mean fractional loss of
tangential momentum of the gas molecules on collision with the walls — typical values
for various gas surface combinations are available in Ref. 1 (see Table 6 in this
reference).
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
These features either accept a voltage measurement from the connecting noncircuit
physics interface and return a current from an Electrical Circuit interface or the other
way around.
• A choice is made in the Settings window for the noncircuit physics interface feature,
which then announces (that is, includes) the coupling to the Electrical Circuit
interface. Its voltage or current is then included to make it visible to the connecting
circuit feature.
• A voltage or current that has been announced (that is, included) is selected in a
feature node’s Settings window.
• Apply the voltage or current from the connecting “External” circuit feature as an
excitation in the noncircuit physics interface.
• Define your own voltage or current measurement in the noncircuit physics interface
using variables, coupling operators and so forth.
• In the Settings window for the Electrical Circuit interface feature, selecting the
User-defined option and entering the name of the variable or expression using
coupling operators defined in the previous step.
1 In the Model Builder, right-click the Study node and select Show Default Solver.
Solving
Some modeling errors lead to the error message The DAE is structurally
inconsistent being displayed when solving. This error typically occurs
from having an open current loop, from connecting voltage sources in
parallel, or connecting current sources in series.
In this respect, the predefined coupling features are also treated as (ideal)
voltage or current sources. The remedy is to close current loops and to
connect resistors in series with voltage sources or in parallel with current
sources.
Postprocessing
The Electrical Circuits interface, unlike most of the other physics interfaces, solves for
a relatively large number of global dependent variables (such as voltages and currents),
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. For
resistors it adds a heat source variable as well.
References
1. V. Kaajakari, Practical MEMS, Small Gear Publishing, Las Vegas, 2009.
2. S.D. Senturia, Microsystem Design, Springer Science and Business Media, New
York, 2001.
3. A.F. Bower, Applied Mechanics of Solids, CRC Press, Boca Raton, FL, 2010
(http://www.solidmechanics.org).
SPICE Import
The circuit definition in COMSOL Multiphysics adheres to the SPICE format
developed at the University of California, Berkeley (Ref. 1). SPICE netlists can be
imported and the corresponding circuit nodes are generated in the COMSOL
Multiphysics model. Most circuit simulators can export to this format or some version
of 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
SPICE Export
The SPICE Export functionality creates a SPICE netlist file containing a description of
the circuit represented by the physics interface. This functionality can be accessed from
the physics interface context menu (right-click the physics interface node and select
Export SPICE Netlist). After specifying a filename, the circuit is exported and messages
from the export process display in the Messages window. During the export process, a
series of operations are performed:
The title of the exported netlist file is the model’s filename, and the time, date, and
version of COMSOL Multiphysics is added as a comment in the netlist file.
AC/DC Interfaces
This chapter describes the physics interfaces found under the AC/DC branch .
In this chapter:
See the COMSOL Multiphysics Reference Manual for details about The Magnetic
Fields Interface and the Theory of Magnetic Fields.
71
The Electrostatics Interface
The Electrostatics (es) interface ( ), found under the AC/DC>Electric Fields and
Currents 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 when it is used 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.
Physics-Controlled Mesh
The physics-controlled mesh is controlled from the Mesh node’s Settings window (if the
Sequence type is Physics-controlled mesh). There, in the table in the Physics-Controlled
Mesh section, find the physics interface in the Contributor column and select or clear
the check box in the Use column on the same table row for enabling (the default) or
disabling contributions from the physics interface to the physics-controlled mesh.
Information from the physics, such as the presence of an infinite elements domain or
periodic condition, will be used to automatically set up an appropriate meshing
sequence.
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.
Select the Use 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.
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
Select the shape order for the Electric potential dependent variable — Linear, Quadratic
(the default), Cubic, Quartic, or Quintic. For more information about the Discretization
section, see Settings for the Discretization Sections in the COMSOL Multiphysics
Reference Manual.
In the COMSOL Multiphysics Reference Manual, see Table 2-4 for links
to common sections and Table 2-5 to common feature nodes. You can
also search for information: press F1 to open the Help window or Ctrl+F1
to open the Documentation window.
n2 D1 – D2 = s
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 3-1 for a list of interior and exterior
boundary conditions, including edge, point, and pair availability.
• Force Calculation
• Ground
1
This feature is available with the Piezoelectricity multiphysics interface.
2
Requires either the Acoustics Module, MEMS Module, or Structural Mechanics
Module.
3
This feature is available with the Ferroelectroelasticity multiphysics interface.
Table 3-1 lists the interior and exterior boundary conditions available with this physics
interface. It also includes edge, point, and pair availability.
TABLE 3-1: INTERIOR AND EXTERIOR BOUNDARY CONDITIONS (INCLUDING EDGE, POINT, AND PAIR
AVAILABILITY) FOR THE ELECTROSTATICS INTERFACE.
In the COMSOL Multiphysics Reference Manual, see Table 2-4 for links
to common sections and Table 2-5 to common feature nodes. You can
also search for information: press F1 to open the Help window or Ctrl+F1
to open the Documentation window.
Charge Conservation
The Charge Conservation node adds the equations for charge conservation according to
Gauss’ law for the electric displacement field. It provides an interface for defining the
constitutive relation and its associated properties such as the relative permittivity.
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
• Relative permittivity (the default) to use the constitutive relation D0rE. Then
the default is to take the Relative permittivity r (dimensionless) values From material.
For User defined, select Isotropic, Diagonal, Symmetric, or Full and enter values or
expressions in the field or matrix. The default is 1.
• Polarization to use the constitutive relation D0E 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 D0rE 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 Full 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.
• Dispersion to use the constitutive relation D 0E PE rS, where the
polarization vector is calculated from the electric field using the dielectric dispersion
model. This option is available when the Material type is set to Solid. You enter the
Relative permittivity rS (dimensionless) values From material. For User defined,
select Isotropic, Diagonal, or Symmetric and enter values or expressions in the field or
matrix. This value of electric permittivity will be used in stationary study, for which
the polarization vector is calculated as P 0rS E.
CONDUCTION CURRENT
By default, the Electrical conductivity for the media is defined From material. You can
also select User defined or Linearized resistivity.
• For User defined select Isotropic, Diagonal, Symmetric, or Full depending on the
characteristics of the electrical conductivity, and then enter values or expressions for
the Electrical conductivity in the field or matrix.
• For Linearized resistivity the default Reference temperature Tref, and Resistivity
temperature coefficient , and Reference resistivity 0 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
appears in the Model Inputs section.
Context Menus
Electrostatics>Charge Conservation>Conduction Loss (Time-Harmonic)
Ribbon
Physics tab with Charge Conservation node selected in the model tree:
Initial Values
The Initial Values node adds an initial value for the electric potential V that can serve
as an initial condition for a transient simulation or as an initial guess for a nonlinear
solver.
Force Calculation
Use the Force Calculation node to define globally available force and torque variables
for the selected domains.
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. The computed force
is made available as a global vector variable with name
<name>.Force<component>_<force name>, where <name> is the Name of the physics
interface.
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, <name>.T<component>_<force name>. The
resulting torque component parallel to the given Torque axis is given as a global scalar
variable, typically <name>.Tax_<force name>.
In frequency domain studies, the cycle-averaged force and torque variables are also
defined, with names <name>.Forceav<component>_<force name>,
ADVANCED SETTINGS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
The Compute both instantaneous and cycle-averaged variables in the frequency domain
check box is selected by default. This 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 More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
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 More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
–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.
SYMMETRY TYPE
Choose between Symmetry and Antisymmetry for the electric field. If Symmetry is
chosen, the boundary condition is
nE = 0
which states that the normal component of the electric field is zero. If Antisymmetry is
chosen, the boundary condition becomes
nE = 0
which states that the tangential component of the electric field is zero.
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 section, which you activate by right-clicking
the Periodic Condition node and select Manual Destination Selection. With that section,
the boundaries which constitute the source and destination surfaces can be manually
specified.
BOUNDARY SELECTION
The software usually automatically identifies the boundaries as either source boundaries or
destination boundaries, as indicated in the selection list. This works fine for cases like opposing
parallel boundaries. In other cases, right-click Periodic Condition and select Manual Destination
Selection to control the destination. By default it contains the selection that COMSOL
Multiphysics identifies.
DESTINATION SELECTION
This section is available for specifying the destination boundaries, if needed, when the
Manual Destination Selection option is selected in the context menu for the Periodic
Condition node. You can only select destination boundaries from the union of all source
and destination boundaries.
• Continuity to make the potential periodic (equal on the source and destination).
• Antiperiodicity to make it antiperiodic.
• 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 More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
ORIENTATION OF SOURCE
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. For information about the
Orientation of Source section, see Orientation of Source and Destination.
ORIENTATION OF DESTINATION
This section appears if the setting for Transform to intermediate map in the Orientation
of Source section is changed from the default value, Automatic, and Advanced Physics
Options is selected in the Show More Options dialog box. For information about the
Orientation of Destination section, see Orientation of Source and Destination.
0 r
n D 1 = ---------- V 1 – V 2
d
0 r
n D 2 = ---------- V 2 – V 1
d
to model a thin gap of a material with a small permittivity compared to the adjacent
domains. The layer has the thickness d and the relative permittivity r. The indices 1
and 2 refer to the two sides of the boundary.
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 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.
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
Isotropic, Diagonal, Symmetric, or Full and enter values or expressions in the field or
matrix. Enter a Surface thickness ds (SI unit: m) of the shielding. The default is 1 m.
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.
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, Terminated (power), or
Terminated (voltage). 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.
• Voltage to specify a fixed electric potential at the terminal. Enter an electric potential
V0. The default is 1 V.
• Circuit to specify a terminal connected to an Electrical Circuit interface that models
an external circuit. The Terminal node provides a current-voltage characteristic to
the circuit.
• Terminated (power) to connect the terminal to an impedance that might represent a
load or a transmission line. When Terminated (power) is selected, the scattering
parameters (S-parameters) are computed. The Terminated (power) setting can only
be used in frequency-domain studies.
- Enter a Terminal power P0 to specify the input power at this terminal from the
transmission line. This excitation can be toggled using a port sweep.
- Select an option from the Characteristic impedance from list to define the value of
the impedance — Physics interface or User defined. For Physics interface, the
Reference impedance Zref defined in the Settings window for the physics under
Sweep Settings is used. For User defined, enter a specific impedance Zref for this
terminal. The default is 50 .
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.
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 More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
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:
–t V = Et = 0
n E = 0.
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 in 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.
The Initial value for charge is available for the Circuit excitation. It specifies the initial
charge on the floating potential when the circuit is connected. Enter an Initial value for
charge Qinit.
Distributed Capacitance
The Distributed Capacitance node adds a distributed capacitance boundary condition
according to the following equations for exterior boundaries:
V ref – V
– n D = 0 r --------------------
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 103 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 in the Physics toolbar. For more
information see Harmonic Perturbation — Exclusive and Contributing Nodes in the
COMSOL Multiphysics Reference Manual.
EDGE SELECTION
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 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 CHARGE
Enter a value or expression to apply a Point charge QP (SI unit: C) to points. This
source represents an electric displacement field flowing out of the point. The default is
0 C.
POINT SELECTION
Use the Change Cross Section node to set the cross-section area for specific geometric
entities.
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 Full and enter values or expressions in the field or
matrix.
Dispersion
The Dispersion subnode to Charge Conservation and Charge Conservation, Piezoelectric
allows you to model possible dielectric losses in the material. The effect can be
accounted for in Eigenfrequency, Frequency Domain and Time Dependent study.
Under the Material model, you can select either Debye (default) or Multipole Debye
dispersion option.
For a single pole Debye dispersion choice, you can enter the Relaxation time and Relative
permittivity contribution. In the Thermal Effects tab, you can select the Shift function
type, which will be used to compute the effective relaxation time accounting for the
temperature effects.
For the Multipole Debye dispersion, you can select two options to enter the Relaxation
Data.
• Relaxation time and Relative permittivity change choice will provide an input table
with dynamically changed size, where you can enter any number of Relaxation time
and the corresponding Relative permittivity contribution values. In the Thermal
Debye Dispersion
For all cases, you can specify how the relative permittivity input, rS, on the parent
node should be interpreted by selecting the Static response (the default value is High
frequency limit). This setting will have effect in Eigenfrequency, Frequency Domain,
and Time Dependent study types.
THERMAL EFFECTS
Dispersion properties depend on the temperature. For many materials, a change in the
temperature can be transformed directly into a change in the time scale. Thus, the
relaxation time is modified to aTTm, where aT(T) is a shift function.
• When the default, None, is kept, the shift function aTT is set to unity and the
relaxation time is not modified.
DISCRETIZATION
To display this section, click the Show More Options button ( ) and select
Discretization from the Show More Options dialog box. Select the element order from
the list box for the auxiliary electric field vector variables en.
You can change the solution algorithm by using the check box Use local time
integration (checked by default). When unchecked, the solution of the equations for
the auxiliary dependent variables in Time Dependent study will be governed by the
time stepper algorithm used in the solver.
Context Menus
Electrostatics>Charge Conservation>Dispersion
Electrostatics>Charge Conservation, Piezoelectric>Dispersion
Attributes>Dispersion
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.
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.
Information from the physics will be used to automatically set up an appropriate mesh
sequence.
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 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 in 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 3D components, from the Condition for the x = x0 plane, Condition for
the y = y0 plane, and Condition for the z = z0 plane lists, choose Off (the
default), Symmetric, or Antisymmetric. Then enter the value for the plane
location x0, y0, or z0 (SI unit: m) as required.
FAR-FIELD APPROXIMATION
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
For more information about the Far Field Approximation settings, see
Far-Field Approximation Settings in the COMSOL Multiphysics
Reference Manual.
INFINITY CONDITION
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.
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).
Charge Conservation
The Charge Conservation node adds the equations for charge conservation according to
Gauss’ law for the electric displacement field. It provides an interface for defining the
constitutive relation defined by the relative permittivity.
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.
Floating Potential
The Floating Potential node is used when modeling a metallic electrode at floating
potential.
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.
–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.
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 as
Secondary interface.
The physics interface solves a current conservation equation based on Ohm’s law using
the scalar electric potential as the dependent variable. The interface is typically used to
model good conductors for which Ohm’s law applies. It is not suitable for modeling
poor conductors like air and insulators.
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.
Physics-Controlled Mesh
The physics-controlled mesh is controlled from the Mesh node’s Settings window (if the
Sequence type is Physics-controlled mesh). There, in the table in the Physics-Controlled
Mesh section, find the physics interface in the Contributor column and select or clear
the check box in the Use column on the same table row for enabling (the default) or
disabling contributions from the physics interface to the physics-controlled mesh.
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.
THICKNESS (2D)
Enter a default value for the Out-of-plane thickness d (SI unit: m) (see Equation 3-9).
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 Use 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
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.
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 share
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
Select the shape order for the Electric potential dependent variable — Linear, Quadratic
(the default), Cubic, Quartic, or Quintic. For more information about the Discretization
section, see Settings for the Discretization Sections in the COMSOL Multiphysics
Reference Manual.
• Domain, Boundary, Edge, Point, and Pair Nodes for the Electric
Currents Interface
• Theory for the Electric Currents Interface
n2 J1 – J2 = 0
1
This feature is available with the Piezoresistivity, Domain Currents multiphysics
interface, which requires the MEMS Module.
In the COMSOL Multiphysics Reference Manual, see Table 2-4 for links
to common sections and Table 2-5 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 3-1 lists the interior and exterior boundary conditions available with this physics
interface. It also includes edge, point, and pair availability.
TABLE 3-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 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
User Defined
For User defined select Isotropic, Diagonal, Symmetric, or Full 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:
• Dielectric losses: uses the constitutive relation D 0' j"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 Full and enter values or
expressions in the field or matrix. The default is 1. Note that the material parameters
Relative permittivity (real part) ' and the Relative permittivity (imaginary part) "
form the complex relative permittivity r ' – j''. The time-harmonic Sign
Convention requires a lossy material to have a positive material parameter ''.
• Loss tangent, loss angle: uses the constitutive relation D 0' jtanE. Specify
the Relative permittivity (real part) ' (dimensionless) and Loss angle (SI unit: rad).
• Loss tangent, dissipation factor: uses the constitutive relation D 0' jtanE.
Specify the Relative permittivity (real part) ' (dimensionless) and the Dissipation
factor tan (dimensionless).
• Dispersion to use the constitutive relation D 0E PE rS, where the
polarization vector is calculated from the electric field using the dielectric dispersion
model. This option is available only when the Material type is set to Solid. You enter
the Relative permittivity rS (dimensionless) values From material. For User defined,
select Isotropic, Diagonal, or Symmetric and enter values or expressions in the field or
matrix. This value of electric permittivity will be used in stationary study, for which
the polarization vector is calculated as P 0rS E.
Under the Material model, you can select either Debye (default) or Multipole Debye
dispersion option.
For a single pole Debye dispersion choice, you can enter the Relaxation time and Relative
permittivity contribution. In the Thermal Effects tab, you can select the Shift function
type, which will be used to compute the effective relaxation time accounting for the
temperature effects.
For the Multipole Debye dispersion, you can select two options to enter the Relaxation
Data.
• Relaxation time and Relative permittivity change choice will provide an input table
with dynamically changed size, where you can enter any number of Relaxation time
and the corresponding Relative permittivity contribution values. In the Thermal
Effects tab, you can select the Shift function type, which will be used to compute the
effective relaxation times accounting for the temperature effects.
• Constant Loss Tangent option. In this case, you enter the Loss Tangent (fc) together
with the Center frequency fc. You also specify the model Bandwidth (decades) that
defines a frequency interval centered at fc, in which the loss tangent will be
approximately constant and equal to (fc). You can also select the Accuracy for the
approximation to be either Normal (default) or High. The software will automatically
deduce the necessary number of Debye poles together with the values of the
corresponding relaxation times and relative permittivity contributions, which will be
used in computations to maintain the requested bandwidth and accuracy.
Debye Dispersion
For all cases, you can specify how the relative permittivity input, rS, on the parent
node should be interpreted by selecting the Static response (the default value is High
frequency limit). This setting will have effect in Eigenfrequency, Frequency Domain,
and Time Dependent study types.
THERMAL EFFECTS
Dispersion properties depend on the temperature. For many materials, a change in the
temperature can be transformed directly into a change in the time scale. Thus, the
relaxation time is modified to aTTm, where aT(T) is a shift function.
• When the default, None, is kept, the shift function aTT is set to unity and the
relaxation time is not modified.
• For Vogel-Fulcher enter values or expressions for these properties:
- Reference temperature T0. The default is 230 K.
- Activation energy Q. The default is 8000 J/mol.
DISCRETIZATION
To display this section, click the Show More Options button ( ) and select
Discretization from the Show More Options dialog box. Select the element order from
the list box for the auxiliary electric field vector variables en.
You can change the solution algorithm by using the check box Use local time
integration (checked by default). When unchecked, the solution of the equations for
the auxiliary dependent variables in Time Dependent study will be governed by the
time stepper algorithm used in the solver.
Context Menus
Electric Currents>Current Conservation>Dispersion
Ribbon
Physics tab with Current Conservation node selected in the model tree:
Attributes>Dispersion
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.
J = E + J e
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.
Electric Insulation
The Electric Insulation node, which is the default boundary condition, adds electric
insulation as the boundary condition:
nJ = 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 = I 0 ,
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
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 current, or connected to a circuit,
it behaves like a terminal.
The Harmonic Perturbation subnode (it is of the exclusive type) is available from the
context menu (right-click the parent node) or in 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 in 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 = Qj
–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 in 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 in 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 = ----
1
r + jC s V – V ref
s
n J 1 – J 2 = ------- V – V ref + 0 r V – V ref
1
dL t
n J 1 – J 2 = ----- V – V ref + C s V – V ref
1
s 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 in 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.
• 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.
• Voltage to specify a fixed electric potential at the terminal. Enter the value of the
electric potential V0 (SI unit: V). The default is 1 V.
• Circuit to specify a terminal connected to an Electrical Circuit interface that models
an external circuit.
• Power to specify a fixed electrical power for the terminal. The power is computed as
the product of the terminal current and voltage and, in the frequency domain, is the
cycle-averaged power. When Power is selected:
- Enter a Terminal power P0 (SI unit: W) to specify the input power at this terminal.
The default is 1 W.
• Terminated (power) to connect the terminal to an impedance that might represent a
load or a transmission line. When Terminated (power) is selected, the scattering
parameters (S-parameters) are computed.
- Enter a Terminal power P0 (SI unit: W) to specify the input power at this terminal
from the transmission line. This excitation can be toggled using a terminal sweep.
The default is 1 W.
- Select an option from the Characteristic impedance from list to define the value of
the impedance — Physics interface (the default) or User defined. For Physics
interface the Reference impedance Zref defined in the Settings window for the
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
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.
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 More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
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.
THICKNESS
Enter a value or expression for the Surface thickness ds (SI unit: m). The default is
1 cm.
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:
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.
+ j 0 r
n J 1 = -------------------------------- V 1 – V 2
ds
+ j 0 r
n J 2 = -------------------------------- V 2 – V 1
ds
n J 1 = ----
1- + jC V – V
s 1 2
s
n J 2 = ----
1
- + jC s V 2 – V 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
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.
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 More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box. The constraint setting can be set
to Pointwise constraints (default), or Weak constraints. Additionally, the Nitsche
constraints option is available for the scalar potential in the Electric Currents, Magnetic
Field, No Currents, and Rotating Machinery, Magnetic physics interfaces. The
displacement current term is excluded from the Nitsche constraints in the Electric
Currents interface and is therefore valid in the quasistationary approximation.
POINT SELECTION
POINT SELECTION
DIPOLE SPECIFICATION
Select a Dipole specification — Magnitude and direction (the default) or Dipole moment.
• 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).
• Specify an 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: m/
Pa).
For a Elastoresistance form, select an Elastoresistance coupling matrix Ml
(SI unit: m).
When this physics interface is added, it adds a default Ground Node feature and
associates that with node zero in the electrical circuit.
Circuit nodes are nodes in the electrical circuit (electrical nodes) and
should not be confused with nodes in the Model Builder tree of the
COMSOL Multiphysics software. Circuit node names are not restricted
to numerical values but can contain alphanumeric characters.
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is cir.
GROUND CONNECTION
Set the Node name for the ground node in the circuit. The convention is to use 0 (zero)
for the ground node. If adding more ground nodes, each must have a unique node
name (number).
Voltmeter
The Voltmeter ( ) feature connects a voltmeter (voltage measurement device)
between two nodes in the electrical circuit. A voltmeter behaves electrically as an open
circuit. The voltmeter node adds a Probe sampling the voltage across it.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
Ampère Meter
The Ammeter ( ) feature connects an ammeter (current measurement device)
between two nodes in the electrical circuit. An ammeter behaves electrically as a short
circuit. The ammeter node adds a Probe sampling the current through it.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
DEVICE PARAMETERS
Enter the Resistance of the resistor.
Resistor
The Resistor ( ) feature connects a resistor between two nodes in the electrical
circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the resistor.
Capacitor
The Capacitor ( ) feature connects a capacitor between two nodes in the electrical
circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the capacitor.
DEVICE PARAMETERS
Enter the Capacitance of the capacitor.
Inductor
The Inductor ( ) feature connects an inductor between two nodes in the electrical
circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the inductor.
DEVICE PARAMETERS
Enter the Inductance of the inductor.
Voltage Source
The Voltage Source ( ) feature connects a voltage source between two nodes in the
electrical circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the voltage source. The first node
represents the positive reference terminal.
• For a General source, the Voltage Vsrc (default value: 1 V). General sources are active
in Stationary, Time-Dependent and Frequency Domain studies.
• For an AC-source: the Voltage Vsrc (default value: 1 V) and the Phase (default
value: 0 rad). AC-sources are active in Frequency Domain studies only.
• For a sine source: the Voltage Vsrc (default value: 1 V), the Offset Voff (default value:
0 V), the Frequency (default value: 1 kHz), and the Phase (default value: 0 rad).
The sine sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
• For a pulse source: the Voltage Vsrc (default value: 1 V), the Offset Voff (default value:
0 V), the Delay td (default value: 0s), the Rise time tr and Fall time tf (default values:
0 s), the Pulse width pw (default value: 1 s), and the Period Tper (default value: 2 s).
The pulse sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
For the AC source, the frequency is a global input set by the solver. AC
sources should be used in Frequency-domain studies only. Do not use the
Sine source unless the model is time dependent.
Current Source
The Current Source ( ) feature connects a current source between two nodes in the
electrical circuit.
NODE CONNECTIONS
Set the two Node names for the connecting nodes for the current source. The first node
represents the positive reference terminal from where the current flows through the
source to the second node.
• For a General source, the Current isrc (default value: 1 A). General sources are active
in Stationary, Time-Dependent and Frequency Domain studies.
• For an AC-source: the Current isrc (default value: 1 A) and the Phase (default
value: 0 rad). AC-sources are active in Frequency Domain studies only.
• For a sine source: the Current isrc (default value: 1 A), the Offset ioff (default value:
0 A), the Frequency (default value: 1 kHz), and the Phase (default value: 0 rad).
The sine sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
• For a pulse source: the Current isrc (default value: 1 A), the Offset ioff (default value:
0 A), the Delay td (default value: 0 s), the Rise time tr and Fall time tf (default values:
0 s), the Pulse width pw (default value: 1 s), and the Period Tper (default value: 2 s).
The pulse sources are active in Time-Dependent studies and also in Stationary
studies, providing that a value for t has been provided as a model parameter or
global variable.
For the AC source, the frequency is a global input set by the solver. AC
sources should be used in frequency-domain studies only. Do not use the
Sine source unless the model is time dependent.
NODE CONNECTIONS
Specify four Node names: the first pair for the connection nodes for the voltage source
and the second pair defining the input control voltage. The first node in a pair
represents the positive reference terminal.
NODE CONNECTIONS
Specify four Node names: the first pair for the connection nodes for the current source
and the second pair defining the input control voltage. The first node in a pair
represents the positive voltage reference terminal or the one from where the current
flows through the source to the second node.
DEVICE PARAMETERS
There are two options to define the relationship between the control voltage and
resulting current. The Use gain method defines the resulting current to be the control
voltage multiplied by the gain (SI units: S). The Custom expression method can define
the relationship with an arbitrary expression.
NODE CONNECTIONS
Set two Node names for the connection nodes for the voltage source. The first node in
a pair represents the positive reference terminal.
DEVICE PARAMETERS
There are two options to define the relationship between the control current and
resulting voltage. The Use gain method defines the resulting voltage to be the control
current multiplied by the gain (SI units: ). The Custom expression method can define
the relationship with an arbitrary expression.
NODE CONNECTIONS
Specify two Node names for the connection nodes for the current source. The first node
in a pair represents the positive reference terminal from where the current flows
through the source to the second node.
DEVICE PARAMETERS
There are two options to define the relationship between the control current and
resulting current. The Use gain method defines the resulting current to be the control
current multiplied by the gain. The Custom expression method can define the
relationship with an arbitrary expression.
Switch
The Switch ( ) feature is used to connect or disconnect the conducting path in a
circuit under specific conditions.
NODE CONNECTIONS
Specify two Node names for the connection nodes for the current source. The first node
in a pair represents the positive reference terminal from where the current flows
through the source to the second node.
SWITCH CONDITIONS
There are three types of conditions, Voltage controlled, Current controlled, and Custom
expressions. For each type of condition there are two conditions, one for turn on and
one for turn off. The on condition is true if the On condition expression is larger than
zero, while the off condition is true if the Off condition is less than zero.
The Initial state list has three options, Use on condition, Use off condition, and Boolean
expression. The two former options mean that the switch will have an initial state
matching to the on or off condition. The third option makes the switch's initial state
match a custom Boolean expression. Separating on, off, and initial states makes the
switch more flexible and can support Schmitt-trigger style switches and various latches.
SWITCH PARAMETERS
When the switch is in the on state it has a nonzero resistance specified by the On
resistance expression. For the off state no current flows through the switch (infinite
resistance). There is also a transition time for the switch to turn on and off set by the
Switching time expression. The switch triggers an implicit event that updates a discrete
state variable (with suffix _state).
Subcircuit Definition
The Subcircuit Definition ( ) feature is used to define subcircuits, which can be
inserted as devices into the main circuit using Subcircuit Instance nodes. Create the
subcircuit by adding subnodes to the Subcircuit Definition node, either by using the
Physics toolbar, or by right-clicking the Subcircuit Definition.
SUBCIRCUIT PINS
Define the Pin names at which the subcircuit connects to the main circuit or to other
subcircuits when referenced by a Subcircuit Instance node. The Pin names refer to
circuit nodes in the subcircuit. The order in which the Pin names are defined is the
order in which they are referenced by a Subcircuit Instance node. The devices
constituting the subcircuit should be connected only to the subcircuit’s pins and to
themselves.
INPUT PARAMETERS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options. Specify input parameters to a subcircuit that can be changed from a
subcircuit instance. These input parameters can be used in all expression-style edit
fields that affect the parameters of a device, for example, resistance, capacitance, and
current gain. In this way, a subcircuit can represent a parameterized custom device
model.
NODE CONNECTIONS
Select the Name of subcircuit link from the list of defined subcircuits in the circuit model
and the circuit Node names at which the subcircuit instance connects to the main circuit
or to another subcircuit if used therein.
INPUT PARAMETERS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options. Specify input parameters to a subcircuit that can be changed from a
subcircuit instance. These input parameters can be used in all expression-style edit
fields that affect the parameters of a device, for example, resistance, capacitance, and
current gain. In this way, a subcircuit can represent a parameterized custom device
model.
NODE CONNECTIONS
Specify three Node names for the connection nodes for the BJT device. These represent
the collector, base, and emitter nodes for the NPN transistor, and the emitter, base, and
collector nodes for the PNP transistor.
NODE CONNECTIONS
Specify four Node names for the connection nodes for the n-Channel MOSFET or
p-Channel MOSFET device. These represent the drain, gate, source, and bulk nodes,
respectively.
MODEL PARAMETERS
Specify the Model Parameters. Reasonable defaults are provided but for any particular
MOSFET, the device manufacturer should be the primary source of information.
Mutual Inductance
The Mutual Inductance allows specifying a coupling between two existing Inductor
features in the circuit. The mutual inductance of the coupling is
where k is the coupling factor and L1 and L2 are the inductances of the inductors.
DEVICE PARAMETERS
Enter values or expressions for the:
• Coupling factor k (dimensionless). The value must be between 0 and 1, and the
default is 0.98.
• First inductance L1 (SI unit: H) and Second inductance L2 (SI unit: H). These must
be set to two different Inductor features in the circuit.
Transformer
The Transformer feature represents either a combination of two Inductor and a Mutual
Inductance features, or an ideal transformer.
NODE CONNECTIONS
Enter or edit the table in the Node names column for the primary and secondary node
connections.
DEVICE PARAMETERS
Choose a Transformer model — Specify inductors (the default) or Ideal transformer.
For Ideal transformer enter values or expressions for the Winding ratio N1/N2
(dimensionless). The default is 10.
Diode
The Diode device model ( ) is a large-signal model for a diode. It is an advanced
device model and no thorough description and motivation of the many input
parameters are attempted here. The interested reader is referred to Ref. 1 for more
details on semiconductor modeling within circuits. Many device manufacturers
NODE CONNECTIONS
Specify two Node names for the positive and negative nodes for the Diode device.
MODEL PARAMETERS
Specify the Model Parameters. Reasonable defaults are provided but for any particular
diode, the device manufacturer should be the primary source of information.
External I vs. U
The External I vs. U ( ) feature connects an arbitrary voltage measurement (for
example, a circuit terminal or circuit port boundary or a coil domain from another
physics interface) as a voltage source between two nodes in the electrical circuit. The
resulting circuit current from the first node to the second node is typically coupled
back as a prescribed current source in the context of the voltage measurement.
NODE CONNECTIONS
Specify the two Node names for the connecting nodes for the voltage source. The first
node represents the positive reference terminal.
EXTERNAL DEVICE
Enter the source of the Voltage. If circuit or current excited terminals or circuit ports
are defined on boundaries or domains or a multiturn coil domains is defined in other
physics interfaces, these display as options in the Voltage list. Also select the User defined
option and enter your own voltage variable, for example, using a suitable coupling
operator. For inductive or electromagnetic wave propagation models, the voltage
measurement must be performed as an integral of the electric field because the electric
External U vs. I
The External U vs. I ( ) feature connects an arbitrary current measurement (for
example, a coil domain from another physics interface) as a current source between
two nodes in the electrical circuit. The resulting circuit voltage between the first node
and the second node is typically coupled back as a prescribed voltage source in the
context of the current measurement.
NODE CONNECTIONS
Specify the two Node names for the connecting nodes for the current source. The
current flows from the first node to the second node.
EXTERNAL DEVICE
Enter the source of the Current. Voltage excited terminals or lumped ports defined on
boundaries in other physics interfaces are natural candidates but do not appear as
options in the Voltage list because those do not have an accurate built-in current
External I-Terminal
The External I-Terminal ( ) feature connects an arbitrary voltage-to-ground
measurement (for example, a circuit terminal from another physics interface) as a
voltage-to-ground assignment to a node in the electrical circuit. The resulting circuit
current from the node is typically coupled back as a prescribed current source in the
context of the voltage measurement. This node does not apply when coupling to
inductive or electromagnetic wave propagation models because then voltage must be
defined as a line integral between two points rather than a single point measurement
of electric potential. For such couplings, use the External I vs. U node instead.
NODE CONNECTIONS
Set the Node name for the connecting node for the voltage assignment.
EXTERNAL TERMINAL
Enter the source of the Voltage. If circuit- or current-excited terminals are defined on
boundaries in other physics interfaces, these display as options in the Voltage list. Also
See SPICE Import and Export about the supported SPICE commands.
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.
• Debye Dispersion
The different physics interfaces involving only the scalar electric potential can be
interpreted in terms of the charge relaxation process. The fundamental equations
involved are Ohm’s law for the conduction current density
J c = E
-----
- + J = 0
t c
E =
By combining these, one can deduce the following differential equation for the space
charge density in a homogeneous medium
–t
t = 0 e
where
= ---
is called the charge relaxation time. For a good conductor like copper, is of the order
of 1019 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 3-3 below,
TABLE 3-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 D0E
P between D and E, you can rewrite Gauss’ law as a variant of Poisson’s equation
– 0 V – P =
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 with very
little influence from the capacitive layers. Then the Magnetic Fields interface is the
appropriate modeling tool.
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.
– 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.
H = J (3-1)
D J
J = Ji + + p (3-2)
t
where
J i = E
D = 0 E
is the electric displacement in the high frequency limit (in the sense that the excitation
duration is smaller than the shortest charge relaxation time).
Jp = Jm
m=1
J m E
m + J m = 0 m (3-3)
t t
where m is the relaxation time, and m is the relative permittivity contribution. Note
that both quantities can be taken as diagonal matrices to cover the anisotropic case.
e m
J m = 0 m
t
Equation 3-3 can be integrated once with respect to time assuming zero initial fields
e m
m + em = E (3-4)
t
Jp = 0 m e m (3-5)
t
m=1
or
–1
Jp = 0 m m E – e m
m=1
By applying the divergence operator, Equation 3-1 yields the current conservation
equation:
J = 0
Assume that the electric conductivity is small, so that the induction current can be
neglected. Using Equation 3-5, the current conservation equation can be written as
N
E
0 + 0 m e m = 0
t t
m=1
which can be integrated once with respect to time assuming zero initial fields to give
the charge conservation equation
N
0 E + 0 m e m = 0 (3-6)
m=1
Equation 3-6 can be solved using the scalar electric potential V as a dependent variable
determine the electric field as E V, and it needs to be solved together with N
vector Equation 3-4 for N vector dependent variables em. In time domain, these
auxiliary field variables can be treated as the local state variables, and each of N
corresponding ODE equations can be integrated locally.
FREQUENCY DOMAIN
In frequency domain, Equation 3-5 for the auxiliary variables can be solved
analytically, so that Equation 3-6 gives
where is the angular frequency. By separating the real and imaginary parts, the
equation can be further rewritten as
0 ' – j'' E = 0
where
N
2 –1
' = + m I + m
m=1
and
N
2 –1
'' = m m I + m
m=1
Each term in the sum in the above equations represents a pole on the complex w plane.
Thus, the model is referred to as Multipole Debye Dispersion model. To determine
the number of poles needed to approximate the lossy material behavior can be a
difficult task. One possible approach is presented in Constant Loss Tangent section.
EIGENFREQUENCY STUDY
Direct use of the complex frequency-dependent permittivity would lead to a nonlinear
eigenvalue problem, for which the solution would depend on the eigenvalue
linearization point (transform point frequency, specified on the solver). The solution
would be accurate only for the eigenfrequency that is close to the transform point
frequency. Thus, the eigenfrequencies would need to be computed one by one using
iterative updates of the transform point frequency value.
Instead, one can solve N algebraic equations for the auxiliary variables em
j m e m + e m = E (3-7)
together with Equation 3-7. Even though this approach requires extra degrees of
freedom, it will produce frequency independent contributions to the damping and
stiffness matrices. The corresponding linear damped eigenvalue problem can be solved
In Ref. 2, a special parameterization was suggested for the multipole Debye model
based on the original result derived in Ref. 3 for equivalent electric circuits. It allows
to fit the model parameters for Multipole Debye Model so that the loss tangent is
nearly constant at a certain frequency range.
The input parameters are the loss tangent (fc) together with the center frequency fc,
at which is has been measured, and the model bandwidth nd that defines a frequency
interval (in decades) centered at fc, in which the loss tangent will be approximately
constant and equal to (fc).
The relaxation times are computed as m 1m, where the corresponding frequencies
are equidistantly spaced in the logarithmic space as
N + 1 – 2m 2
m = c 10
m–1
m = k – 1 k
where k > 1 is the spacing parameter. Note that the relative electric permittivity in the
N
low and high frequency limits are related simply as r f 0 = k .
f = f c 6%
for
–nd 2 nd 2
10 f f c 10
2
which requires N = n d + 3 poles, and the spacing parameter k = 10 with the
loss angle d computed from (fc).
f = f c 2%
in the same interval, which will require N = 3n d 2 + 5 poles, and the spacing
4 3
parameter k = 10 .
TEMPERATURE EFFECTS
For many dielectric materials such as polymers, the charge relaxation properties have a
strong dependence on the temperature. A common assumption is that a change in the
temperature can be transformed directly into a change in the time scale. The reduced
time is defined as
t
dt'
tr = ------------------------
T T t'
0
The implication is that the problem can be solved using the original material data,
provided that the time is transformed into the reduced time. Think of the shift
function T(T) as a multiplier to the relaxation time in Equation 3-4 and
Equation 3-8.
Vogel-Fulcher Shift
The shift function is given by
Q
log T = -------------------------
-
R T – T0
here, a base-e logarithm is used, Q is the activation energy (SI unit: J/mol), and R is
the universal gas constant. Temperature T0 is the temperature at which the relaxation
time would becomes infinite.
Arrhenius Shift
The shift function is assumed to follow an Arrhenius law
here, a base-e logarithm is used, Q is the activation energy (SI unit: J/mol), and R is
the universal gas constant.
Williams–Landel–Ferry Shift
One commonly used shift function is defined by the WLF (Williams–Landel–Ferry)
equation:
– C 1 T – T ref
log T = -------------------------------------
-
C 2 + T – T ref
where a base-10 logarithm is used. This shift is only valid over a certain temperature
range, typically around the glass transition temperature.
The first step to compute the shift function T(T) consists of building a master curve
based on experimental data. To do this, the curves of the dielectric properties versus
time or frequency are measured at a reference temperature Tref. Then, the same
properties are measured at different temperatures.
Since the master curve is measured at an arbitrary reference temperature Tref, the shift
factor T(T) can be derived with respect to any temperature, and it is commonly taken
as the shift with respect to the glass transition temperature.
Tool–Narayanaswamy–Moynihan Shift
Charge relaxation in glass can be modeled using the so-called Tool–Narayanaswamy–
Moynihan shift factor which is given as
Q 1–
log T = ---- --- 1 -
- + ------------ – --------
RT Tf T ref
here, a base-e logarithm is used, Q is the activation energy (SI unit: J/mol), R is the
universal gas constant, T is the current temperature, Tref is a reference temperature,
is a dimensionless activation energy fraction, and Tf is the so-called fictive temperature.
The fictive temperature is given as the weighted average of partial fictive temperatures.
Tf = wi Tfi with wi = 1
i i
Here, wi are the weights and Tfi are the partial fictive temperatures. The partial fictive
temperatures are determined from a system of coupled ordinary differential equations
(ODEs) which follow Tool’s equation
d T – T fi
T = -----------------
d t fi T
0i
The Electric Currents Interface assumes that there is only one type of charge carrier.
When there are multiple charge carriers, it is often necessary to solve the transport
equation for each charge carrier as well as the Poisson’s equation with Electrostatics
interface.
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
– d V – J e = dQ j (3-9)
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.
E + J e + j = 0
D =
E + J e + = 0
t
– 0 V – P – V – J e = Q j (3-10)
t
In planar 2D the dynamic formulations also involves the thickness d in the z direction:
• 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 3-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
-------------- --------------
NF VT NR VT
–1
= ----------------------------------------------- 1 + 1 + 4I S ----------------------- + ------------------------
1 e – 1 e
f bq
v bc v be I A I KR A
2 1 – ---------- - – ---------- - KF
V V
AF 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
BR
v be v bc
I S -------------
N V
- --------------
N V
i ce = A ------- e F T + e C T
f bq
k B T NOM
V T = -----------------------
-
q
There are also two capacitances that use the same formula as the junction capacitance
of the diode model. In the parameter names below, replace x with C for the
base-collector capacitance and E for the base-emitter capacitance.
v bx – MJx
1 – ---------
-
V Jx v bx F C V Jx
C jbx
= AC Jx
v bx v bx F C V Jx
1 – F – 1 – M Jx 1 – F 1 + M + M ---------
-
C C Jx Jx V Jx
The following equations are used to compute the relations between currents and
voltages in the circuit.
C gd = C gd0 W
C gs = C gs0 W
–MJ
1 – v 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 JP
B
The following equations are used to compute the relations between currents and
voltages in the circuit.
vd –M
1 – ------
- vd FC VJ
V J
C j = C J0
– 1 – M vd
1 – FC 1 – F C 1 + M + M ------- v d F C V J
V J
k B T NOM
V T = -----------------------
-
q
This chapter describes the physics interfaces found under the Fluid Flow
branch ( ).
The list below details the documentation for the MEMS Module physics interfaces
on the Fluid Flow branch:
• The Thin-Film Flow Interfaces (in the CFD Module User’s Guide)
• Theory for the Thin-Film Flow Interfaces (in the CFD Module User’s Guide)
• The Fluid–Solid Interaction Interface (in the Structural Mechanics Module
User’s Guide)
• The Solid–Thin-Film Damping Interface (in the Structural Mechanics Module
User’s Guide)
• The Piezoelectricity Interface, Solid (in the Structural Mechanics Module
User’s Guide)
See The Single-Phase Flow, Laminar Flow Interface in the COMSOL Multiphysics
Reference Manual for other Fluid Flow interface and feature node settings.
181
182 | CHAPTER 4: FLUID FLOW INTERFACES
5
The MEMS module includes many of the features found in the Solid Mechanics
interface in the Structural Mechanics Module.
The list below details the documentation for the MEMS Module physics interfaces
on the structural mechanics, acoustics, and fluid flow branches in the Add Physics
dialog:
| 183
• The Elastic Waves, Time Explicit Interface (in the Elastic Waves Interfaces chapter
of the Acoustics Module User’s Guide)
• The Thermal Stress, Solid Interface (in the Multiphysics Interfaces and Couplings
chapter of the Structural Mechanics Module User’s Guide)
• The Piezoelectricity Interface, Solid (in the Multiphysics Interfaces and Couplings
chapter of the Structural Mechanics Module User’s Guide)
• The Electrostriction Interface (in the Multiphysics Interfaces and Couplings chapter
of the Structural Mechanics Module User’s Guide)
• The Ferroelectroelasticity Interface (in the Multiphysics Interfaces and Couplings
chapter of the Structural Mechanics Module User’s Guide)
• The Fluid–Solid Interaction Interface (in the Multiphysics Interfaces and Couplings
chapter of the Structural Mechanics Module User’s Guide)
• The Fluid–Solid Interaction, Fixed Geometry Interface (in the Multiphysics
Interfaces and Couplings chapter of the Structural Mechanics Module User’s
Guide)
• The Fluid–Solid Interaction, Conjugate Heat Transfer Interface (in the
Multiphysics Interfaces and Couplings chapter of the Structural Mechanics
Module User’s Guide)
• Hygroscopic Swelling Coupling (in the Multiphysics Interfaces and Couplings
chapter of the Structural Mechanics Module User’s Guide)
When this multiphysics interface is added using the Model Wizard, a Solid Mechanics
interface and an Electrostatics interfaces are added to the Model Builder. In addition, the
Multiphysics Couplings node is added, which automatically includes the multiphysics
coupling feature Electromechanical Forces.
The full functionality of the Solid Mechanics and Electrostatics interfaces is accessible
under the corresponding interfaces at the domain, boundary line, or point level in the
geometry. The participating Electrostatics interface by default adds the Charge
Conservation node with its selection set to all domains, and the material type set to
Solid. This node serves to represent solid dielectric material domains. Any number of
Charge Conservation nodes can be also added manually to the interface and configured
to represent either solid or nonsolid (for example, air or free space) materials.
The presence of such multiphysics node in the model will enforce the
Include geometric nonlinearity option on all applicable study nodes.
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers, and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is eme1.
DOMAIN SELECTION
The domain selection is locked, so that all applicable domains are selected
automatically. Such domains represent an intersection of the applicable domains under
the corresponding Electrostatics and Solid Mechanics interfaces selected in the coupling
feature.
• Charge Conservation, if its material type input is set to Solid. Use this domain feature
for solid dielectric materials.
• Charge Conservation, Ferroelectric
• Terminal. You can use the domain terminal feature to represent conductive solid
materials such as, for example, doped silicon.
In Solid Mechanics interface, the following domain material feature are applicable:
For nonsolid dielectric domains, remove them from the Solid Mechanics
interface selection but keep them selected in the Electrostatics interface
selection.
ELECTROMAGNETIC STRESS
To display this section, click the Show More Options button ( ) and select Advanced
Physics Options in the Show More Options dialog box.
Here, you can select the Electromagnetic Stress Tensor to be used in the coupling. The
following options are available:
• Minkowsky (default)
• Einstein-Laub
• Chu
COUPLED INTERFACES
This section defines the physics interfaces involved in this Electromechanical Forces
coupling. The Solid mechanics and Electrostatics lists include all applicable physics
interfaces.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(macOS and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when the physics interface is chosen in the Model Wizard
or Add Physics window, then the two participating physics interfaces are selected.
You can also select None from either list to uncouple the coupling from a physics
interface. If the physics interface is removed from the Model Builder, for example Solid
If a physics interface is deleted and then added to the model again, then
in order to reestablish the coupling, you need to choose the physics
interface again from the lists. This is applicable to all multiphysics
coupling nodes that would normally default to the once present physics
interface. See Multiphysics Modeling Workflow in the COMSOL
Multiphysics Reference Manual.
Deforming Domain
When an Electromechanics multiphysics interface is added using the Model Wizard, a
Deforming Domain and a Symmetry nodes are added automatically in Model Builder
under the Definitions node. These nodes serve to represent domains and boundaries
where the mesh can deform. By default, both nodes have empty selection. Under the
Deforming Domain, you can select those domain that represent nonsolid dielectric
materials such as air, fluid, and free space. However, this is only needed if the
boundaries of such domains experience significant changes due to the deformation or
rotation of the adjacent solid domains. Otherwise, the moving mesh computations
could just introduce unnecessary overhead and nonlinearity in the model (see
Deformed Mesh Fundamentals). Domains selected in the Solid Mechanics interface
are not applicable.
Under the Deforming Domain, you can also choose the Mesh smoothing type, by default
set to Hyperelastic. More information of the smoothing type can be found in the
Deforming Domain chapter in the COMSOL Multiphysics Reference Manual.
By default, the mesh is fixed at all external boundaries of the geometry. Use the
Symmetry node to select the external boundaries, where the mesh can slide along the
boundary. You can also add other types of boundary conditions for the mesh motion,
for details see Moving Mesh Features in the COMSOL Multiphysics Reference
Manual.
2
d u
2
= + fv
dt
where is the density, u denotes the structural displacement vector, is the stress
tensor, and fv is an external volume force (such as gravity).
In the stationary case, there is no acceleration, and the equation representing the force
balance is
0 = + fv
The stress tensor must be continuous across a stationary boundary between two
materials. This corresponds to the equation
2 – 1 n = 0 (5-1)
Material 2
Material 1
For more information on the mechanical stress tensor for elastic materials,
see the documentation for the physics interfaces. For example, Structural
Mechanics Interfaces.
Maxwell stresses exists in all dielectric materials, including air and even free space.
However, the force magnitude is rather small, and it usually can cause significant
deformations only at small dimension.
The electromagnetic stress tensor in a vacuum (in the absence of magnetic fields) is
given by (Ref. 1):
1
EM, V = 0 E E – --- 0 E E I (5-2)
2
Where E is the electric field, I is the identity tensor, 0is the permittivity of free space,
and
E E ij = E i E j
where the electric susceptibility can be a function of the mechanical strain in the
material
= 0 I – 2T + tr T I
1 T
T = --- u + u
2
1
EM = 0 r E E – --- 0 r E E I (5-3)
2
where the relative permittivity of the material with the absence of deformation is
introduced as r = 1 + 0 . It can be written equivalently as
1
EM = EM, V + 0 0 E E – --- E E I
2
where the two terms represent, respectively, the contributions from: the underlying
free space and the material polarization.
LARGE DEFORMATIONS
For finite deformations, the expressions for the electromagnetic stress and material
polarization can be derived using the following thermodynamic potential called
electric enthalpy:
1 –1
H eme = W s C – --- 0 r JC : E m E m
2
where the subscript m indicates that the vector components must be taken on the
material frame, and the right Cauchy-Green deformation tensor is
T
C = F F
H eme
S = 2
C
H eme
Dm = –
Em
The electromagnetic stress tensor can be used to compute the forces acting on a
dielectric body. From Equation 5-1 the balance of forces at the surface of a dielectric
body (material 1) in vacuum or air (material 2) implies:
2 – 1 n = 0
1 = EM + mech
2 = EM V – p 2 I
where m is a mechanical component of the total stress and the ambient pressure p2 in
to surrounding air or other medium (if present). Using Equation 5-2:
1 n = – p 2 n – --- 0 E E n + 0 n E E
1
2
1
EM = D E – --- D E I
2
1
f = EM = e E – --- E E :
2
There are two most often used alternatives. The first one is called Einstein–Laub stress
tensor:
1
EM = D E – --- 0 E E I
2
This form is widely accepted in modern electroelasticity and material science, see
Ref. 4. The corresponding body force can be written as:
f = EM = e E + P E
which is also called Kelvin force. Note that it is nonzero as soon as there are electric
field variation and polarization within the material:
P = 0 E
1
EM = 0 E E – --- 0 E E I (5-4)
2
f = EM = e E – P E
–1
D m = 0 JC E m + 0 m E m (5-5)
1
H eme = W s C – --- D m E m
2
Then,
H eme
Dm = –
Em
H eme –1 –T
S EM = 2 = JF EM F (5-6)
C
where EM is given by Equation 5-4, which shows the consistency of the theory. The
key assumption behind the derivation of Equation 5-6 is that the material components
of the electric susceptibilitym are independent of deformations.
When this multiphysics interface is added using the Model Wizard, a Solid Mechanics
interface and an Electrostatics, Boundary Elements interface are added to the Model
Builder. In addition, the Multiphysics Couplings node is added, which automatically
includes the multiphysics coupling feature Electromechanical Interface.
Electromechanical Interface
Use an Electromechanical Interface ( ) boundary multiphysics coupling to model
electromechanical interaction between deformable solid domains and nonsolid voids,
where the electric field is computed using the boundary element method. The
backward coupling to Electrostatics is due to the deformations of the boundaries.
The presence of such multiphysics node in the model will enforce the
Include geometric nonlinearity option on all applicable study nodes.
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers, and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is eme1.
For nonsolid dielectric domains, remove them from the Solid Mechanics
interface selection but keep them selected in the Electrostatics, Boundary
Elements interface selection.
The coupling is a boundary load on a structure caused by the Maxwell Stress at the
selected boundaries. The backward coupling to Electrostatics is due to the change of
the spatial position of the boundaries caused by the deformation of the structure.
COUPLED INTERFACES
This section defines the physics interfaces involved in this Electromechanical Forces
coupling. The Solid mechanics and Electrostatics, Boundary Elements lists include all
applicable physics interfaces.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(Mac and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when the physics interface is chosen in the Model Wizard
or Add Physics window, then the two participating physics interfaces are selected.
You can also select None from either list to uncouple the coupling from a physics
interface. If the physics interface is removed from the Model Builder, for example Solid
If a physics interface is deleted and then added to the model again, then
in order to reestablish the coupling, you need to choose the physics
interface again from the lists. This is applicable to all multiphysics
coupling nodes that would normally default to the once present physics
interface. See Multiphysics Modeling Workflow in the COMSOL
Multiphysics Reference Manual.
The Piezoelectricity interface is available for 3D, planar 2D, and axisymmetric 2D
geometries.
When this multiphysics interface is added using the Model Wizard, an Electrostatics
interface and a Heat Transfer in Solids interface are added to the Model Builder.
The Heat Transfer in Solids will use its default domain feature, Solid, where you enter
the thermal properties of the material. All heat transfer functionality for modeling is
also accessible to include surrounding solid or nonsolid domains. To model heat
transfer by convection in nonsolid domains, you can add a Fluid domain feature.
You enter the pyroelectric coupling data on the Pyroelectricity multiphysics coupling
node.
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers, and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is pye1.
COUPLED INTERFACES
This section defines the physics involved in the Pyroelectricity multiphysics coupling.
The Electrostatics and Heat Transfer drop-down menu lists include all applicable physics
interfaces.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(macOS and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when a Pyroelectricity multiphysics interface has been
selected in the Model Wizard or Add Physics window, then the participating
Electrostatics and Heat Transfer in Solids interfaces are selected.
You can also select None from either list to uncouple the Pyroelectricity node from a
physics interface. If the physics interface is removed from the Model Builder, for
If a physics interface is deleted and then added to the model again, then
in order to reestablish the coupling, you need to choose the physics
interface again from the lists. This is applicable to all multiphysics
coupling nodes that would normally default to the once present interface.
See Multiphysics Modeling Workflow in the COMSOL Multiphysics
Reference Manual.
DOMAIN SELECTION
The domain selection is set by default to all domains, so that all applicable domains are
selected automatically. Such domains represent an intersection of the applicable
domains under the corresponding Electrostatics and Heat Transfer in Solids interfaces
selected in the coupling feature.
• Charge Conservation, if its material type input is set to Solid. Use this domain feature
for solid dielectric materials, for which a linear dependency can be assumed for the
electric polarization with respect to the applied electric field.
• Charge Conservation, Piezoelectric. Use this domain feature for solid linear
piezoelectric materials.
• Charge Conservation, Ferroelectric. Use this domain feature for solid ferroelectric or
nonlinear piezoelectric materials.
• Solid. Use this domain feature for solid pyroelectric and dielectric materials.
To model solid nonpyroelectric domain, remove such domains from the coupling
feature selection.
COUPLING TYPE
From the list, choose one of these coupling types:
• Pyroelectric effect to include only the change in the material polarization caused by
temperature variations.
• Electrocaloric effect, also known as inverse pyroelectric effect, to include only the
heat source due to the polarization variation in time.
Q p = – T ----- p ET E
t
• Fully coupled (default) to include both the direct and inverse pyroelectric effects.
COUPLING SETTINGS
You enter the Total pyroelectric coefficient pET which is a vector (SI unit: C/m2/K)
specified in the selected coordinate system. The default is to take the values From
material.
You also can enter the Reference temperature Tref with the default value of 293.15 K.
p ET = p ES + e ES
The Piezoelectricity and Pyroelectricity interface is available for 3D, planar 2D, and
axisymmetric 2D geometries.
When this multiphysics interface is added using the Model Wizard, a Solid Mechanics
interface, an Electrostatics interface, and a Heat Transfer in Solids interface are added to
the Model Builder.
The participating Solid Mechanics interface includes the default Piezoelectric Material
feature with its selection set to all domains. The Electrostatics interface has a default
Charge Conservation, Piezoelectric feature with similar settings. Such features can be also
added manually to their corresponding interfaces similar to any other material model
therein. The Piezoelectricity multiphysics coupling node can be active only on the
selection, where both features Piezoelectric Material and Charge Conservation,
Piezoelectric are active. You specify the piezoelectric coupling data together with both
the structural and electrical material data under the Piezoelectric Material node.
All solid mechanics and electrostatics functionality for modeling is also accessible to
include surrounding elastic solids or air domains. For example, add any solid
mechanics material for other solid domain, a dielectric model for air (via Charge
Conservation feature), or a combination. Note that in order to model a nonsolid
dielectric domain, you need to remove such domain from the domain selection for the
entire Solid Mechanics interface. This is because all material models under that
interface represent solid materials (with the Linear Elastic Material node being always
present and active in all those domains, where it is not explicitly overridden by any
other material model).
You enter the pyroelectric coupling data on the Pyroelectricity multiphysics coupling
node, and the thermal expansion data on the Thermal Expansion multiphysics
coupling node.
When this multiphysics interface is added, the Model Builder tree structure is
populated with two physics interfaces (Solid Mechanics and Heat Transfer in Solids), and
the Thermal Expansion multiphysics coupling.
In the Settings window for the Solid Mechanics interface, the Structural Transient
Behavior is set to Include inertial terms by default. In the Settings window for Thermal
Expansion, under the section Heat Sources, the check box for Thermoelastic damping is
enabled by default.
The Heat Transfer in Solids interface provides features for modeling heat transfer by
conduction, convection, and radiation. A Solid model is active by default on all
domains. All functionality for including other domain types, such as a fluid domain, is
also available. The temperature equation defined in solid domains corresponds to the
differential form of the Fourier’s law that may contain additional contributions like
heat sources.
On the Couplings
The coupling occurs on the domain level, where two physical effects occur:
• The temperature from the Heat Transfer interface acts as a thermal load for the Solid
Mechanics interface, causing thermal expansion.
• The tension and compression from the Solid Mechanics interface acts as a heat
source for the Heat Transfer interface, causing local cooling and heating.
COUPLING FEATURES
When physics interfaces are added from the Model Wizard or using the Add Physics
window then the appropriate coupling feature is automatically added to the physics
interface. However, the physics interfaces can be added one at a time, followed by the
coupling features.
Coupling features are available from the context menu (right-click the
Multiphysics node) or from the Physics toolbar, Multiphysics Couplings
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
Thermal Expansion node is described in the section The Thermal Stress, Solid
Interface in the Structural Mechanics User’s Guide.
The equations of thermoelasticity are derived from the first law of thermodynamics,
which can be stated as follows:
where dU is the change in internal energy, dQ' is the heat flow into the system (the
prime indicates an inexact differential in this case) and dW' is the work done on the
system. For a small part of a solid (sufficiently small that the stresses and strains are
uniform), with an initial reference density, 0, the first law can be rewritten in the
following form (assuming that the differential changes occur between equilibrium
states):
where Ta is the absolute temperature, s is the entropy per unit mass, is the elastic
part of the second Piola-Kirchhoff stress (in general a rank 2 tensor), is the material
strain (also a tensor). In general the second Piola-Kirchhoff stress tensor, p, must be
split into elastic () and inelastic () parts such that:
p = +
The elastic part of the stress tensor, , does work :d during a change in the strain.
The inelastic part of the stress tensor, , generates heat at a rate :(d/d) when the
strain is changing and is identified with internal or material damping. These internal
damping mechanisms are associated with microscopic phenomena such as dislocation
movement.
From Equation 5-8 it is possible to make the following identifications for Ta and
u u
Ta = = 0
S S
Next the entropy balance equation must be derived. Since thermoelasticity involves
irreversible processes, the assumption of equilibrium required to derive Equation 5-8
is no longer valid. Instead an assumption of “local” equilibrium is made. It is assumed
that although the system is not in equilibrium, there exists within small elements a state
of local equilibrium, for which the local entropy per unit mass, s, is the same function
of the internal energy, strain, and particle number as it was in equilibrium. This
assumption is commonly employed in the modeling of transport phenomena and is
justified only by the validity of conclusions derived from it and by results obtained from
specific microscopic models, for near-equilibrium situations. For a small volume
element in the material frame Equation 5-8 can then be written as:
1 1
0 ds = ------- 0 du – ------- : d
Ta Ta
ds 1 du 1 d
0 ------ = 0 ------- ------- – ------- : ------ (5-9)
dt T a dt T a dt
From the first law (Equation 5-7) the rate of change of internal energy is given by:
where w is the work done per unit volume and q is the heat accumulated per unit
volume. The heat accumulated can be written as the sum of the heat sources and the
divergence in the material frame heat flux:
dq d
------- = – q + Q + :
dt dt
where Q represents the heat source per unit volume and is the inelastic part of the
stress tensor. The rate of doing work (per unit reference volume) by a linear elastic
material is given by the elastic part of the second Piola–Kirchhoff stress contracted with
the rate of material strain. Per unit volume the following equation is obtained:
dw d
-------- = : ------
dt dt
ds 1 1 1 d
0 ------ = – ------- q + ------- Q + ------- : (5-10)
dt Ta Ta Ta d t
q = – T a
ds d
T a 0 ------ = T a + Q + : (5-11)
dt dt
It is now necessary to derive an expression for the rate of change of entropy with
respect to time. In order to do this an assumption of local equilibrium is used once
again. Using Equation 5-8 the equation is written:
which defines a new thermodynamic potential, the Gibbs free energy per unit mass,
given by:
Changes in the Gibbs free energy per unit mass take the form:
1
dg = – sdT a – ------ :d
0
g g
s = – = –0
T a T a
By differentiating each of the above equations a second time, it is possible to derive the
following Maxwell relation:
2
s = -----
1- g
(5-12)
T a 0 T a = – T
a
It is now possible to derive an expression for the entropy of the solid. Assuming that
the elastic stress is an invertible function of the strain, we can write s=s(,Ta). Thus:
s s
ds = :d + dT
T a T a a
1 s
ds = ------ :d + dT
0 T a T a a
Therefore:
1 d s dT
------ = ------ : ------- +
ds ----------a-
dt 0 T a dt T a dt
By definition the heat capacity of the solid at constant stress is given by:
q s
cp = = Ta
T T
Thus:
Substituting Equation 5-13 into Equation 5-11 gives the following equation for
thermoelasticity:
dT a ------
0 c p ----------- = T a + Q + : – T a
d d
: - (5-14)
dt dt T a dt
An additional heat source term is present in Equation 5-14, compared to the standard
heat transfer equations in solids. This term couples the structural problem with the
heat transfer problem. In turn the heat transfer equation couples back into the
structural problem through the constitutive relationship.
In the particular case of a linear elastic material (in the absence of damping) the stress
and strain are related by Duhamel–Hooke’s law:
– i = C: – i – T a – T ref (5-15)
where C is the elasticity tensor, i is the initial stress, i is the initial strain, and Tref is
the reference temperature at which the strain and stresses take the initial values.
This equation couples the heat transfer equation to the structural problem. Given a
temperature independent thermal expansivity, and no material damping,
Equation 5-14 takes the form:
dT a d
0 c p ----------- = T a + Q – T a : -------
dt dt
This is the usual form of the equation for linear thermoelasticity. In the
COMSOL Multiphysics Material Library, a slightly unusual definition of
the thermal expansivity is used, so that is not temperature independent.
COMSOL Multiphysics therefore implements a version of
Equation 5-14, which automatically includes these effects. In brief, an
expression for is obtained by solving Equation 5-15 and the expression
is differentiated with respect to Ta while keeping constant.
2
u = u 0 + au 1 + a u 2 +
where a is a constant that is subsequently set to 1 and that is used to keep track of the
order of the term in the expansion.
The corresponding temperatures, stresses, and strains in the device are given by:
2
T a = T 0 + aT 1 + a T 2 +
2
= 0 + a 1 + a 2 +
2
= 0 + a 1 + a 2 +
2
= a1 1 + a 2 +
where a is a constant that is subsequently set to 1. a can be used to track the order of
the deviation from T0.
Substituting these relations into Equation 5-14 and keeping terms up to the second
order gives:
dT 1 2 dT 2 2
a 0 c p ---------- + a 0 c p ---------- = T 0 + a T 1 + a T 2 + Q 0
dt dt
d 1 ---------
d 1
– aT 0
2
+ a 1 : : -
dt T a dt
---------
d 1 d
–a T1 : - – a T0
2 2 2
: ----------
T a dt T a dt
2 ---------
d 1
–a T0 T1 : -
T a T a dt
Note that the heat source term, Q, is assumed to be time independent (Q=Q0). The
final term in the above equation results from consideration of the Taylor expansion of
the strain at constant stress, which to first order can be represented by a Taylor series
of the form:
Consequently:
-
--------- 0 0
T a = ---------- + aT 1 ---------- ---------- +
T a T a T a
0 – i = C: 0 – i – T 0 – T ref
2 2 2
a 1 + a 2 = C: a 1 + a 2 – aT 1 – a T 2
2 2 2
a 1 + a 2 = iC: a 1 + a 2 – aT 1 – a T 2
where the material damping has been introduced into the equation system by means
of an anisotropic loss factor, .
Grouping together terms of the same order in a gives the following equations:
– T 0 = Q 0
dT 1 0 d 1
a 0 c p ---------- = a T 1 – aT 0 : ----------
dt T a dt
2 dT 2 2 2 d 1 2 0 d 1
a 0 c p ---------- = a T 2 + a 1 : – a T1 : ----------
dt dt T a dt
2 0 d 1 2 0 d 2
–a T0 T1 : ---------- – a T 0 : ----------
T a T a dt T a dt
dT 0 d
0 c p -------- = T – T 0 : ------- (5-16)
dt T a dt
Where a=1 and the subscript has been dropped from the dependent variable T
(previously T1) and the stress deviation (previously 1). Equation 5-16 is solved by
the Thermoelasticity interface. Writing the equation in the frequency domain gives:
This is the frequency domain form of the equations solved by COMSOL Multiphysics.
Considering the quadratic terms in Equation 5-12, there are four second-order heat
source terms in the equation system:
d 1
Q visc = 1 :
dt
0 d 1
Q TED1 = – T 1 : ----------
T dt
0 ---------
d 1
Q TED2 = – T 0 T 1 : -
T T dt
0 d 2
Q TED3 = – T 0 : ----------
T dt
In the frequency domain, the first three terms involve the products of small harmonic
deviations with other similar terms (that is, terms that vary in time in a similar manner
to cos2 t. The trigonometric identity cos2 t=(1cos 2 t)/2 makes it explicit that
these terms have a time-independent heating effect in addition to the heating that
occurs at double the frequency of the driving signal). These second-order terms have
a small effect on the time-varying solution and lead to a second harmonic term in the
frequency domain. However, if the time average is nonzero, this can, over many
oscillatory periods, result in a change in the background temperature. In COMSOL
Multiphysics the (nonzero) time-averaged heat sources are available as variables, so
that the self-heating of the structure can be computed.
First note that the time average of the term QTED3 must be zero, since although the
second-order stress term 2 has a nonzero time average, its time derivative cannot. The
other terms in the expression for QTED3 are all time independent. Each of the
remaining terms can have a nonzero time average as these involve the multiplication
of two separate first-order terms (potentially with phase differences between them).
The following result for the multiplication of two general oscillatory complex
numbers, a1ei(t+1) and, a2ei(t+2) are useful:
The time average heat sources above can be written in the following form in the
frequency domain:
1
Q visc = --- Re : i
2
1 0
Q TED1 = – --- Re T : i
2 T
1 0
Q TED2 = – --- Re T 0 T : i
2 T T
The derivation of the last term <QTED2> strongly relies on the prerequisite that the
temperature dependence of the thermal expansion coefficient alpha is
thermodynamically consistent. Since this prerequisite is usually not satisfied due to
measurement or formulation errors, we have decided to remove the variable from the
physics interface.
COMSOL Multiphysics also defines a total heat source term (solid.Qh) of the form:
Q h = Q visc + Q TED1
The realdot operator is used to evaluate these quantities in a manner that produces
the correct terms in the Jacobian for sensitivity analyses. To compute the slow time
dependence of the background temperature, the steady-state equation can be modified
dT 0
0 cp = T 0 + Q h + Q 0
dt
i 0 c p T = T – iT 0 :
T
With this formulation, the transient part of the problem can be solved on
a time scale that is much longer than the period of the oscillations.
Piezoresistivity Interfaces
This chapter describes the physics interfaces found under the Structural
Mechanics>Piezoresistivity subbranch ( ).
In this chapter:
| 217
T he P i e z o r e s i st i v i t y, D omai n Cu rren t s
Interface
The Piezoresistivity, Domain Currents predefined multiphysics coupling interface ( ),
found under the Structural Mechanics>Piezoresistivity branch ( ) when adding a
physics interface, combines the Solid Mechanics and Electric Currents interfaces to
simulate the piezoresistive effect. This multiphysics interface is appropriate for
situations when the thicknesses of the conducting and piezoresistive layers are resolved
by the mesh.
When this multiphysics interface is added using the Model Wizard, a Solid Mechanics
interface and an Electric Currents interface are added to the Model Builder.
In addition, the Multiphysics Couplings node is added, which automatically includes the
multiphysics coupling feature Piezoresistivity, Domain Currents.
The participating Electric Currents interface by default adds the Current Conservation,
Piezoresistive node with its selection set to all domains. This node can also be added
manually to the interface similar to any other material model therein. All
piezoresistivity properties are entered in this node.
The Piezoresistivity, Domain Currents multiphysics coupling node can be active only on
the selection where both features Current Conservation, Piezoresistive and Linear Elastic
Material (or Hyperelastic Material) are active.
For example, if both Solid Mechanics and Electric Currents interfaces are added,
COMSOL Multiphysics adds an empty Multiphysics Couplings node. You can choose
the available coupling feature Piezoresistivity, Domain Currents, but the modified
settings are not included. In this example, the Current Conservation, Piezoresistive node
needs to be added manually.
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.
Coordinate Systems
When setting up a model, the coordinate systems used for the participating physics
interfaces should be consistent. In particular, the selected item for the Coordinate
system menu under the Solid Mechanics>Linear Elastic Material>Coordinate System
Selection section should match the one for the Coordinate system menu under the
Electric Currents>Current Conservation, Piezoresistive>Coordinate System Selection
section.
Coupling Feature
The following multiphysics coupling feature node is described in this section:
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers, and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is pzrd1.
DOMAIN SELECTION
Only domains that have both Current Conservation, Piezoresistive selected in the Electric
Currents interface and Linear Elastic Material selected in the Solid Mechanics interface are
applicable.
COUPLED INTERFACES
This section defines the physics involved in the Piezoresistivity, Domain Currents
multiphysics coupling. The Solid mechanics and Electric Currents, shell lists include all
applicable physics interfaces.
The default values depend on how the Piezoresistivity, Domain Currents node is created.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(Mac and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when a Piezoresistivity, Domain Currents multiphysics
coupling interface is selected in the Model Wizard or Add Physics window, then the
participating Solid Mechanics and Electric Currents interfaces are selected.
You can also select None from either list to uncouple the Piezoresistivity, Domain
Currents node from a physics interface. If the physics interface is removed from the
If a physics interface is deleted and then added to the model again, then
in order to reestablish the coupling, you need to choose the physics
interface again from the lists. This is applicable to all multiphysics
coupling nodes that would normally default to the once present interface.
See Multiphysics Modeling Workflow in the COMSOL Multiphysics
Reference Manual.
When this multiphysics interface is added using the Model Wizard, a Solid Mechanics
interface and an Electric Currents in Shells interface are added to the Model Builder.
In addition, the Multiphysics Couplings node is added, which automatically includes the
multiphysics coupling feature Piezoresistivity, Boundary Currents.
The participating Electric Currents in Shells interface by default adds the Piezoresistive
Shell node with its selection set to all domains. This node can also be added manually
to the interface similar to any other material model therein. All piezoresistivity
properties are entered in this node.
The Piezoresistivity, Boundary Currents multiphysics coupling node can be active only
on the selection where both features Piezoresistive Shell and Linear Elastic Material (or
Hyperelastic Material) are active.
For example, if both Solid Mechanics and Electric Currents in Shells interfaces are added,
COMSOL Multiphysics adds an empty Multiphysics Couplings node. You can choose
the available coupling feature Piezoresistivity, Boundary Currents but the modified
settings are not included. In this example the Piezoresistive Shell node needs to be
added manually.
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.
Coordinate Systems
When setting up a model, the coordinate systems used for the participating physics
interfaces should be consistent. In particular, the selected item for the Coordinate
system menu under the Solid Mechanics>Linear Elastic Material>Coordinate System
Selection section should match the one for the Coordinate system menu under the
Electric Currents in Shells>Piezoresistive Shell>Coordinate System Selection section.
Coupling Feature
The following multiphysics coupling feature node is described in this section:
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers, and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is pzrb1.
BOUNDARY SELECTION
Only boundaries that have both Piezoresistive Shell selected in the Electric Currents in
Shells interface and are adjacent to domains with Linear Elastic Material selected in the
Solid Mechanics interface are applicable.
COUPLED INTERFACES
This section defines the physics involved in the Piezoresistivity, Boundary Currents
multiphysics coupling. The Solid mechanics and Electric Currents, shell lists include all
applicable physics interfaces.
The default values depend on how the Piezoresistivity, Boundary Currents node is
created.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(Mac and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when a Piezoresistivity, Boundary Currents multiphysics
coupling interface is selected in the Model Wizard or Add Physics window, then the
participating Solid Mechanics and Electric Currents, Shell interfaces are selected.
You can also select None from either list to uncouple the Piezoresistivity, Boundary
Currents node from a physics interface. If the physics interface is removed from the
If a physics interface is deleted and then added to the model again, then
in order to reestablish the coupling, you need to choose the physics
interface again from the lists. This is applicable to all multiphysics
coupling nodes that would normally default to the once present interface.
See Multiphysics Modeling Workflow in the COMSOL Multiphysics
Reference Manual.
When this multiphysics interface is added using the Model Wizard, a Shell interface and
an Electric Currents, Shell interface are added to the Model Builder.
In addition, the Multiphysics Couplings node is added, which automatically includes the
multiphysics coupling feature Piezoresistivity, Shell.
The participating Electric Currents in Shells interface by default adds the Piezoresistive
Shell node with its selection set to all domains. (Such a node can also be added
manually to the corresponding interface.) All piezoresistivity properties are entered in
this node.
The Piezoresistivity, Shell multiphysics coupling node can be active only on the
selection where both a Piezoresistive Shell and a Linear Elastic Material feature are active.
For example, if both Shell and Electric Currents in Shells interfaces are added, COMSOL
adds an empty Multiphysics Couplings node. You can choose the available coupling
feature Piezoresistivity, Shell but the modified settings are not included. In this example
the Piezoresistive Shell node needs to be added manually.
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.
Coordinate Systems
When setting up a model, the coordinate systems used for the participating physics
interfaces should be consistent. In particular, the selected item for the Coordinate
system menu under the Shell>Linear Elastic Material>Shell Local System>Coordinate
System Selection section should match the one for the Coordinate system menu under
the Electric Currents, Shell>Piezoresistive Shell>Coordinate System Selection section.
Coupling Feature
The following multiphysics coupling feature node is described in this section:
• Piezoresistivity, Shell
SETTINGS
The Label is the default multiphysics coupling feature name.
The Name is used primarily as a scope prefix for variables defined by the coupling node.
Refer to such variables in expressions using the pattern <name>.<variable_name>. In
order to distinguish between variables belonging to different coupling nodes or physics
interfaces, the name string must be unique. Only letters, numbers, and underscores (_)
are permitted in the Name field. The first character must be a letter.
The default Name (for the first multiphysics coupling feature in the model) is pzrs1.
BOUNDARY SELECTION
Only boundaries that have both Piezoresistive Shell selected in the Electric Currents in
Shells interface and Linear Elastic Material selected in the Shell interface are applicable.
COUPLED INTERFACES
This section defines the physics involved in the Piezoresistivity, Shell multiphysics
coupling. The Shell and Electric currents, shell lists include all applicable physics
interfaces.
The default values depend on how the Piezoresistivity, Shell coupling node is created.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(Mac and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when a Piezoresistivity, Shell multiphysics coupling
interface is selected in the Model Wizard or Add Physics window, then the newly
added Shell and Electric Currents in Shells interfaces are selected.
You can also select None from either list to uncouple the Piezoresistivity, Shell node
from a physics interface. If the physics interface is removed from the Model Builder, for
example if Shell is deleted, then the list defaults to None as there is nothing to couple to.
The Piezoresistivity, Layered Shell interface ( ) combines structural layered shells and
electric currents in shells together with the constitutive relationships required to model
the piezoresistive effects. The piezoresistive coupling can be formulated using either
the piezoresistance or elastoresistance form.
When a predefined Piezoresistivity, Layered Shell interface is added from the Structural
Mechanics branch ( ) of the Model Wizard or from Add Physics windows, Layered Shell
and Electric Currents in Layered Shells interfaces are added to the Model Builder.
The participating Electric Currents in Layered Shells interface will include by default a
Piezoresistive Shell feature. (Such a feature can be also added manually to the
corresponding interface.) You this feature to select those boundaries and layers, where
the material is supposed to experience piezoresistive coupling. All piezoresistivity
properties are entered in this feature. Note that the Piezoresistivity, Layered
multiphysics coupling node can be active only on the selection, where Piezoresistive
Shell is active. When used without the coupling feature, the Piezoresistive Shell node
acts as an ordinary Conductive Shell feature with its material data input limited to the
electric conductivity only.
All structural and electric functionality for modeling is also accessible to include
surrounding elastic solids or air domains. For example, add any material model for
other solid domains, a conductive or dielectric model, or a combination of these. Note
that in order to model a nonsolid conductive or dielectric domain, such as air, you need
to remove such domains from the domain selection for the entire Layered Shell
interface. This is because all material models under that interface represent solid
materials (with the Linear Elastic Material node always being present and active in all
domains where it is not explicitly overridden by any other material model).
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 both Layered Shell and Electric Currents in Layered Shells interfaces are
added, COMSOL Multiphysics adds an empty Multiphysics node. You can add the
available coupling feature Piezoresistivity, Layered but the modified settings are not
included.
Coupling features are available from the context menu (right-click the
Multiphysics node) or from the Physics toolbar, Multiphysics Couplings
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.
Piezoresistivity, Layered
The Piezoresistivity, Layered multiphysics coupling node ( ) computes the
appropriate conductivity tensor according to the mechanical stress or strain from the
Linear Elastic Material node in the Layered Shell interface and the piezoresistivity
properties from the Piezoresistive Shell node in the Electric Currents in Layered Shells
interface.
SETTINGS
The Label is the default multiphysics coupling feature name.
The default Name (for the first multiphysics coupling feature in the model) is pzrl1.
COUPLED INTERFACES
This section defines the physics involved in the Piezoresistivity, Layered Shell
multiphysics coupling. The Structure and Electric currents in shells lists include all
applicable physics interfaces.
The default values depend on how the Piezoresistivity, Layered coupling node is
created.
• If it is added from the Physics ribbon (Windows users), Physics contextual toolbar
(Mac and Linux users), or context menu (all users), then the first physics interface
of each type in the component is selected as the default.
• If it is added automatically when a Piezoresistivity, Layered Shell multiphysics
coupling interface is selected in the Model Wizard or Add Physics window, then the
newly added Layered Shell and Electric Currents in Layered Shells interfaces are
selected.
You can also select None from either list to uncouple the Piezoresistivity, Layered node
from a physics interface. If the physics interface is removed from the Model Builder, for
example if Layered Shell is deleted, then the list defaults to None as there is nothing to
couple to.
If a physics interface is deleted and then added to the model again, then
in order to reestablish the coupling, you need to choose the physics
interface again from the lists. This is applicable to all multiphysics
coupling nodes that would normally default to the once present interface.
See Multiphysics Modeling Workflow in the COMSOL Multiphysics
Reference Manual.
• Piezoresistive Materials
• Electric Currents Formulation
• Tensor Versus Matrix Formulations
• Material Property Definitions
• Coordinate Systems and Material Orientation
• References for the Piezoresistivity Interfaces
The Piezoresistivity, Shell interface uses features from the Shell interface,
which requires the Structural Mechanics Module. Both the theory and
features are described in the Structural Mechanics Module User’s Guide.
Piezoresistive Materials
A piezoresistive material responds to an applied stress with a change in its resistivity.
Normally piezoresistive materials are semiconductors, and the effect is associated with
band structure changes that alter the carrier mobility and number density much more
significantly than an applied stress would change, for example, the resistance of a
similar metal. The equivalent elastoresistance effect relates the strain in the material to
a change in resistivity. Note that the dependence between the structural mechanics and
the electrical properties of the material is unidirectional, so that applied currents do not
induce changes in the material stress or strain (assuming that other effects, such as
heating of the material, are negligible).
E = J + J (6-1)
where is the resistivity and is the induced change in the resistivity. In the general
case both and are rank 2 tensors (matrices). The change in resistance is related to
the stress S (for the piezoresistance form of the equations), or the strain (for the
elastoresistance form of the equations), by the constitutive relationship:
= S piezoresistance form
(6-2)
= M elastoresistance form
where is the piezoresistance tensor (SI unit: mPa1) and M is the elastoresistance
tensor (SI unit: m). Both of these quantities are material properties. and M are in
this case rank 4 tensors; however, they can be represented as matrices if the resistivity,
stress, and strain are converted to vectors within a reduced subscript notation.
J – J e + D = Qj
t
where J is the induced electric current, Je includes external current sources that do
not result from the material conductivity (such as thermoelectric currents), D is the
electric displacement field, and Qj is a current source term. Both Qj and Je are usually
zero for piezoresistive devices.
In the stationary or quasistatic case, the electric current is related to the electric field,
E, by the constitutive relation:
J = c E (6-3)
– c V – J e + 0 r V = Q j (6-4)
t
If current flow is restricted to a thin layer on the surface of a material this equation
takes the form:
– t d s c t V – J e + 0 r t V = Q j (6-5)
t
where the operator t represents the tangential derivative along the thin layer surface.
Equation 6-4 and Equation 6-5 form the basis of the Electric Currents interface.
–1 –1
J = c + E (6-6)
Although this expression appears cumbersome, it is required in this form for the
general case where the conductivity and change in resistivity are nonsymmetric tensor
quantities. For convenience, the tensor (c, eff) is defined as:
–1 –1
c ,eff = c + (6-7)
– c ,eff V – J e + 0 r V = Q j
t
while The Piezoresistivity, Boundary Currents Interface and The Piezoresistivity, Shell
Interface solve a modified form of equation Equation 6-5 on the boundary level:
– t d s c ,eff t V – J e + 0 r t V = Q j
t
j j kl
E i = ij J + ijkl J S piezoresistance form
(6-8)
j kl j
E i = ij J + M ij J kl elastoresistance form
Here the Einstein summation convention is used. Note that this equation implies that
ij=ijklSkl and ij=Mijklkl. Tensors and M have a number of symmetry
properties. Tensors S, , and are symmetric and can be represented as six
component vectors. There are two conventions for ordering the components in these
vectors, and the Piezoresistivity interfaces all use the Voigt notation prevalent in the
literature.
Cubic materials, such as silicon, possess symmetries that mean the two notations are,
in any case, equivalent when defining material properties. Within Voigt notation, the
vectors that correspond to S, , and are:
This notation is known as the reduced subscript notation, since the number of
subscripts required to represent the stress, strain, or change in resistivity is reduced
from two to one. The relationship between the stress or strain and the resistivity
change becomes:
r = s piezoresistance form
r = m e elastoresistance form
Where and m are now matrices related to the components of the corresponding
and M tensors in the following manner:
11 22 33 23 13 12
M 11 M 11 M 11 M 11 M 11 M 11
11 22 33 23 13 12
M 22 M 22 M 22 M 22 M 22 M 22
11 22 33 23 13 12
M 33 M 33 M 33 M 33 M 33 M 33
m =
11 22 33 23 13 12
M 23 M 23 M 23 M 23 M 23 M 23
11 22 33 23 13 12
M 13 M 13 M 13 M 13 M 13 M 13
11 22 33 23 13 12
M 12 M 12 M 12 M 12 M 12 M 12
For the most general triclinic crystal groups, the and m matrices are not
necessarily symmetric, as is the case for the elasticity matrix in structural
mechanics.
This is because the rank four tensors and M lack one of the symmetry properties of
the elasticity tensor (ijkl klij and Mijkl Mklij see Ref. 1 and Ref. 2 for details).
Finally the relationship between and m is given by:
m = D (6-9)
where D is the elasticity matrix (s=De), which is defined in Voigt notation as:
= 0 '
(6-10)
m = 0 m'
where 0 is the scalar resistivity of the undeformed material. Here has the units Pa1
and m is dimensionless. For comparison has the units mPa1 and m has units m.
Reduced subscript matrices and vectors do not obey the usual tensor
transformation rules for coordinate transforms. The COMSOL
Multiphysics software uses the quadratic transformation laws to transform
these quantities between coordinate systems where this is necessary. See
Ref. 1 for further details.
The local coordinate system for a material is set within the Linear Elastic Material and
Charge Conservation, Piezoresistive features (or their equivalents). These features
default to using the global system, but alternative systems can be defined by adding
one of the coordinate system nodes to the Definitions node level of the Model Builder.
Often the Rotated System feature is used to specify the new coordinate system in terms
of Euler angles.
3. M. Bao, Analysis and Design Principles of MEMS Devices, Elsevier B. V., 2005.
5. A. A. Barlian; W.-T. Park; J.R. Mallon, Jr.; A.J. Rastegar; and B.L. Pruitt, “Review:
Semiconductor Piezoresistance for Microsystems”, Proceedings of the IEEE, vol. 97,
no. 3, 2009.
Glossary
243
Glossary of Terms
accelerometer A sensor that measures acceleration or gravitational force. See also
sensor.
aspect ratio The ratio of the dimensions of a device in different directions, for
example, the ratio of height to width. See also lateral aspect ratio and vertical aspect
ratio.
biosensor A general term for sensor devices that either detect biological substances or
use antibodies, enzymes, or other biological molecules in their operation. Biosensors
are a subcategory of chemical sensors.
bonding A process by which one substrate is firmly attached to the surface of another.
cantilever beam A beam with one end fixed and one end free.
capacitive sensor A sensor that produces a signal due to a change in its capacitance.
Cauchy stress The most fundamental stress measure defined as force/deformed area
in fixed directions not following the body.
compliance matrix The inverse of the elasticity matrix. See elasticity matrix.
= D
Eulerian Model described and solved in a coordinate system that is fixed (spatial). See
also Lagrangian and arbitrary Lagrangian–Eulerian method.
initial strain The strain in a stress-free structure before it is loaded. See also strain and
residual strain.
initial stress The stress in a nondeformed structure before it is loaded. See also stress
and residual stress.
Kn = ----
L
Lagrangian Model described and solved in a coordinate system that moves with the
material. See also Eulerian and arbitrary Lagrangian–Eulerian method.
large deformation The deformations are so large so the nonlinear effect of the change
in geometry or stress stiffening need to be accounted for.
lateral aspect ratio The ratio of the length of a structure in the plane of a wafer to its
width in that plane. See also aspect ratio and vertical aspect ratio.
plane strain An assumption on the strain field where all out-of-plane strain
components are assumed to be zero.
plane stress An assumption on the stress field where all out-of-plane stress
components are assumed to be zero.
principle of virtual work States that the variation in internal strain energy is equal to
the work done by external forces.
proof mass A predetermined test mass in a measurement device or machine that serves
as the reference mass for the quantity to be measured.
residual strain Strain remaining in a structure after some operation, for example, the
strain resulting from cooling a system after high-temperature bonding. It often appears
as the initial strain for any consequent tasks. See also strain and initial strain.
residual stress Stress remaining in a structure after some operation, for example, stress
resulting from cooling a system after high-temperature bonding. It often appears as the
initial stress for any consequent tasks. See also stress and initial stress.
sensor A device that measures a physical variable such as temperature or pressure and
converts it (usually) to an electrical signal.
slide-film damping The damping effect of a fluid between two solid surfaces when the
distance between them is small compared to their area and there is mostly tangential
movement (lubrication). Slide-film damping is a type of thin-film damping.
squeeze-film damping The damping effect of a fluid between two solid surfaces when
the distance between them is small compared to their area and there is mostly normal
movement. Squeeze-film damping is a type of thin-film damping.
stress Internal forces in a material. Normal stresses are defined as forces/area normal
to a plane, and shear stresses are defined as forces/area in the plane. A fundamental
concept in structural mechanics.
thermoelastic damping When an elastic rod is stretched reversibly and adiabatically its
temperature drops. The drop in temperature compensates for the increase in entropy
caused by the stress in the rod (since the process is reversible the entropy remains
constant). In a general deformation, the strain in a structure is nonuniform and so the
temperature changes in a nonuniform manner throughout the structure. The
irreversible flow of heat that accompanies these temperature changes result in energy
loss and damping of the structure.
transducer A device that converts one type of energy to another. It often refers to a
device that responds to a physical parameter and converts it to an electrical signal.
vertical aspect ratio The ratio of the height of a structure perpendicular to a wafer’s
surface to its depth in the wafer’s plane.
S-parameters 52 squeeze-film 61
INDEX| 249
domain nodes theory 155
electric currents 113 electrostatics, boundary elements inter-
electrostatics 75 face 100
dopant concentrations, silicon 239 emailing COMSOL 27
drain node 148 emitter node 147
entropy balance equation 208
E edge nodes
equation of continuity 155
electric currents 113
equivalent circuit 37
electrostatics 75
error message, electrical circuits 63
Einstein summation 236
Euler angles 240
elastoresistance effect 233
exporting
elastoresistance tensor 236
SPICE netlists 68
electric currents interface 110
external current density (node) 121
theory 168
external I vs. U (node) 150
electric displacement field (node) 84
external I-terminal (node) 152
electric fields theory 155
external surface charge accumulation
electric insulation (node) 122
(node) 83
electric point dipole (node)
external U vs. I (node) 151
electric currents 135
electric point dipole (on axis) (node) 136 F FBAR filters 34
electric potential (node) 82, 106 fixed current, ports 50
electric scalar-scalar potential coupling floating potential (node)
(node) 108 electric currents 122
electric shielding (node) 130 electrostatics 89, 106
electrical circuit interface 138 fluid flow 60
theory 171 fluid-structure interaction 60, 195
electrical circuits force calculation (node) 80
modeling techniques 63 forced voltage, port 50
electromechanical coupling factor 40 forces, continuum mechanics and 189
electromechanical forces (node), multi- frequency domain study 169
physics 186 force calculation and 80
electromechanics interface 185, 195 theory, electric currents 157
theory 189
G gate node 148
electroquasistatic criterion 47
Gauss’ law and charge relaxation theory
electrostatic forces 32
155
electrostatic point dipole (node) 99
Gauss’ law equation 158
electrostatically actuated structures 185,
Gibbs free energy 209
195
ground (node) 82, 106
electrostatics interface 72
ground node (node) 140
250 | I N D E X
H heat transfer, in MEMS devices 56 multiphysics
line current source (on axis) (node) 134 oscillatory complex numbers 214
INDEX| 251
face 218 S-parameter calculations
theory 233 electric field, and 51
piezoresistivity, shell interface 226 SPICE
theory 233 exporting 68
Piola-Kirchhoff stress 208 SPICE netlists 67, 153
PNP BJT (node) 147 squeeze-film damping 61
point charge (node) 93 standard notation, structural mechanics
point charge (on axis) (node) 94 interfaces 236
point current source (node) 134 standard settings 21
point current source (on axis) (node) study types
135 electric currents 169
point nodes subcircuit definition (node) 146
electric currents 113 subcircuit instance (node) 147
electrostatics 75 surface charge density (node) 83, 107
predefined couplings, electrical circuits switch (node) 145
64 symmetry plane (node) 84
pyroelectricity (node) 199 symmetry properties, piezoresistivity in-
pyroelectricity multiphysics interface 198 terfaces 236
252 | I N D E X
thermoelasticity interface 204
theory 207
thin low permittivity gap (node) 86
thin-film flow, definition 60
time dependent study 157, 169
TMAC 62
transformer (node) 149
V variables
S-parameters 52
Voigt notation 236
voltage source (node) 141
voltage-controlled current source
(node) 144
voltage-controlled voltage source (node)
143
voltmeter (node) 140
INDEX| 253
254 | I N D E X