AC/DC Module Updates

For users of the AC/DC Module, COMSOL Multiphysics® version 5.6 offers a new Magnetic Fields, Currents Only interface, empirical loss models, and various improvements to the existing Electrical Circuit interface. Learn about these and other AC/DC Module updates in more detail below.

New Physics Interface for Inductance Matrix Computations

The new Magnetic Fields, Currents Only interface is designed to efficiently compute the lumped inductance matrix of complex circuits in 3D, such as those commonly found on printed circuit boards (PCBs). It can be used to calculate the partial contributions from magnetic fields generated by open conductors, also known as nonsolenoidal conductors. The ability to handle open conductors significantly reduces modeling complexity for users needing to solve for mutual inductance matrices.

The interface uses the magnetic vector potential as the dependent variable, and computes the magnetic fields generated by currents under the assumption that all regions are nonmagnetic (in other words; that they have a uniform relative magnetic permeability of "one"). This new formulation does not require the conductor currents to be divergence-free. In free space, it returns the value of the Biot–Savart integral. The interface supports the Stationary and the Stationary Source Sweep study steps. The latter is an efficient way of sweeping over many terminals. You can see this interface in use in the new Inductance Matrix Calculation of PCB Coils tutorial model and in the updated Magnetic Field of a Helmholtz Coil model.

An array of copper PCB coils showing the magnetic flux density norm in a dark blue to white color gradient and magnetic fields as arrows. PCB coils An inductively coupled array of PCB coils, subjected to a coil sweep using the Magnetic Fields, Currents Only interface. The resulting L-matrix can be used to construct a lumped circuit representation of the PCB.

Loss Models for Laminated Cores and Yokes in Motors and Transformers

The new Loss Calculation feature uses empirical models (Steinmetz, Bertotti, or User defined) to calculate iron loss caused by magnetic hysteresis, eddy currents, and other effects in (laminated) cores or yokes in, for example, electric motors and transformers. As a fallback option for materials like copper, it is able to provide the cycle average resistive losses.

In the Model Builder, Loss Calculation is available as a subnode to the Ampère's Law node, domain Coil node, Faraday’s Law node, and Magnetic Flux Conservation node. It is available for both Time Dependent and Frequency Domain studies. For Time Dependent studies, the Loss Calculation feature is designed to be used in combination with the dedicated Time to Frequency Losses study step. This study step determines the different harmonics present in the time-dependent solution, and inserts them into one of the frequency-dependent empirical loss models. You can see this new feature in these updated models:

The COMSOL Multiphysics version 5.6 UI showing the Loss Calculation settings with the Steinmetz loss model selected and a PM motor model in the Graphics window. Demonstrating the Steinmetz loss model The empirical loss model Steinmetz, used to determine the cycle-average (frequency-dependent) losses in the laminated iron of a stator yoke. The different harmonics can be analyzed individually.

Wide Support for Eigenfrequency Analysis

The Eigenfrequency study is now supported for most of the AC/DC Module interfaces: Electric Currents, Electric Currents in Shells, Electric Currents in Layered Shells, Electrical Circuit, Electrostatics, and Magnetic Fields. In addition to supporting full-wave cavity mode analysis in the Magnetic Fields interface, it is possible to run eigenfrequency analyses with models involving electrical circuits. The eigenfrequency support is primarily developed for the AC/DC Module, but other modules that provide one of the affected physics interfaces will benefit from it too.

The COMSOL Multiphysics version 5.6 UI showing the Model Builder, Global Evaluation settings with the Data and Expressions sections expanded, and a probe plot for the eigenfrequency analysis of an RLC circuit. Eigenfrequency analysis example The resonance peak of a simple RLC circuit. The eigenfrequency and Q-factor are analyzed and compared to analytically determined values.

New and Enhanced Functionality for the Electrical Circuit Interface

For Time Dependent studies, the Electrical Circuit interface has been equipped with an "event-based" Switch feature. This allows you to model the "instantaneous" on-off switching of certain connections in the circuit. The switch can be current controlled, voltage controlled, or controlled by user-defined Boolean expressions.

