6. Magnetostatics and Magnetic Materials

Introduction and Context

Magnetostatics concerns magnetic fields produced by steady currents and permanent magnetization distributions. It is the magnetic counterpart of electrostatics, but the analogy is not complete: magnetic field lines form closed loops, isolated magnetic charges have not been observed, and magnetic media exhibit strongly nonlinear and history-dependent behavior. Magnetostatic analysis is therefore central to the design of inductors, transformers, electric machines, relays, actuators, magnetic sensors, shielding systems, and permanent-magnet devices.

The principal field quantities are the magnetic flux density B, measured in tesla, the magnetic field intensity H, measured in ampere per metre, and the magnetization M, also measured in ampere per metre. In a material,

B = μ0(H + M),

where μ0 is the permeability of free space. For a linear, isotropic medium, M = χmH and hence B = μH, where μ = μ0(1 + χm). Real ferromagnetic materials, however, cannot generally be represented by one constant permeability over all operating conditions.

This chapter develops the laws used to calculate magnetic fields, introduces vector magnetic potential and inductance, and then connects field solutions to magnetic circuits, mechanical force, and the behavior of magnetic materials. Particular attention is given to saturation, hysteresis, anisotropy, eddy-current loss, and permanent magnets because these effects determine whether an analytical design remains reliable in practice.

Learning Objectives

After studying this chapter, the reader should be able to:

1. Apply the Biot–Savart law and Ampère’s circuital law to symmetric current distributions.

2. Distinguish B, H, M, magnetic flux, and magnetic scalar or vector potential.

3. Formulate magnetic fields using the vector potential A and appropriate boundary conditions.

4. Calculate inductance, magnetic energy, coenergy, and electromagnetic force.

5. Construct and assess magnetic-circuit models, including leakage and fringing effects.

6. Interpret hysteresis loops, saturation, remanence, coercivity, and magnetic losses.

7. Explain anisotropy and select materials for soft-magnetic and permanent-magnet applications.

8. Identify modeling assumptions, design risks, and appropriate computational tools.

The Magnetostatic Field Equations

For steady currents, charge conservation requires ∇·J = 0, and Maxwell’s magnetostatic equations reduce to

∇·B = 0,

∇×H = Jf,

where Jf denotes free-current density. The first equation states that magnetic flux has no beginning or end. Consequently, magnetic flux lines form closed paths or extend through a mathematical boundary and return elsewhere. The second equation is the local form of Ampère’s law.

The integral form of Ampère’s law is

∮ H·dl = Ienc,

where Ienc is the free current enclosed by the contour. In a vacuum or nonmagnetic medium this may be written as ∮B·dl = μ0Ienc. In magnetic materials, using H is usually more useful because bound-current effects are included through M.

The magnetic flux through a surface S is

Φ = ∫S B·dS.

Since ∇·B = 0, the net flux through any closed surface is zero. This is a major distinction from electric flux, for which a nonzero closed-surface integral may result from enclosed charge.

Biot–Savart Law

The Biot–Savart law gives the magnetic field generated by a current element:

B(r) = (μ0/4π) ∫ [I dl' × (r − r')]/|r − r'|³.

For a volume current distribution, I dl' is replaced by J(r')dv'. The law is particularly useful when geometry is known but Ampère’s law does not provide a simple symmetry.

For an infinitely long straight conductor carrying current I, integration gives

B = μ0I/(2πρ) aφ,

where ρ is the radial distance from the conductor and aφ is the azimuthal unit vector. For a circular loop, the field on its axis at distance z from its centre is

Bz = μ0IR²/[2(R² + z²)^(3/2)].

At the centre, B = μ0I/(2R). These expressions illustrate how geometry determines both magnitude and direction.

Choosing Between the Two Laws

Ampère’s law is efficient when symmetry makes H constant or piecewise constant along a chosen path. It is therefore ideal for infinitely long solenoids, toroids, and straight conductors. Biot–Savart is more general but often leads to more difficult integrations. A common error is to use Ampère’s law as though it automatically implies a uniform field. The law is always valid, but a convenient algebraic extraction of H from the integral is possible only when symmetry justifies it.

