Skip to content

Planar Inductor

A planar inductor whose two-layer PCB winding is imported from CAD into an E core: an inductance sweep at 5 and 50 kHz with the copper traces meshed as solids, so skin and proximity effects are part of the result.

  • Component: E core in Ferroxcube 3C90, no gap; two copper layers on a 1.6 mm FR4 board.
  • Simulation: Inductance, frequency sweep at 5 and 50 kHz with 1 V RMS on the winding (Without Circuits); full model.
  • Winding model: Massive (solid).
  • You get: inductance, AC resistance, Q factor, capacitance, core and copper loss per frequency.
  • Files: the winding, coil1.iges, in the CAD folder of the example project.
  • In TRAFOLO: Welcome screen › Examples › Planar Inductor is this model with its results.
  • Run time: about 4 min on 8 CPU cores.

Model description

The winding is the copper of a planar inductor design: two layers on the board, joined by a via. Two E core halves enclose the board around the winding.

Planar inductor with an E ferrite core and a two-layer copper winding on a PCB
The planar inductor: E core and two-layer copper winding on the PCB.

Modelling instructions

Setup

  1. Under DOMAIN TYPE, select Full: the winding is not symmetric about z = 0.
  2. Under ELECTROMAGNETICS, select Inductance.
  3. Tick Without Circuits. TRAFOLO then applies the voltage directly to the winding terminals, and HEAT TRANSFER stays Off.
Setup tab with Full, Inductance and Without Circuits marked
Full (1), Inductance (2), Without Circuits (3).

Materials

  1. Under Global Materials, select Ferroxcube 3C90 100C and click Copy.
  2. Copy Copper and Fiberglass Reinforced Epoxy (FR4) the same way. All three appear under Project Materials.
  3. Under Default Materials, set the material of each component:
Core Coil Other Bobbin Gaps
Ferroxcube 3C90 100C Copper Not set Fiberglass Reinforced Epoxy (FR4) Not set
Materials tab with Copy, the three project materials and the default material per component
Copy (4), Project Materials (5), Default Materials (6).

Core

  1. On Core, click New Group. Group1core appears.
  2. Click Configure next to Group1core.
  3. Set Source to Template, Scale to mm and Gap Type to None.
  4. In Geometry Builder, set Type to E, enter the dimensions in mm, keep Mirror ticked and click Apply:

    A B C D E F
    32 5.25 20.5 2.2 24.9 6.4
  5. Under Translation and Rotation, set x to −0.7685, y to 0.9 and z to 33 mm and click Apply. This places the core over the imported winding.

Group1core settings with Source, Scale, Gap Type None, the E core dimensions and the translation
Group1core: Source, Scale and Gap Type (7), E dimensions and Apply (8), translation and Apply (9).
E core in an angled 3D view with the dimension sketch
The two E core halves with the dimension sketch.

Coil

  1. On Coil, click New Group. Group1coil appears.
  2. Click Configure next to Group1coil.

    Coil tab with New Group and Configure for Group1coil
    Coil tab: New Group (10), Configure (11).

  3. Set Source to CAD and click Import.

    Group1coil settings with Source CAD, Import and Wire Type Massive solid
    Group1coil: Source CAD and Import (12), Wire Type Massive (solid) (16).

  4. In Import CAD geometry, click Import and open coil1.iges from the CAD folder of the Planar Inductor example. Keep Scale mm, set Geometry Handling Algorithm to Unmodified and click Ok.

    Import CAD geometry dialog with the winding file, Scale mm, Unmodified and the two-layer winding
    Import CAD geometry: Import the file (13), Scale and Geometry Handling Algorithm (14), Ok (15).

  5. Set Wire Type to Massive (solid).

Two-layer planar winding in an angled 3D view
The imported two-layer winding in the Coil tab 3D view.

Bobbin

The bobbin is the PCB.

  1. On Bobbin, click New Group. Group1bobbin appears.
  2. Click Configure next to Group1bobbin.
  3. Set Source to Template and Scale to mm.
  4. In Geometry Builder, set Type to Box with X 28, Y 1.6 and Z 60 mm, and click Apply.
  5. Under Translation and Rotation, set x to −0.5, y to 0.9 and z to 30 mm and click Apply.
Group1bobbin settings with the Box dimensions of the board and its translation
Group1bobbin: Source and Scale (17), Box dimensions and Apply (18), translation and Apply (19).

Assembly

  1. Under Air Domain, set Min Z to 0.1 mm, click its lock so that it turns red, and click Apply. The air domain then cuts the winding's two leads, so their ends lie on the boundary and act as terminals.

    Air Domain panel with Min Z 0.1 locked and Apply
    Air Domain: Min Z 0.1 mm, locked (20), then Apply (21).

  2. Click Assemble Geometry.

    Assemble Geometry button
    Assemble Geometry (22).

Assembly table with air, coil, core and board and their materials
Assembly: the bodies and their materials.
Assembled planar inductor in an angled 3D view
The assembled model: core, winding and board.

ElMag

  1. Under Primary Winding, Excitation is Urms [V] because of Without Circuits. Set Coil group connection type to Series.

    Primary Winding with Urms excitation and Series connection
    Primary Winding: Urms [V], set by Without Circuits, and Series (23).

  2. In the coil group table, set Group1coil to Winding Primary and Connection Series.

    Coil group table with Group1coil Primary Series
    Coil group (24): Group1coil, Primary, Series.

