(adsbygoogle = window.adsbygoogle || []).push({}); 1. The problem statement, all variables and given/known data

I have to build a 2D model of simple magnetic circuit, that is iron core (2) with coil(1) and air gap(3) applied to it. Model has to be solved on numerical basis, by finite difference or element method (whichever is easier to implement). Number 4 on drawing represents boundary condition, that would be magnitude of magnetic field density B = 0.

http://img813.imageshack.us/img813/4049/trafon.th.jpg [Broken]

2. Relevant equations

What equations to use and how properly discretize them?

3. The attempt at a solution

I started with differential form of Ampere law, which is:

curl(H) = J, where H is vector of magnetic field strength and J is vector of current density.

Because B has only 2 space components (third iz zero), curl(H) translates to:

(dHy/dx - dHx/dy) *ez = J * ez, where Hy is y component of H, Hx is x component of H and ez is the direction in z-axis.

ez falls out, so final form is:

dHy/dx -dHx/dy = J

First derivatives are substituted with differences, i'll use central difference approximation:

du/dx = (u_{(i+1)}- u_{(i-1)})/(2*deltax)

Hy_{(i+1,j)}- Hy_{(i-1,j)}- Hx_{(i,j+1)}+ Hx_{(i ,j-1)}= J*2*deltah

Grid is equidistant, so for both dimensions deltah is same.

I used the formula above for the coil region, elsewhere I set equation to (where no current density):

Hy_{(i+1,j)}- Hy_{(i-1,j)}- Hx_{(i,j+1)}+ Hx_{(i ,j-1)}= 0

I got some strange and most possible, wrong results.

thank you in advance

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Numerical model of simple magnetic circuit

Can you offer guidance or do you also need help?

**Physics Forums | Science Articles, Homework Help, Discussion**