For an ideal toroid with N turns carrying current I, the field inside the core at radius ρ is approximately

Hφ = NI/(2πρ),

and Bφ = μHφ. Outside the toroid, the idealized field is zero. Finite cores, leakage, winding spacing, and air gaps produce departures from this result.

Vector Magnetic Potential and Field Formulation

Because ∇·B = 0, B can be represented as the curl of a vector potential A:

B = ∇×A.

This representation automatically satisfies the divergence condition. The vector potential is not unique: replacing A by A + ∇ψ leaves B unchanged. This freedom is called gauge freedom. A frequently used choice is the Coulomb gauge, ∇·A = 0.

In a linear homogeneous medium, combining B = ∇×A with Ampère’s law gives

∇²A = −μJ,

under the Coulomb gauge. In free space, the solution is analogous to the electrostatic potential integral:

A(r) = (μ/4π) ∫ J(r')/|r − r'| dv'.

For filamentary currents, the corresponding volume integral becomes an integral over the current path. The vector potential is especially valuable in two-dimensional and three-dimensional numerical formulations, where solving directly for A can enforce ∇·B = 0 automatically.

At an interface between two media, the normal component of B is continuous:

B1n = B2n.

The tangential component of H changes according to the surface free-current density Kf:

a_n × (H2 − H1) = Kf.

If no free surface current exists, tangential H is continuous. These conditions are fundamental in analyzing air gaps, laminated cores, shielding layers, and permanent-magnet interfaces.

Magnetic Materials: From Constitutive Laws to Physical Behavior

A material responds to an applied magnetic field through magnetization. In a linear isotropic material, B = μH. The relative permeability μr = μ/μ0 indicates how strongly the material increases flux density compared with vacuum. This simple description is adequate for many nonmagnetic materials and for limited operating regions of some magnetic cores, but it is inadequate near saturation or under cyclic excitation.

Materials are often classified as diamagnetic, paramagnetic, or ferromagnetic. Diamagnetic materials develop a weak magnetization opposing the applied field. Paramagnetic materials develop a weak magnetization in the direction of the applied field. Ferromagnetic materials possess strong cooperative magnetic ordering and may retain magnetization after the external field is removed.

Ferromagnetic behavior arises from magnetic domains. Within a domain, atomic magnetic moments are substantially aligned. An applied field causes domain walls to move and domains favorably oriented with the field to grow. At higher fields, rotation of magnetization and progressive alignment occur until saturation is approached.

Hysteresis and Saturation

The B–H curve of a ferromagnetic material depends on its magnetic history. Starting from an unmagnetized state, increasing H produces an initial magnetization curve. When the field is reduced, B does not retrace the same path. At H = 0, the remaining flux density is the remanence Br. A reverse field of magnitude Hc is required to reduce B to zero; Hc is the coercive field.

The area enclosed by a complete hysteresis loop equals the energy dissipated per unit volume per cycle:

Wh = ∮ H dB.

For a sinusoidal excitation, hysteresis loss per unit volume per second is approximately proportional to frequency and loop area. Practical core loss also includes eddy-current loss and excess or anomalous loss. Thin laminations, high-resistivity alloys, and appropriate operating flux densities reduce eddy-current loss.

Saturation occurs when increasing H produces only a small additional increase in B. Near saturation, incremental permeability falls sharply. This causes large magnetizing currents, waveform distortion, increased acoustic noise, and excessive heating in transformers and machines. Designers must distinguish initial, secant, and differential permeability because the relevant value depends on whether a total operating point or a small perturbation is being considered.

Soft and Hard Magnetic Materials

Soft magnetic materials have low coercivity, narrow hysteresis loops, high permeability, and low loss. They are used in transformer cores, inductors, relays, and machine stators. Electrical steel, ferrites, nickel–iron alloys, and cobalt–iron alloys are selected according to frequency, flux density, mechanical requirements, and cost.

Hard magnetic materials have high coercivity and substantial remanence. They are used as permanent magnets because their magnetization is resistant to external disturbance. Alnico, ferrite, samarium–cobalt, and neodymium–iron–boron represent different compromises among energy density, temperature capability, corrosion resistance, mechanical robustness, and cost.

Magnetic Circuits and Their Limits

A magnetic circuit is an approximate network model that represents magnetic flux as flowing through reluctances. For a uniform section,

Reluctance ℜ = l/(μA),

and the magnetomotive force is

F = NI.

The magnetic-circuit analogue of Ohm’s law is

Φ = F/ℜ.

For series sections, reluctances add. For parallel paths, flux divides according to reluctance. This model is valuable for preliminary design because it quickly relates coil turns, current, core dimensions, air gaps, and flux.

Consider a core with one air gap of length lg and cross-sectional area Ag. If leakage and fringing are neglected, the required magnetomotive force is approximately

NI = Hcorelc + Hg lg

= Bcore lc/μcore + Bg lg/μ0.

Because μcore is much larger than μ0, the air gap commonly dominates the required magnetomotive force. This explains why a modest physical gap can substantially reduce inductance and increase coil current.

Magnetic-circuit calculations become unreliable when the core saturates, the air gap is large, fringing is significant, leakage flux is substantial, or the geometry is strongly two- or three-dimensional. Effective gap area can be larger than the physical area because flux spreads at the gap edges. Empirical correction factors may be used, but finite-element analysis is preferable for high-performance designs.

Inductance, Energy, and Coenergy

The flux linkage of a coil is λ = NΦ. For a linear magnetic system, inductance is

L = λ/I = NΦ/I.

For a long solenoid of length l, area A, and N turns,

L = μN²A/l.

A magnetic field stores energy. The energy density in a linear medium is

w = (1/2)B·H.

The total stored energy is

W = ∫V (1/2)B·H dv.

For a linear inductor, W = (1/2)LI². In a gapped magnetic circuit, much of the energy is stored in the air gap rather than in the high-permeability core. This fact is exploited in energy-storage inductors and is also important in calculating actuator force.

For nonlinear systems, energy and coenergy must be distinguished. With current as the independent variable, magnetic coenergy is defined by

W' = ∫0^I λ(i,x) di,

where x is a mechanical coordinate. At constant current, force can be found from

F = (∂W'/∂x)I.

At constant flux linkage, the corresponding energy relation is

F = −(∂W/∂x)λ.

These expressions are not competing physical laws; they are alternative formulations suited to different independent variables. For a linear system, W = W' and both methods give the same force.

For a simple air-gap actuator with approximately uniform gap field, the attractive force is approximately

F ≈ B²A/(2μ0).

The approximation assumes a well-defined gap field and neglects fringing, leakage, saturation, and mechanical deformation. In precision actuator design, the coenergy method with a field solution is more reliable.

Permanent Magnets and Demagnetization

A permanent magnet operates on a portion of its second-quadrant demagnetization curve, where the internal magnetic field opposes the magnet’s remanent flux. Important parameters include Br, Hc, and the maximum energy product (BH)max. A high-energy-product magnet can produce strong fields with relatively small volume, but it may be sensitive to temperature, reverse fields, corrosion, or mechanical damage.

The operating point of a permanent magnet is determined jointly by its material curve and the magnetic circuit surrounding it. Increasing the air gap generally moves the operating point along the demagnetization curve. A magnet can suffer reversible temperature-related reduction in magnetization or irreversible demagnetization if its operating point crosses a critical region.

Permanent-magnet design therefore requires more than selecting the largest Br. The designer must consider coercivity, temperature coefficient, thermal path, short-circuit conditions, fault currents, armature reaction, mechanical retention, and environmental exposure. In rotating machines, rotor stresses and magnet eddy-current heating may be as important as the nominal magnetic field.

Anisotropy and Directional Design

An isotropic material has the same magnetic response in every direction. An anisotropic material does not. In a uniaxial material, the constitutive relation may be represented by a permeability tensor rather than a scalar:

B = [μ]H.

Crystal structure, rolling processes, residual stress, grain orientation, and deliberate manufacturing treatments can create anisotropy. Grain-oriented electrical steel has high permeability and low loss in the rolling direction, making it valuable in transformer limbs. In permanent magnets, texture can produce a preferred easy axis and improve remanence in the desired direction.

Ignoring anisotropy can lead to errors in predicted flux, torque ripple, demagnetization margin, and loss. In numerical models, the principal material axes must be assigned correctly, and tensor data should be used where available.

Applied Scenario: Designing a Gapped Electromagnet

Background

Consider an illustrative electromagnet intended to lift a steel component. The device has a coil around a high-permeability core and a working air gap. The engineering team must estimate current, force, heating, and the risk of saturation.

Situation and Problem

A preliminary magnetic-circuit model predicts adequate lifting force using a low current. However, the core has a narrow section near the coil, and the component surface is not perfectly flat. The team must determine whether the prediction is credible.

Analysis

The first estimate uses

NI = Bcorelc/μcore + Bgaplg/μ0.

Because the gap term dominates, the nominal force is estimated from F ≈ Bgap²Ag/(2μ0). The narrow core section is then checked using B = Φ/A; its smaller area produces a larger local flux density and may approach saturation. If saturation occurs, the assumed constant μcore becomes invalid and the required current increases sharply.

The imperfect contact introduces an additional effective gap. Surface roughness and tilt also reduce the useful area. A thermal calculation must evaluate copper loss I²R and the duty cycle. Finally, the release condition is considered: residual magnetism may prevent the load from separating when current is removed.

Approach and Solution

The team first uses a reluctance model for sensitivity studies, varying gap length, current, and core area. It then performs a nonlinear finite-element analysis using the measured B–H curve, including the real pole geometry and fringing. A transient thermal model checks continuous and intermittent operation. A demagnetizing or reverse-voltage circuit is added so that the magnet can release reliably.

Outcome

The refined model predicts less force than the ideal calculation but identifies a modest increase in pole area and a controlled current boost as sufficient. The design is accepted only after verifying temperature rise, mechanical retention, fault behavior, and release time experimentally.

Learning Points

The example illustrates that the air gap usually controls magnetomotive force, while local core constrictions control saturation. Ideal magnetic circuits are excellent for first estimates, but fringing, nonlinear material data, thermal limits, and operational transients determine the final design.

Practical and Professional Applications

Magnetostatic methods support transformer and inductor sizing, motor and generator torque prediction, magnetic bearings, solenoid valves, contactors, MRI-gradient components, Hall-effect sensors, current transformers, magnetic separators, and electromagnetic compatibility studies. In industry, engineers combine analytical estimates with material databases, prototype measurements, and finite-element simulation.

A sound workflow is to define the current sources and geometry, exploit symmetry, select the simplest valid analytical model, identify regions likely to saturate, estimate energy and force, and then validate critical assumptions numerically or experimentally. Measurements should include not only nominal flux but also temperature, waveform distortion, inductance versus current, and behavior under fault or demagnetizing conditions.

Challenges, Limitations, and Common Mistakes

Several errors recur in magnetostatic design. Treating μr as constant through the entire operating range can greatly overestimate flux. Applying Ampère’s law without adequate symmetry produces unjustified field values. Confusing B and H leads to incorrect force, energy, and material calculations. Ignoring winding leakage and air-gap fringing causes inductance and force predictions to be optimistic.

Another common mistake is to use a closed hysteresis-loop loss estimate for a minor-loop operating condition. Core loss depends on waveform, frequency, temperature, DC bias, and magnetic history. Permanent magnets are sometimes selected solely by remanence, while coercivity and temperature stability are neglected. Finally, a two-dimensional model may miss end effects, skew, leakage, or axial force.

Professional, Safety, and Governance Considerations

Magnetic devices must be designed for abnormal as well as nominal conditions. Stored magnetic energy can create sudden mechanical motion, high short-circuit currents, and strong attraction of ferromagnetic tools. Permanent magnets can damage instruments, pinch fingers, and interfere with medical implants or data-storage equipment. High-current coils require insulation coordination, thermal protection, and controlled discharge paths.

Engineering documentation should state material data sources, temperature limits, tolerances, model assumptions, mesh or convergence criteria, and validation results. Where a design affects safety, lifting, medical equipment, or grid infrastructure, independent review and applicable standards are essential. Material substitution should not be approved merely because two materials have similar nominal permeability; their loss, saturation, coercivity, temperature response, and manufacturing condition may differ substantially.

Computational Tools and Emerging Directions

Finite-element magnetostatic solvers calculate fields in complex geometries and can incorporate nonlinear B–H curves, anisotropic tensors, permanent-magnet remanence, motion coupling, and thermal feedback. Automated parameter sweeps and optimization can evaluate pole shapes, winding placement, or magnet dimensions. Reduced-order models and surrogate models are increasingly used to accelerate design exploration, but they must be trained and validated against trustworthy field solutions or measurements.

Digital monitoring of inductance, current, temperature, and vibration can reveal saturation, insulation deterioration, or mechanical misalignment. Data-driven methods may assist fault diagnosis, yet physical constraints remain important: an algorithm that fits nominal measurements may fail under a changed air gap, altered material batch, or abnormal temperature. The best practice is therefore a hybrid approach combining Maxwell-based models, measured material data, uncertainty analysis, and engineering judgment.

Best Practices for Magnetostatic Design

1. Define whether the problem is linear, nonlinear, transient, or effectively magnetostatic before selecting equations.

2. Use symmetry only when the sources, geometry, and material properties share that symmetry.

3. Distinguish B, H, M, flux, flux linkage, energy, and coenergy in calculations.

4. Use measured B–H and loss data over the actual temperature and frequency range.

5. Check saturation at narrow sections, corners, teeth, and pole tips rather than only in the main core.

6. Include fringing, leakage, manufacturing tolerances, and temperature in final verification.

7. Validate critical simulations with analytical bounds, prototype measurements, or both.

8. Evaluate normal operation, start-up, short circuit, demagnetization, release, and fault conditions.

Chapter Summary

Magnetostatics describes fields generated by steady currents and permanent magnetization. Biot–Savart law provides a general field integral, while Ampère’s law is especially efficient for highly symmetric geometries. The divergence-free nature of B permits the introduction of the vector magnetic potential A, which is central to analytical and finite-element formulations.

Magnetic circuits approximate flux paths using magnetomotive force and reluctance, but their accuracy is limited by leakage, fringing, saturation, and complex geometry. Inductance follows from flux linkage, while magnetic energy and coenergy provide systematic methods for calculating stored energy and force. Magnetic materials are characterized not only by permeability but also by hysteresis, remanence, coercivity, saturation, loss, anisotropy, and temperature dependence. Soft magnets are optimized for low-loss reversible magnetization, whereas hard magnets are designed to resist demagnetization. Reliable engineering practice combines field theory, material data, thermal and mechanical analysis, numerical computation, and experimental validation.

Key Terms

Magnetic flux density B: Magnetic flux per unit area, measured in tesla.

Magnetic field intensity H: Field quantity related to free current and measured in ampere per metre.

Magnetization M: Magnetic dipole moment per unit volume of a material.

Permeability μ: Ratio relating B to H in a linear material.

Magnetic flux Φ: Surface integral of B through a specified area.

Biot–Savart law: Integral law giving the magnetic field produced by current elements.

Ampère’s law: Relation stating that the circulation of H equals enclosed free current.

Vector magnetic potential A: Vector field whose curl equals B.

Magnetic scalar potential: A potential representation applicable in current-free regions.

Magnetomotive force: Quantity NI that drives magnetic flux around a magnetic circuit.

Reluctance: Opposition to magnetic flux, equal to l/(μA) for a uniform section.

Inductance: Ratio of flux linkage to current under the specified operating conditions.

Saturation: Region in which increases in H produce relatively small increases in B.

Hysteresis: Dependence of magnetic state on the previous history of applied field.

Remanence: Flux density remaining when applied H is returned to zero.

Coercivity: Reverse field required to reduce magnetization or B to a specified zero condition.

Anisotropy: Directional dependence of magnetic properties.

Fringing: Spreading of magnetic flux near an air-gap boundary.

Magnetic coenergy: Energy-like quantity useful for force calculations at constant current.

Maximum energy product: Maximum value of the product BH on a permanent-magnet demagnetization curve.

Review Questions

Short Answer Questions

1. State the differential and integral forms of Ampère’s law for magnetostatics.

2. Why is the vector potential defined by B = ∇×A?

3. Distinguish remanence from coercivity.

4. What is magnetic reluctance, and how is it related to magnetomotive force?

5. Why does an air gap usually dominate the magnetomotive-force requirement?

Long Answer Questions

1. Derive the magnetic field around an infinitely long straight conductor using Ampère’s law and explain the symmetry used.

2. Explain the physical origin of hysteresis, saturation, and core loss in a ferromagnetic material.

3. Develop the magnetic-circuit model of a gapped core and discuss the effects neglected by the ideal model.

4. Explain magnetic energy and coenergy and derive their use in calculating electromagnetic force.

5. Compare soft-magnetic and hard-magnetic materials with reference to applications, material parameters, and design constraints.

Analytical / Application Questions

1. A coil has N turns around a core containing a dominant air gap. Formulate the approximate current required for a specified gap flux density and identify the assumptions behind the result.

2. A finite-element model predicts a high flux density at a narrow pole tip while the average core flux density is moderate. Explain the likely consequences and propose design changes.

3. A permanent-magnet actuator loses force at elevated temperature and fails to release after a reverse-current fault. Analyze the possible magnetic causes and recommend a verification plan.

22. Multiple Choice Questions with Answers

Question 1

Which law is most directly suited to finding the magnetic field inside an ideal long toroid?

A. Gauss’s law for electric fields

B. Ampère’s circuital law

C. Faraday’s law of induction

D. Coulomb’s law

Correct Answer: B. Ampère’s circuital law

Question 2

What condition is automatically satisfied when magnetic flux density is written as B = ∇×A?

A. ∇·B = 0

B. ∇×B = 0

C. ∇·H = 0

D. ∇×A = 0

Correct Answer: A. ∇·B = 0

Question 3

For a uniform magnetic-circuit section, which expression represents reluctance?

A. μA/l

B. l/(μA)

C. NIμA/l

D. B/(μH)

Correct Answer: B. l/(μA)

Question 4

The area enclosed by a complete B–H hysteresis loop represents

A. remanent flux density

B. coercive field strength

C. energy dissipated per unit volume per cycle

D. magnetic susceptibility

Correct Answer: C. energy dissipated per unit volume per cycle

Question 5

Why can a narrow core section become saturated even when the average core flux density is acceptable?

A. Flux is zero at narrow sections

B. The same flux passes through a smaller cross-sectional area

C. Permeability becomes infinite there

D. Air-gap reluctance disappears

Correct Answer: B. The same flux passes through a smaller cross-sectional area

Question 6

For a linear inductor carrying current I, the stored magnetic energy is

A. LI

B. L/I

C. (1/2)LI²

D. (1/2)I²/L

Correct Answer: C. (1/2)LI²

Question 7

Which property is generally desirable in a transformer core material?

A. High coercivity and a wide hysteresis loop

B. Low permeability and high remanence

C. Low coercivity, high permeability, and low loss

D. Very high electrical conductivity to increase eddy currents

Correct Answer: C. Low coercivity, high permeability, and low loss

Question 8

At an interface without free surface current, which field component is continuous?

A. The tangential component of H

B. The normal component of H always

C. The tangential component of B always

D. The magnitude of B always

Correct Answer: A. The tangential component of H

Question 9

The force on a simple uniform air gap is approximately proportional to

A. B²A/(2μ0)

B. BAμ0/2

C. B/(2μ0A)

D. 2μ0A/B²

Correct Answer: A. B²A/(2μ0)

Question 10

Why is a permanent magnet not selected solely on the basis of high remanence?

A. Remanence is unrelated to magnetic field

B. Coercivity, temperature stability, and demagnetization resistance are also critical

C. High remanence always causes zero force

D. Permanent magnets operate only in the first quadrant

Correct Answer: B. Coercivity, temperature stability, and demagnetization resistance are also critical