Halbach Array
Theory and Physics
What is a Halbach Array?
Professor, what's the difference between a Halbach array and a regular magnet arrangement?
It's an arrangement that concentrates magnetic flux on one side by gradually rotating the magnetization direction. In an ideal Halbach array, the magnetic flux density on one side doubles, while it becomes zero on the opposite side.
$d$: magnet thickness, $\lambda$: wavelength of magnetization pattern, $M$: number of magnet segments per wavelength.
M=4 is common (rotating 90 degrees each), right?
Correct. For M=4, $\sin(\pi/4)/(\pi/4) \approx 0.90$. Increasing M brings it closer to ideal, but increases assembly precision requirements and part count. Used in linear motors, magnetic levitation, and Wiggler magnets.
Summary
- One-Side Concentration — Enhances magnetic flux density on one side up to 2x
- $M$ Segments — M=4 (90-degree rotation) is the standard configuration
- Applications — Linear motors, SPM-type motors, particle accelerators
Halbach Array—The Magic of Doubling the Magnetic Field on One Side Just by "Rotating" Magnet Arrangement
The Halbach array is a permanent magnet arrangement invented by Klaus Halbach (1980, Lawrence Berkeley Laboratory) as an undulator for accelerators. By arranging magnets with their orientations rotated 90° each, the magnetic field on one side reinforces (theoretically doubling), while on the opposite side it cancels out to nearly zero. This "self-shielding" property eliminates the need for a back iron yoke, allowing linear motors, MRI, and magnetic levitation systems to be lighter with the same magnet volume. In CAE, FEM magnetic field analysis is used to optimize the number of divisions and magnetization angles of the Halbach array, evaluating its approximation to the theoretical value.
Physical Meaning of Each Term
- Electric Field Term $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$: Faraday's law of electromagnetic induction. A time-varying magnetic flux density generates an electromotive force. [Everyday Example] A bicycle dynamo (generator) produces voltage in a nearby coil by rotating a magnet—a direct application of this law where a changing magnetic field induces an electric field. Induction cooking (IH) heaters also use the same principle, where high-frequency magnetic field changes induce eddy currents in the pot bottom, heating it via Joule heat.
- Magnetic Field Term $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$: Ampère-Maxwell's law. Electric current and displacement current generate a magnetic field. [Everyday Example] Passing current through a wire creates a magnetic field around it—this is Ampère's law. Electromagnets operate on this principle, creating a strong magnetic field by passing current through a coil. Smartphone speakers also apply this law: current → magnetic field → force on the diaphragm. At high frequencies (e.g., GHz-band antennas), the displacement current $\partial D/\partial t$ cannot be ignored, describing electromagnetic wave radiation.
- Gauss's Law $\nabla \cdot \mathbf{D} = \rho_v$: Indicates that electric charge is the divergence source of electric flux. [Everyday Example] Rubbing hair with a plastic sheet creates static electricity, making hair stand up—the charged sheet (electric charge) radiates electric field lines outward, exerting force on the light hair. Capacitor design uses this law to calculate the electric field distribution between electrodes. ESD (electrostatic discharge) countermeasures are also based on electric field analysis following Gauss's law.
- Magnetic Flux Conservation $\nabla \cdot \mathbf{B} = 0$: Expresses the non-existence of magnetic monopoles. [Everyday Example] Cutting a bar magnet in half does not create a magnet with only an N pole or only an S pole—N and S poles always exist as a pair. This means magnetic field lines form "closed loops with no start or end points." In numerical analysis, the formulation using vector potential $\mathbf{B} = \nabla \times \mathbf{A}$ is used to satisfy this condition automatically, ensuring magnetic flux conservation.
Assumptions and Applicability Limits
- Linear material assumption: Permeability and permittivity are independent of magnetic/electric field strength (nonlinear B-H curve needed in saturation region)
- Quasi-static approximation (low frequency): Displacement current term can be ignored ($\omega \varepsilon \ll \sigma$). Common in eddy current analysis.
- 2D assumption (cross-section analysis): Effective when current direction is uniform and end effects can be ignored.
- Isotropic assumption: For anisotropic materials (e.g., silicon steel rolling direction), direction-specific property definitions are needed.
- Non-applicable cases: Additional constitutive relations are needed for plasma (ionized gas), superconductors, nonlinear optical materials.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Magnetic Flux Density $B$ | T (Tesla) | 1T = 1 Wb/m². Permanent magnets: 0.2〜1.4T |
| Magnetic Field Strength $H$ | A/m | Horizontal axis of B-H curve. Conversion from CGS Oe (Oersted): 1 Oe = 79.577 A/m |
| Current Density $J$ | A/m² | Calculated from conductor cross-section and total current. Note non-uniform distribution due to skin effect. |
| Permeability $\mu$ | H/m | $\mu = \mu_0 \mu_r$. In vacuum $\mu_0 = 4\pi \times 10^{-7}$ H/m |
| Electrical Conductivity $\sigma$ | S/m | Copper: approx. 5.96×10⁷ S/m. Decreases with temperature rise. |
Numerical Methods and Implementation
Modeling with FEM
How do you model a Halbach array in FEM?
Assign a different magnetization direction vector $\mathbf{M}_0$ to each magnet segment. In JMAG or Maxwell, define a material coordinate system for each segment and specify the direction of remanent magnetization.
Should the gaps between segments be included in the model?
In actual machines, there is an adhesive layer (0.05〜0.2 mm). This gap reduces magnetic flux density by a few percent, so it should be included for precise design. Using 2D periodic boundary conditions allows modeling just one period.
Summary
- Set magnetization direction individually for each segment — Specify via material coordinate system
- Adhesive layer gap — Affects performance by a few percent
- Periodic Boundary Conditions — Reduces computational cost
The "Magnetization Direction Setting Mistake" Often Made When Implementing Halbach Arrays in FEA
When modeling a Halbach array in FEA, the most common mistake is setting the direction of the magnetization vector for each magnet. While each magnet's orientation rotates slightly within the array, confusing local and global coordinate systems quickly scrambles the magnetic field. The correct procedure is: "Calculate the angle of each magnet relative to the array's symmetry axis beforehand, create a list of magnetization vectors as (Mx, My) = Br×(cos θ, sin θ), and then input them into the model." Creating an angle table in Excel makes it easier to handle parametric changes.
Edge Elements (Nedelec Elements)
Elements specialized for electromagnetic field analysis. Automatically guarantee continuity of tangential components and eliminate spurious modes. Standard for 3D high-frequency analysis.
Nodal Elements
Used for scalar potential formulations. Effective for scalar potential methods in magnetostatics and electrostatic field analysis.
FEM vs BEM (Boundary Element Method)
FEM: Handles nonlinear materials and inhomogeneous media. BEM: Naturally handles infinite domains (open boundary problems). Hybrid FEM-BEM is also effective.
Nonlinear Convergence (Magnetic Saturation)
Nonlinearity of B-H curve handled by Newton-Raphson method. Residual criterion: $||R||/||R_0|| < 10^{-4}$ is common.
Frequency Domain Analysis
Reduced to a steady-state problem using time-harmonic assumption. Requires complex number operations, but broadband characteristics are obtained via time-domain analysis.
Time Domain Time Step
Time step less than 1/20 of the highest frequency component is required. Implicit time integration allows larger steps but requires attention to accuracy.
Choosing Between Frequency Domain and Time Domain
Frequency domain analysis is like "tuning a radio to a specific frequency"—it can efficiently calculate the response at a single frequency. Time domain analysis is like "recording all channels simultaneously"—it can reproduce transient phenomena containing all frequency components, but at a higher computational cost.
Practical Guide
Design in Practice
Related Topics
なった
詳しく
報告