9
$\begingroup$

I'm trying to simulate Maxwell equations by finite elements method for a current carrying wire. My 3d geometry consists of a cylinder and a box containing it. I will use a mixed formulation and Nedelec's elements introducing a vector potential. I'm in a magnetostatic regime. At the end of the simulation I would like to plot the lines of the magnetic field around the wire (the cylinder) and probably compute forces and see that the numerical method agrees with classical results. My troubles concern the boundary conditions I have to impose on the surfaces of the cylinder (side, top and bottom) and on the box. I think I will use a value of permittivity of $\mu_0$ for the box and a copper relative one for the wire. I hope someone of you can help me seeing that my physics background is not so good.

$\endgroup$
3
  • 2
    $\begingroup$ Maybe better suited to be posted in Computational Science $\endgroup$ Feb 27, 2014 at 13:42
  • 1
    $\begingroup$ @ja72 I think this topic is fine for Physics SE. $\endgroup$
    – DanielSank
    May 23, 2017 at 0:29
  • 3
    $\begingroup$ The way this site works is that you ask a question and other users answer. This post does not seem to ask a question. Phrase edit it. $\endgroup$
    – DanielSank
    May 23, 2017 at 0:30

2 Answers 2

0
$\begingroup$

It is easier to answer if you have a sketch of the problem you want to solve.

I think that good results can be obtained only by setting the outer space section large enough and giving no boundary conditions at the outer boundary (which is equivalent to giving $\bf{n}\times\bf{H}=\bf{0}$ at the outer boundary).

[Edit #1] A similar problem was solved numerically. Centered cubic iron (assumed linear material having relative permeability 1500) and circular coil. Boundary conditions ($\vec{n}\times \vec{A}=0$) are applied to the x=0 and y=0 planes to meet symmetry. Boundary conditions

Numerically calculated magnetic B fields and vector A potentials are shown. Magnetic B fields Vector potential A

$\endgroup$
0
$\begingroup$

Some updates....I found out I will have to impose boundary conditions like: $$\mu^{-1}\nabla\times \textbf{A} \times \textbf{n} =\textbf{0}$$ in the outer boundary and in the interface. Those correspond to the Neumann natural conditions of my weak formulation and those imply that in the two boundary pieces, $\textbf{H} \times \textbf{n} = \textbf{0}$. Furthermore, I need to impose $$ \textbf{B} \cdot \textbf{n} =0 $$ just at the interface and this in terms of the vector potential formulation corresponds to $\textbf{A} \times \textbf{n} = \textbf{0}$. This last corresponds to a Dirichlet boundary condition.

The problem is that even if is true, I don't know how to implement the Dirichlet boundary conditions with FreeFem++.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.