Introduction To A CDC Module
Introduction To A CDC Module
AC/DC Module
Introduction to the AC/DC Module
© 1998–2018 COMSOL
Protected by patents listed on www.comsol.com/patents, and U.S. Patents 7,519,518; 7,596,474; 7,623,991; 8,457,932;
8,954,302; 9,098,106; 9,146,652; 9,323,503; 9,372,673; and 9,454,625. Patents pending.
This Documentation and the Programs described herein are furnished under the COMSOL Software License
Agreement (www.comsol.com/comsol-license-agreement) and may be used or copied only under the terms of the
license agreement.
COMSOL, the COMSOL logo, COMSOL Multiphysics, COMSOL Desktop, COMSOL Server, and LiveLink are either
registered trademarks or trademarks of COMSOL AB. All other trademarks are the property of their respective
owners, and COMSOL AB and its subsidiaries and products are not affiliated with, endorsed by, sponsored by, or
supported by those trademark owners. For a list of such trademark owners, see www.comsol.com/trademarks.
Version: COMSOL 5.4
Contact Information
Visit the Contact COMSOL page at www.comsol.com/contact to submit general inquiries, contact
Technical Support, or search for an address and phone number. You can also visit the Worldwide
Sales Offices page at www.comsol.com/contact/offices for address and contact information.
If you need to contact Support, an online request form is located at the COMSOL Access page at
www.comsol.com/support/case. Other useful links include:
Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
The Use of the AC/DC Module . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
The AC/DC Module Physics Interfaces. . . . . . . . . . . . . . . . . . . . 10
Physics Interface Guide by Space Dimension and Study Type . . . 14
Tutorial Example: Modeling a 3D Inductor. . . . . . . . . . . . . . . . . 17
|3
4|
Introduction
|5
benchmark models for verification and validation of the AC/DC interfaces. See
Tutorial Example: Modeling a 3D Inductor to start exploring now.
This guide will give you a jump start to your modeling work. It has examples of
the typical use of the module, a list of all the AC/DC interfaces including a short
description for each, and a tutorial example that introduces the modeling
workflow.
6|
Figure 1: Electric potential distribution and magnetic flux lines for a MEMS inductor.
In Figure 2, Figure 3, and Figure 4, the dynamics of a magnetic brake is shown. The
brake consists of a disc of conductive material and a permanent magnet. The
magnet generates a constant magnetic field, in which the disc is rotating. When a
conductor moves in a magnetic field it induces currents, and the Lorentz forces
from the currents counteracts the spinning of the disc.
The braking torque is computed from the magnetic flux and eddy current
distributions. It is then fed into a rigid body model for the dynamics of the
spinning disc that is solved simultaneously with the electromagnetic model.
|7
Figure 2: The eddy current magnitude and direction in the spinning disc of a magnetic brake.
8|
Figure 4: The time evolution of the dissipated power.
The AC/DC Module includes a range of tools to evaluate and export results, for
example, the evaluation of force, torque, and lumped circuit parameters like
resistance, capacitance, inductance, impedance, and scattering matrices
(S-parameters). Lumped parameters can also be exported in the Touchstone file
format.
The combination of fully distributed electromagnetic field modeling and
simplified circuit based modeling is an ideal basis for design, exploration, and
optimization. More complex system models can be exploited using circuit-based
modeling while maintaining links to full field models for key devices in the circuit
allowing for design innovation and optimization on both levels.
|9
The AC/DC Module Physics Interfaces
The AC/DC interfaces are based upon Maxwell’s equations or subsets and special
cases of these together with material laws. Different subsets of Maxwell’s
equations in combination with special material laws result in different electric,
magnetic, or electromagnetic formulations. In the module, these laws of physics
are translated by the AC/DC interfaces to sets of partial differential equations with
corresponding initial and boundary conditions.
An AC/DC interface defines a number of features. Each feature represents a term
or condition in the underlying Maxwell-based formulation and may be defined in
a geometric entity of the model, such as a domain, boundary, edge (for 3D
components), or point.
Figure 5 on the next page shows the Model Builder window and the Settings
window for the selected Charge Conservation 1 feature node in the Electrostatics
interface. The Charge Conservation 1 node adds the terms representing
Electrostatics to the model equations for a selected geometrical domain in the
model.
Furthermore, the Charge Conservation 1 node may link to the Materials node to
obtain physical properties, in this case the relative permittivity of a user-defined
dielectric. The properties, defined by the Dielectric material node, may be
functions of the modeled physical quantities, such as temperature. The Zero
10 |
Charge 1 boundary condition feature adds the natural boundary conditions that
limit the Electrostatics domain.
Figure 5: The Model Builder window (to the left), and the Settings window for Charge Conservation 1 for
the selected feature node (to the right). The Equation section in the Settings window shows the model
equations and the terms added to the model equations by the Charge Conservation 1 node. The added
terms are underlined with a dotted line. The text also explains the link between the material properties
represented by the Dielectric node and the values for the Relative permittivity.
The AC/DC Module has several AC/DC interfaces ( ) for different types of
electric and magnetic modeling. Figure 6 shows the AC/DC interfaces as well as
those under the Heat Transfer branch, which contains the Joule Heating and
| 11
Induction Heating interfaces. Also see Physics Interface Guide by Space Dimension
and Study Type. A brief overview of the AC/DC interfaces follows.
Figure 6: The AC/DC Module physics interfaces as displayed in the Model Wizard.
E LECTRIC C URRENTS
The Electric Currents interface ( ) is used to model DC, AC, and transient
electric current flow in conductive and capacitive media. This physics interface
solves a current conservation equation for the electric potential. Some examples of
its use are designing busbars for DC current distribution and designing AC
capacitors.
12 |
E LECTRICAL C IRCUIT
The Electrical Circuit interface ( ) has the equations to model electrical circuits
with or without connections to a distributed fields model, solving for the voltages,
currents, and charges associated with the circuit elements. Circuit models can
contain passive elements like resistors, capacitors, and inductors as well as active
elements such as diodes and transistors.
E LECTROSTATICS
The Electrostatics interface ( ) solves a charge conservation equation for the
electric potential given the spatial distribution of the electric charge. It is used
primarily to model charge conservation in dielectrics under static conditions. This
physics interface applies also under transient conditions but then it is usually
combined with a separate transport model for single species or multispecies charge
transport. Such transport models can be found in the Chemical Reaction
Engineering Module and in the Plasma Module. AC simulations are supported
and used especially together with the Piezoelectric Devices interface in the
Acoustics Module, Structural Mechanics Module, and MEMS Module. Some
typical applications that are simulated in the Electrostatics interface are capacitors,
dielectric sensors, and bushings for high voltage DC insulation.
M AGNETIC F IELDS
The Magnetic Fields interface ( ) solves Ampère’s law for the magnetic vector
potential. It is used to model magnetostatics, AC, and transient
| 13
magnetodynamics. Magnets, magnetic actuators, electric motors, transformers,
induction based nondestructive testing, and eddy current generation are typical
applications for this physics interface. It supports both linear media, media with
magnetic saturation and time-domain hysteresis.
14 |
I NDUCTION H EATING
The Induction Heating interface ( ) combines all features from the Magnetic
Fields interface in the time-harmonic formulation with the Heat Transfer interface
to model induction and eddy current heating. The predefined multiphysics
couplings add the electromagnetic power dissipation as a heat source, and the
electromagnetic material properties can depend on the temperature. This physics
interface is based on the assumption that the magnetic cycle time is short
compared to the thermal time scale (adiabatic assumption).
J OULE H EATING
The Joule Heating interface ( ) combines all features from the Electric Currents
interface with the Heat Transfer interface to model resistive heating and heating
due to dielectric losses. The predefined multiphysics couplings add the
electromagnetic power dissipation as a heat source, and the electromagnetic
material properties can depend on the temperature.
AC/DC
| 15
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
Electromagnetic Heating
16 |
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
| 17
Tutorial Example: Modeling a 3D Inductor
Inductors are used in many applications for low-pass filtering or for impedance
matching of predominantly capacitive loads. They are used in a wide frequency
range from near static up to several MHz. An inductor usually has a magnetic core
to increase the inductance while keeping its size small. The magnetic core also
reduces the electromagnetic interference with other devices as the magnetic flux
tends to stay within it. Because there are only crude analytical or empirical
formulas available to calculate impedances, computer simulations or
measurements are necessary during the design stage. In general, inductor
modeling is more complex than modeling resistors and capacitors but similar
principles apply. Using an external CAD software to design and draw the model,
the geometry is imported into the AC/DC Module for static and frequency
domain analysis. The inductor geometry is shown in Figure 7.
Core
Copper
winding
96 mm
Feed gap
18 |
frequency, capacitive effects become significant. The equivalent circuit model
involves connecting an ideal capacitor in parallel with the DC circuit. The circuit
parameters are obtained by analyzing the frequency dependent impedance
obtained from a frequency domain simulation. In this tutorial, the AC analysis is
done up to the point when the frequency dependent impedance is computed.
These step-by-step instructions guide you through the detailed modeling of an
inductor in 3D. The module also has a physics interface to model electrical
circuits, which is detailed in the magnetostatic part of this model. The first step is
to perform a DC simulation.
Model Wizard
Note: These instructions are for the user interface on Windows but apply, with
minor differences, also to Linux and Mac.
1 To start the software, double-click the COMSOL icon on the desktop. When
the software opens, you can choose to use the Model Wizard to create a new
COMSOL model or Blank Model to create one manually. For this tutorial, click
the Model Wizard button.
If COMSOL is already open, you can start the Model Wizard by selecting
New from the File menu and then click Model Wizard .
The Model Wizard guides you through the first steps of setting up a model. The
next window lets you select the dimension of the modeling space.
2 In the Space Dimension window click 3D .
3 In the Select Physics tree under AC/DC, select Magnetic Fields (mf) .
4 Click Add then click Study .
5 In the Select Study window under General
Studies, click Stationary .
6 Click Done .
G e o me t r y 1
The main geometry is imported from a file. Air domains are typically not part of a
CAD geometry, so these are added to the model. For convenience, three
additional domains are defined in the CAD file and then used to define a narrow
feed gap where an excitation is applied.
| 19
Note: The location of the files used in this exercise varies based on the
installation. For example, if the installation is on your hard drive, the file path
might be similar to:
C:\Program Files\COMSOL\COMSOL53a\Multiphysics\applications\.
Import 1
1 On the Home toolbar, click Import .
Note: On Linux and Mac, the Home
toolbar refers to the specific set of controls
near the top of the Desktop.
2 In the Settings window under Import,
click Browse. Then navigate to your
COMSOL Multiphysics installation
folder, locate the subfolder
ACDC_Module\Inductive_Devices_and
_Coils, select the file
inductor_3d.mphbin and click Open.
3 Click Import.
Sphere 1
1 On the Geometry toolbar, click
Sphere .
2 Go to the Settings window for Sphere.
Under Size in the Radius text field, enter
0.2.
20 |
Form Union
1 In the Model Builder under Geometry 1,
click Form Union .
2 In the Settings window for Form
Union/Assembly, click Build All .
3 On the Graphics window toolbar, click the
Zoom Extents button and then the
Wireframe Rendering button.
The geometry should match this figure.
Definitions - Selections
Use the Selection feature available under the Definitions node to make the set up
of materials and physics interfaces easier. Start by defining the domain group for
the inductor winding and continue by adding other useful selections. These steps
illustrate how to set up the Explicit nodes and rename the geometric selections
accordingly.
1 On the Definitions toolbar, click Explicit .
| 21
2 Repeat Step 1 and add a total of six (6) Explicit nodes . Or right-click the
first Explicit node and select Duplicate .
The following steps are done one at a time for each node using the table below.
3 In the Model Builder click an Explicit node to open the Settings window.
Note: There are many ways to select geometric entities. When you know the
domain to add, such as in this exercise, you can click the Paste Selection
button and enter the information in the Selection text field. In this example
for the Explicit 1 node, enter 7,8,14 in the Paste Selection window. For more
information about selecting geometric entities in the Graphics window, see the
COMSOL Multiphysics Reference Manual.
4 In the Settings window under Input Entities, select Domain from the
Geometric entity level list. Use the table as a guide, and in the Graphics window,
select the domains as indicated. Then rename the node by changing its Label in
the Settings window.
DEFAULT NODE SELECT THESE DOMAINS NEW NAME FOR THE NODE
NAME
22 |
5 After renaming all the nodes under Definitions, the sequence of nodes should
match this figure.
Use the Infinite Element feature available under the Definitions node to emulate
an infinite open space surrounding the inductor. The infinite element coordinate
scaling is recognized by any physics interface supporting it - other physics
interfaces are automatically disabled in the infinite element domains.
1 On the Definitions toolbar click Infinite Element Domain .
2 Go to the Settings window for Infinite Elements. Under Domain Selection
select Infinite Elements from the Selection list.
3 Under Geometry from the Type list, select Spherical.
Mate rials
Now define the material settings for the different parts of the inductor. Use the
Selections defined in the previous section.
Copper
1 On the Home toolbar click Add Material .
| 23
2 Go to the Add Material window.
In the Materials tree under
AC/DC, right-click Copper
and choose Add to
Component 1 from the menu.
3 In the Model Builder, click
Copper .
Air
1 Go to the Add Material window.
In the Materials tree under
Built-In, right-click Air and
choose Add to Component 1
from the menu.
2 Close the Add Material window.
3 In the Model Builder click Air
.
4 Go to the Settings window for
Material. Under Geometric Entity
Selection select Non-conducting
from the Selection list.
24 |
3 Under Geometric Entity Selection select Core from
the Selection list.
4 Under Material Contents in the table, enter these settings in the Value column:
- 1e3 for the Relative permeability
- 0 for the Electrical conductivity
- 1 for the Relative permittivity
Ma gn e t ic F ie l d s ( mf )
The model is solved everywhere except in the feed gap region. By keeping this
void, its default boundary conditions support surface currents closing the current
loop. In inductor modeling it is crucial to respect current conservation as this is a
law of nature (Ampère’s law).
| 25
1 In the Model Builder click Magnetic Fields .
2 In the Settings window for Magnetic Fields, select
domains 1–8 and 10–14 only (all domains except
9). Or first click the Clear Selection button ,
then click the Paste Selection button and enter
1–8, 10–14 in the Selection text field.
Coil 1
A dedicated feature is added for the solid copper winding of the inductor. This
feature makes it easier to excite the model in various ways.
1 On the Physics toolbar click Domains and choose Coil .
2 In the Settings window for Coil, under Domain Selection, choose Winding
from the Selection list.
3 In the Model Builder, expand the Coil 1 node.
4 Right-click the Geometry Analysis 1 subnode and choose Output . There
is already a default Input node.
5 Click the Input 1 node.
6 In the Settings window for Input 1, click the Paste Selection button and
enter 58 in the Selection text field.
7 Click OK.
8 Click the Output 1 node.
9 In the Settings window for Output 1, click the Paste Selection button and
enter 79 in the Selection text field.
10Click OK.
The Coil with default settings will treat it as a solid conductor and apply a current
of 1 A flowing from the input to the output.
The node sequence under Magnetic Fields should at this point match this figure.
26 |
Me sh 1
The infinite elements region requires a swept mesh to maintain a good effective
element quality despite the steep radial scaling. To enable this meshing method,
perform the following steps.
1 In the Model Builder click Magnetic Fields .
2 In the Settings window for Magnetic Fields, locate the Physics-Controlled
Mesh section.
Study 1
Before the magnetostatic model is ready to solve an additional study step needs to
be added for computing the exciting coil current.
1 Right-click the Study 1 node and select Study steps>Other>Coil Geometry
Analysis .
2 Right-click Coil Geometry Analysis and select Move Up . The Study
node should at this stage match this figure.
| 27
Results
By adding Selections and adjusting the Data Sets under the Results node you can
create more interesting plots.
1 In the Model Builder, expand the Results>Data Sets node. Right-click
Study 1/Solution 1 and choose Duplicate .
2 Right-click Study 1/Solution 1 (3) and choose Selection.
28 |
3 Go to the Settings window for Selection. Under Geometric entity selection:
- From the Geometric entity level list, select Domain.
- From the Selection list, select Winding.
3D Plot Group 2
1 On the Results toolbar, click 3D Plot
Group .
An additional toolbar containing Plot
Tools for the 3D Plot Group appears when
the 3D Plot Group is selected in the
Model Builder.
| 29
2 On the 3D Plot Group 2 toolbar, click Streamline .
30 |
7 Go to the Settings window for Volume.
Under Data from the Data set list, select
Study 1/Solution 1 (4).
8 Click the Plot button.
On the Graphics window toolbar, click the
Zoom In button twice.
| 31
Derived Values
Next, evaluate the coil inductance and resistance and compare them respectively
to the inductance estimated from the total magnetic energy and to the resistance
defined using Ohm’s law.
1 On the Results toolbar click Global Evaluation .
2 In the Settings window for Global Evaluation, locate the Expression section.
3 In the table, enter the following settings:
Ma gn e ti c F ie l d s ( mf )
As shown in the Magnetic Flux Density plot, this is effectively contained within
the core. Thus, there should be little need to use infinite elements in this model.
Solve the model without the infinite elements to test this hypothesis.
32 |
1 In the Model Builder click Magnetic Fields .
2 In the Settings window for Magnetic Fields, select domains 5–8 and 14 only.
Study 1
Derived Values
2 In the Settings window for Global Evaluation, click the Evaluate button.
The Resistance and Inductance values display in the Table window. The results
change very little.
| 33
Model Wizard
Ma gn e ti c F ie l d s ( MF )
Change the field excitation to connect to the circuit part of the model.
1 In the Model Builder, under
Component 1 (comp1)>Magnetic Fields
(mf), click Coil 1.
2 In the Settings window for Coil ,
change the Coil excitation to
Circuit (current), as in the figure.
34 |
Electrical Circuit (cir)
Voltage Source 1
1 In the Model Builder select the Electrical
Circuit (cir) node.
2 On the Electrical Circuit toolbar click
Voltage Source .
Go to the Settings window for Voltage
Source. The default settings correspond
to a DC source of 1 V. Under Node
Connections enter 1 and 0 in the Node
names table as in the figure to the right.
| 35
Resistor 1
1 On the Electrical Circuit toolbar click
Resistor .
2 Go to the Settings window for Resistor.
Under Node Connections enter 1 and 2
in the Node names table as in the figure to
the right.
3 Under Device Parameters in the R text
field, enter 100[mohm].
External I Vs. U 1
There is a special feature for connecting the
circuit to the finite elements model.
1 On the Electrical Circuit toolbar click
External I Vs. U .
2 In the Settings window for External I Vs.
U under Node Connections, enter Node
names as in to the figure to the right.
3 Under External Device, from the Electric
potential V list choose Coil voltage
(mf/coil1).
The sequence of nodes under Electrical
Circuit should match this figure.
Study 1
36 |
Derived Values
| 37
This can also be seen in the 3D Plot Group.
38 |
1 On the Home toolbar click Add
Study .
2 In the Add Study window click
Frequency Domain .
3 In the Physics interfaces in study
table, under Solve, click to clear the
check mark , which removes the
Electrical Circuit (cir) interface.
4 Click Add Study.
5 Close the Add Study window.
A Study 2 node with a Frequency
Domain study step is added to the
Model Builder.
Definitions
At high frequencies, the skin depth of the copper conductor in the winding
becomes difficult to resolve. The solution is to exclude the interior of the copper
domain from the simulation and instead represent the winding by a boundary
condition, which introduces surface losses via an equivalent surface impedance.
For this purpose add a boundary selection and an associated surface material.
Conductor Boundaries
1 On the Definitions toolbar click Explicit .
| 39
2 Select Domains 7, 8, and 14 only.
3 Go to the Settings window for Explicit.
Under Output Entities from the
Output entities list, select Adjacent
boundaries.
4 In the Label text field enter Conductor
Boundaries.
Mate rials
Copper 1
1 On the Home toolbar click Add Material .
2 Go to the Add Material window.
In the Materials tree under
AC/DC, right-click Copper and
choose Add to Component 1
from the menu.
3 Close the Add Material window.
4 In the Model Builder, click
Copper 1.
40 |
5 Go to the Settings window for Material. Under Geometric Entity Selection:
- From the Geometric entity level list, select Boundary.
- From the Selection list, select Conductor Boundaries.
Ma gn e t ic F ie l d s ( mf )
| 41
2 Go to the Settings window for Magnetic Fields. Under Domain Selection from
the Selection list, select Non-conducting without IE.
Disable Coil 1
The Coil does not apply to an active domain anymore so it needs to be disabled.
1 In the Model Builder under Magnetic Fields, right-click Coil 1 and choose
Disable . The node is grayed out.
42 |
Lumped Port 1
With the Coil replaced by the Impedance Boundary Condition a different
boundary feature, which is more appropriate for high frequency modeling, must
be used to excite the model.
1 On the Physics toolbar click Boundaries and choose Lumped Port .
A Lumped Port node is added to the sequence in the Model Builder.
Ampère's Law 2
Apart from the surface losses in the copper
conductor, there are also AC losses in the
core. The core loss is introduced by
specifying an effective complex permeability
for the magnetic iron in the core, a quantity
which is often available from the
manufacturer. There is a specific constitutive
relation for this.
| 43
1 On the Physics toolbar click Domains and choose Ampère's Law .
A second Ampère's Law node is added to the Model Builder.
2 Go to the Settings window for Ampère's Law. Under Domain Selection, from
the Selection list, select Core.
3 In the Settings window for Ampère’s
Law, locate the Magnetic Field section.
4 From the Constitutive relation list,
choose Magnetic losses.
5 From the μ′ list, choose User defined. In
the μ′ text field, type 1200.
6 From the μ′′ list, choose User defined.
In the μ′′ text field, type 100.
Study 2
44 |
Results
Plot the surface current distribution at the lowest frequency solved for. This
frequency is below resonance and the device is of an inductive nature.
Data Sets
1 In the Model Builder under Results>Data Sets, right-click
Study 2/Solution 3 (5) and choose Duplicate .
2 Right-click Study 2/Solution 3 (6) and choose Selection.
| 45
3 Go to the Settings window for Selection. Under Geometric Entity Selection
from the Geometric entity level list, choose Boundary.
4 From the Selection list, choose Conductor Boundaries.
3D Plot Group 4
1 On the Results toolbar click 3D Plot Group .
A 3D Plot Group 4 node is added to the sequence.
2 Go to the Settings window for 3D Plot Group. Under Data from the Data set
list, choose Study 2/Solution 3 (6).
46 |
name, as in this exercise, you can also enter mf.normJs directly in the
Expression text field).
7 Click the Plot button.
The plot shows the surface current distribution and should match this figure.
Change to the highest frequency solved for and compare the results. Finish the
modeling session by plotting the real and imaginary parts of the coil impedance.
1D Plot Group 5
1 On the Results toolbar click 1D Plot Group .
2 In the Settings window for 1D Plot Group under Data, select
Study 2/Solution 3 (5) from the Data set list.
3 On the 1D Plot Group 5 toolbar click Global .
| 47
4 Go to the Settings window for Global. Under y-Axis Data in the Expression
column, enter real(mf.Zport_1). In the Description column, enter Real
part of impedance.
Note: In general, you can click the Replace Expression to see what variables
are available in the predefined lists. In this case, you can enter the information
directly into the table.
5 Click the Plot button.
The resistive part of the coil impedance peaks when passing through the resonance frequency.
48 |
6 Go to the Settings window for Global and modify it according to the figure
below.
| 49
The reactive part of the coil impedance changes sign when passing through the resonance frequency,
going from inductive to capacitive.
50 |