<br><font size=2 face="sans-serif">Dear GetDP users:</font>
<br>
<br><font size=2 face="sans-serif">I just started learning about magnetism.
 I have a couple of basic questions that I am not able to figure out.
 It will be appreciated if someone here can shed some light on this.</font>
<br>
<br><font size=2 face="sans-serif">I have a permanent magnet - no electric
current in the domain.  I would like to calculate the magnetic flux
density.</font>
<br>
<br><font size=2 face="sans-serif">I have found a 3D code posted by Jasper
in 2006.  The questions are:</font>
<br>
<br><font size=2 face="sans-serif">1. The governing equation for magnetostatics
filed is:</font>
<br><font size=2 face="sans-serif">    Curl ( nu * Curl A) =
J, where A is vector potential and J is electric current density.</font>
<br><font size=2 face="sans-serif">    Now, since there is no
current, I will think that J is zero.</font>
<br><font size=2 face="sans-serif">2. In the Formulation of getdp code
(2D Magnetostatics example on getdp wiki),</font>
<br><font size=2 face="sans-serif">    it showed:</font>
<br><font size=2 face="sans-serif">       Equation
{</font>
<br><font size=2 face="sans-serif">           
        Galerkin { [nu[] * Dof{ Curl a}, {Curl a}];
In Domain; Jacobian JVol; Integration I1; }</font>
<br><font size=2 face="sans-serif">           
        Galerkin { [hc[], {Curl a}; In Domain_M; Jacobian
JVol; Integration I1;}</font>
<br><font size=2 face="sans-serif">        }</font>
<br>
<br><font size=2 face="sans-serif">How does {hc[], {Curl a}} come into
the equation? I believe that hc is Coercitivity of the magnet?</font>
<br>
<br><font size=2 face="sans-serif">Thanks in advance for your reply!</font>
<br>
<br><font size=2 face="sans-serif">Best Regards,</font>
<br>
<br><font size=2 face="sans-serif">Pei</font>