There is no Circuits tab: Without Circuits removes it.

Waveform

  1. Under SWEEP BY, select Frequency.
  2. Set Fixed Urms to 1 V.
  3. In Primary Frequency [Hz] List, enter 5000.0 50000.0.
Waveform tab with Frequency sweep, Fixed Urms 1 V and the frequency list
Frequency (25), Fixed Urms 1 V (26), the frequencies (27).

Heat

  1. Under GLOBAL SETTINGS, set Steady-State Temperature to 80 °C. Without a thermal solve it still sets the temperature of the materials, for example the copper conductivity.

    Heat tab global settings with Steady-State Temperature 80 degrees
    Steady-State Temperature 80 °C (28).

Mesh

  1. Under GLOBAL SETTINGS, set Mesh Algorithm to Netgen and Mesh Refinement to Coarse.

    Mesh global settings with Netgen and Coarse
    Mesh Algorithm Netgen and Mesh Refinement Coarse (29).

  2. Click Compute Mesh. TRAFOLO sizes the elements itself, 0.5 mm on the copper; the mesh has about 92,000 elements.

    Compute Mesh button
    Compute Mesh (30).

Mesh of the planar inductor with the mesh information
The mesh: 91,996 elements, finest on the copper.

Solve

  1. Set Parallel cores to the number of performance cores of your CPU.
  2. Set Solver Preset to Custom and Nonlinear Max Iterations to 20.

    Solve tab with Parallel cores and Solver Preset Custom with 20 nonlinear iterations
    Parallel cores (31); Solver Preset Custom with Nonlinear Max Iterations 20 (32).

  3. Click Simulation Setup, then Run Simulation.

    Simulation Setup and Run Simulation buttons
    Simulation Setup, then Run Simulation (33).

Results

At 1 V RMS (Results › Summary and Inductance):

Frequency (kHz) Circuit Inductance (µH) Energy Inductance (µH) Rac Loss (mΩ) Coil loss (mW) Core loss (mW)
5 159.4 159.4 4.85 0.193 1.04
50 156.4 153.7 6.85 0.0028 0.516

The capacitance of the winding is 11.2 pF.

Circuit or Energy Inductance

Circuit Inductance comes from the voltage and current at the winding's terminals, L = Im(V/I)/(2πf): the inductance the circuit sees. Energy Inductance comes from the magnetic energy W stored in the whole model, L = 2W/I². Both describe the same inductance and agree in a converged solution; a difference of more than a few percent points to a convergence or mesh problem. Which one to trust depends on the model:

  • Circuit Inductance with a saturating core or with several windings carrying current: the energy integral then needs a fully converged nonlinear field, and the stored energy belongs to all the currents together.
  • Energy Inductance with one winding, a linear core and the current crowded into a layer thinner than the elements: the energy, integrated over the whole field, is less sensitive to that than the terminal current.

Here the two agree at 5 kHz. At 50 kHz the circuit value is 1.8 % higher. The skin depth in copper is then about 0.3 mm, finer than the 0.5 mm elements on the traces, and the current crowds to the trace edges (loss density picture below), so the terminal current is computed slightly less accurately than the stored energy. With one winding and the core below saturation, Energy Inductance, 153.7 µH, is the better value at 50 kHz; a 1.8 % difference is small for design.

Inductance plot with circuit and energy inductance at 5 and 50 kHz
Circuit and Energy Inductance per frequency.

Losses and Q

At the fixed 1 V, the current falls from 0.20 A at 5 kHz to 0.020 A at 50 kHz as the reactance rises from 5.0 to 49 Ω, so both losses fall. Skin and proximity effects raise the winding's AC resistance (Rac Loss) by 41 %. The core causes 84 % of the loss at 5 kHz and 99 % at 50 kHz: the 3C90 ferrite, not the copper, limits Q.

Results Summary with coil, core and total loss at 5 and 50 kHz and the computation time
Summary: losses at 5 and 50 kHz, 1 V.

Q = 2πf·L/R, with L the Energy Inductance and R = (coil loss + core loss)/I², is 162 at 5 kHz and 38 at 50 kHz. Q Factor in Results shows 223 and 35 because it divides the reactance by Rac Circuit, the real part of V/I: 22 mΩ at 5 kHz and 1.4 Ω at 50 kHz. Rac Circuit holds every loss of the field solution, the core's included, so it is about 5 and 200 times the winding's Rac Loss. The core loss inside the field solution, 0.70 mW at 5 kHz and 0.58 mW at 50 kHz, differs from the Summary's Core Loss, 1.04 and 0.52 mW, which the core-loss model computes after the solve. Use the losses for Q: they are the values the rest of the results use.

Inductance results table at 5 and 50 kHz
Inductance table: L, XL, Rac Circuit, Rac Loss, Q Factor and Capacitance per frequency.
Magnetic flux density in the core at 5 kHz
Flux density at 5 kHz: mostly 0.1 to 0.15 T, up to 0.3 T at the inner corners (Magnetic Flux).
Copper loss density in the two-layer winding at 50 kHz
Copper loss density at 50 kHz (Coil Losses, frame 2).

The example ships with these views but without 3D field data: run it before Open in ParaView (ParaView (post-processing).

Coil and core loss and flux density of the planar inductor in ParaView
Coil and core loss and flux density in ParaView.
Loss density in the planar winding on a logarithmic scale, concentrated along the trace edges
Loss density on a logarithmic scale: skin and proximity effects push the current to the trace edges.