[Getdp] Nonlinear BH function

[Christophe Geuzaine]

janvrot at infolink.com.br wrote:
> 
> Hi, All
> 
> I'm using a function to express a nonlinear BH curve.
> 
> So, I have mu as a function of H in a Magnetic Static Scalar Potential
> formulation. Is that ok, or it shoul be the other way? Here is the code,

Yes, that's the way to do it. Are you sure that Domain, Steel and Magnet
do not overlap? It would be better to write:

Domain = Region[{NonSteel,Steel}];
mu[NonSteel] = mu0;
mu[Steel] = 1./(as + bs * Norm[$1]);

and simply

Galerkin { [ - mu[{d phi}] * Dof{d phi} , {d phi} ];
             In Domain; Jacobian MyJac; Integration MyInt; }

> and the error message.
> 
> Function {
>   mu0 = 4.e-7 * Pi ;
>   murMagnet = 1.1;
>   as = 417.497;
>   bs = 38.3881;
>   mu[Steel] = 1./(as + bs * Norm[$1]);
> 
>   mu[Magnet]  = murMagnet * mu0;
>   mu[Air]     = mu0;
>   hc[ Magnet1 ] = Vector[0., 0.,  10000.];
>   hc[ Magnet2 ] = Vector[0., 0., -10000.];
> }
> 
> Formulation {
>   { Name MagSta_phi; Type FemEquation;
>     Quantity {
>       { Name phi; Type Local; NameOfSpace Hgrad_phi; }
>     }
>     Equation {
>       Galerkin { [ - mu[] * Dof{d phi} , {d phi} ];
>                  In Domain; Jacobian MyJac; Integration MyInt; }
>       Galerkin { [ - mu[{d phi}] * Dof{d phi} , {d phi} ];
>                   In Steel; Jacobian MyJac; Integration MyInt; }
>       Galerkin { [ - mu[] * hc[] , {d phi} ];
>                  In Magnet; Jacobian MyJac; Integration MyInt; }
>     }
>   }
> }
> 
> I get the message:
> P r o c e s s i n g . . .
> Operation : Generate[A]
> Solver    : Loading parameter file 'solver.par'
> Error     : Wrong type of argument for function 'Norm'
> 
> Thanks for any help.
> --
> Janvrot IVM
> janvrot at infolink.com.br
> _______________________________________________
> GetDP mailing list
> GetDP at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp

-- 
Christophe Geuzaine

Tel: (626) 395-4552    http://www.geuz.org
Fax: (626) 578-0124    mailto:geuzaine at acm.caltech.edu