## Quick Intro to Permanent Magnet Modeling

##### Fanny Littmarck | June 21, 2013

I’ve written several blog entries involving permanent magnets, in one way or another. Reading those may have raised the question “what about a more simple introduction to permanent magnet simulations?” Fair enough, here’s how to model a permanent magnet and its surrounding magnetic field.

*Horseshoe permanent magnet and an iron rod.*

### About Magnetostatics

Before we dive into the permanent magnet modeling, let’s quickly go over the general concept of magnetostatics. Simply put, magnetostatics involves studying magnetic fields in systems where the magnetic currents do not change with respect to time (they’re steady). This is contrary to the eddy current brakes we discussed a while back, where the magnetic field does change with time. The particular permanent magnet model we will go through today turns out to be a great example of a magnetostatic problem.

### Setting up the Permanent Magnet Model

Let’s suppose we want to model the magnetic field that surrounds a horseshoe permanent magnet. In accordance with the image at the top of this post, we also want to have a look at how the magnet interacts with a metal rod. When modeling the magnetic field and computing the force that acts on the rod, we can take advantage of the geometry’s symmetry and the magnetic field’s antisymmetry. Although this isn’t a computation-heavy problem, it is a good example for introducing the concept of symmetry. Knowing when to take advantage of symmetry can save you a lot of time, and in the case of our simple magnet example, only a quarter of the geometry then needs to be modeled. As an advocate for efficiency, I’d say that’s very nice, indeed.

Using COMSOL Multiphysics along with the AC/DC Module, we can get started with defining the magnetic scalar potential, *V*_{m}, by selecting the *Magnetic Fields, No Currents* interface. Doing so prompts the model to use the appropriate equation: –∇ ⋅ (μ_{0}∇*V*_{m} – μ_{0}M_{0}) = 0.

Next, we want to import the permanent magnet geometry from the Model Library. The geometry includes both the horseshoe magnet and the metal rod, but as you can see, it shows the components in their entirety:

I mentioned we could take advantage of our problem’s symmetry and antisymmetry to speed up our simulations, and now is the time to do so. We can use the *xy*- and *xz*-planes to form our geometry’s exterior boundaries, as they correspond with the magnetic field’s symmetry and antisymmetry, respectively. To do that, we will need to create an air box, and then delete the parts of the geometry we no longer want to work with (complete instructions on how to perform this and all other steps can be found in the permanent magnet’s model documentation).

The magnetic field is tangential to the boundary on the *xy*-plane, and perpendicular to the boundary on the *xz*-plane, which is described by the *Magnetic Insulation* and *Zero Magnetic Scalar Potential* conditions, respectively.

### Modeling the Magnetic Field and Calculating the Force on the Metal Rod

If we plot the magnetic flux density, we can visualize the magnetic field that surrounds our permanent magnet. By also adding an arrow plot, we can view the direction of the field.

Finally, we can compute the force the permanent magnet exerts on the rod, by using a Global Evaluation. Upon doing that, we learn that the magnet exerts a force of 1.53 N on the metal rod. Remember, we only modeled a quarter of our geometry, so the actual force on the rod is therefore four times this value, or 6.11 N.

### Further Reading

- Certified Consultants 29
- Chemical 44
- COMSOL Now 130
- Conference 98
- Core Functionality 55
- Electrical 113
- AC/DC 42
- MEMS 13
- Plasma 4
- RF 29
- Semiconductor 2
- Wave Optics 6

- Fluid 85
- CFD 39
- Microfluidics 8
- Mixer 2
- Molecular Flow 8
- Pipe Flow 5
- Subsurface Flow 10

- Interfacing 34
- Mechanical 135
- Multipurpose 14
- Tips & Tricks 31
- Trending Topics 29
- User Perspectives 78
- Video 62

## Comments

Hi, Fanny Littmarck. Thanks for sharing. I am a beginner on Comsol, and I have build a first comsol model on permanent magnet referencing ‘The Magnetic Field from a Permanent Magnet’ . My interest is to check the magnetic density of the magnetic circuit I built with a permanent magnet and several pole pieces. However, the magnetic field is restricted in the permanent domain I have defined. Can you help figure it out? I really appreciate it. Cause I cannot post image or file here, may I post a post in the community or I send you an email? My email is lqwei7@sjtu.edu.cn

Thanks!!!

Hi Charve,

Happy to hear you are using COMSOL Multiphysics. For help with your specific case, please send your file and question(s) to support@comsol.com.