Furthermore, Parameterized Subcircuit Definitions are added. Together with the Subcircuit Instance, these allow you to create your own building blocks containing smaller circuits, and use multiple parameterized variants of those in your larger circuit. Finally, the state, event, and solver machinery has been improved, especially the transient modeling of nonlinear (semiconductor) devices, which has become more robust.

The circuit improvements are primarily developed for the AC/DC Module, but other modules that provide access to the Electrical Circuit interface will benefit too. You can view the new functionality in these updated models:

New Default Plots for the Magnetic Physics Interfaces

New default plots have been developed specifically for physics interfaces that support magnetic fields: Magnetic Fields, Magnetic Fields, No Currents, and the Rotating Machinery, Magnetic multiphysics interface. The magnetic flux density norm is plotted as a multislice or a surface plot, together with streamlines or contours to indicate the field direction. In 2D and 2D axisymmetric models, the streamlines show a divergence-free magnetic field. Whenever possible, a contour plot of the out-of-plane magnetic vector potential is included, which is more exact than the streamlines, in the sense that the contours are closed, and the contour density is directly proportional to the local magnetic flux density. In 3D, the streamlines follow the magnetic field after it has been projected onto the slice planes. The resulting pattern is not necessarily divergence free, but still provides a clear, intuitive indication of what the 3D field will look like. You can view this new plot functionality in the following models:

A topology optimized magnetic circuit model showing the magnetic flux density norm in a Volume plot in white to dark blue color gradient. Demonstrating the new default plot for loudspeakers The new default plots, as seen for loudspeakers in the Topology Optimization of a Magnetic Circuit model.
A 2D generator model showing the magnetic flux density norm in a Surface plot in white to dark blue color gradient. Demonstrating the new default plot for motors and generators The new default plots, as seen for motors and generators in the Generator in 2D model.

Ferroelectric Material Model in Electrostatics

For the Charge Conservation feature in the Electrostatics interface, the list of dielectric material models has been extended with a Ferroelectric material option. The ferroelectric material model is the dielectric equivalent of the Jiles Atherton magnetic hysteresis model; instead of a time-varying magnetization, it models a time-varying polarization. You can view this feature in the new Hysteresis in Ferroelectric Material tutorial model.

This material model is the electrostatic part of a new Ferroelectroelasticity multiphysics interface, intended for analysis of ferroelectric materials, which exhibit nonlinear piezoelectric properties. The Ferroelectric material option is included in the AC/DC Module. However, the Ferroelectroelasticity multiphysics interface additionally requires either the MEMS Module or the Structural Mechanics Module.

The COMSOL Multiphysics version 5.6 UI showing the Charge Conservation feature settings with the Constitutive Relation D-E and Ferroelectric Material Properties sections expanded and a hysteresis curve to the right. Demonstrating the ferroelectric material model A screenshot of the Hysteresis in Ferroelectric Material tutorial model, showing both the used material model and the resulting hysteresis curve.

Magnetic Materials from Bomatec

The AC/DC material library has been extended with magnetic materials from Bomatec. This includes all the NdFeB standard grades ("normal", M, H, SH, UH, EH, AH, BH, /S, M/S, H/S, SH/S, UH/S, EH/S, AH/S, BH/S, H/ST, SH/ST, UH/ST, EH/ST, AH/ST, and BH/ST); NdFeB (bonded, injection molded, and extruded); ferrite (isotropic, anisotropic, and injection molded); SmCo (injection molded); SmFeN (injection molded); SmCo5 grades; Sm2Co17 grades; and AlNiCo (cast and sintered).

The materials include electromagnetic properties — relative permittivity , recoil permeability , and remanent magnetic flux density norm — as well as several properties related to thermal modeling: thermal conductivity, ; density, ; and heat capacity at constant pressure .

Enhanced Feature: Coil with Slanted Cut

The Coil domain feature (both for the Single Conductor model and the Homogenized Multiturn case) has been equipped with support for a Slanted Cut. The Slanted Cut setting allows you to relax the constraints used for determining the coil geometry (it applies the input/output boundary constraint "on average" rather than locally). Relaxing the constraint is needed for cases where the natural direction of the current (or the wires) is at an angle with respect to the input or output normal. A typical example is a helical coil, where the input/output boundaries represent a periodicity plane — as is the case for twisted three-phase cables, or fully resolved Litz wires or Milliken conductors. You can see this new feature in the updated Cable Tutorial Series tutorial model.

The COMSOL Multiphysics version 5.6 UI showing the Input settings with the slanted cut option selected for a submarine cable model whose results appear in the Graphics window. Demonstrating the slanted cut option in a cable model The slanted cut setting, demonstrated in a three-phase cable model.
The COMSOL Multiphysics version 5.6 UI showing the Input settings with the slanted cut option selected for a spiral wire model whose results appear in the Graphics window. Demonstrating the slanted cut option in a wire model The effect of the slanted cut setting on the current direction is illustrated in the wire to the right.

Faster Cable Analysis: New Tutorial Model

The Cable Tutorial series has been updated to demonstrate short-twisted periodicity, a kind of periodicity that allows you to reduce the computational effort required for fully resolved 3D twisted cable models, from hours on large cluster systems to about a minute on a laptop. The logic behind the used periodicity is explained in detail, including cases that involve more than two lay lengths, cables with fully resolved Milliken conductors, stranded screens, or a double armor, for example.

The theory for twisted strands and the use of appropriate meshing follows the latest developments in numerical research, and applies to both large-scale insulated submarine cables and umbilical cables, EV charging cables, fully-resolved Litz wires, and the bare cables used in overhead power lines. Additionally, the series discusses the resistive, capacitive, inductive, and thermal properties of cable systems, as well as different bonding schemes: single-point bonding, solid bonding, and cross bonding. You can see this new feature in the updated Cable Tutorial Series tutorial model.

The COMSOL Multiphysics version 5.6 UI showing the 3D Plot Group settings with the Data, Plot Settings, and Color Legend sections expanded, and a double-armor submarine cable model in the Graphics window. Submarine cable model with a double armor A demonstration model showing a double armor, based directly on the design and theory discussed in the Cable Tutorial Series. (Note: The double armor model is not included in the series explicitly, rather, it is directly derived from it.)

New Tutorial Models

COMSOL Multiphysics® version 5.6 brings several new tutorial models to the AC/DC Module.

Hysteresis in Ferroelectric Material

A point graph showing hysteresis loops for different Vmax: blue denotes 600 V, green 800 V, red 1000 V, and turquoise 1600 V. Hysteresis loops model A new model demonstrates how to analyze the nonlinear polarization behavior of ferroelectric materials. The image shows the hysteresis loops fully established after three cycles for different values of the maximum applied voltage.

Application Library Title:

ferroelectric_hysteresis

Download from the Application Gallery

High-Voltage Insulator

A rainbow model of the electric potential of a high-voltage insulator. High-voltage insulator model The image shows the distribution of the electric potential of a typical high-voltage insulator. The insulator is widely used in high-voltage power transmission lines to support the weight of suspended conductors without allowing the current to flow through the tower to the ground.

Application Library Title:

high_voltage_insulator

Download from the Application Gallery

Position Optimization of Grading Rings

A line graph showing tangential electric field without gradient rings, blue; with grading rings, green; and with grading rings at optimized positions, red. Grading rings position optimization model A new model that extends the High-Voltage Insulator model and shows how to find the optimum position of grading rings using the Optimization interface. The figure compares the z-component of the tangential electric field along the surface of the first six sheds from the line end in different cases.

Application Library Title:

grading_ring_optimization

Download from the Application Gallery

Inductance Matrix Calculation of PCB Coils

Inductance matrix calculation in color gradient from dark blue to white. Inductance matrix calculation of PCB coils model A new model shows how to compute the inductance matrix of a series of conductors with the new Magnetic Fields, Currents Only interface. The image shows the computed inductance matrix of a PCB composed of 12 independent conductors.

Application Library Title:

inductance_matrix_pcb_coils

Download from the Application Gallery