[Getdp] Re: (no subject)

[Johan Gyselinck]

Hi Janvrot,

> Thank you very much for the files. I would not be able to do it by my
> own.

You're welcome.

>
> I took some time to figure out how to compute the values for another
> magnetization curve, and I guess I got it. I ran the pro putting a
> magnet on, taking the inductor off and results seems ok.

  hc[ Magnet ] = Vector[0., 0., -920000.];  ???

In a 2D problem, hc must be in the xy plane!

>
> Nevertheless, I'm having problems running in 3D, my target problem. Is
> this ok for 3D?
> I get
> Warning   : Null right hand side in linear system
> Error     : Zero row encountered in ILU

The source term (due to the permanent magnet?) is not properly written, I
guess. This may produce a zero RHS.


>
> Attached you may find a complete pro file (without includes) and
> solver.par

This is still for the 2D problem, I guess. Can you send the new geo file, so
that I can run your stuff?

Johan



>
> Thanks again
> Janvrot
>
> Johan Gyselinck wrote:
> >
> > Hello Janvrot,
> >
> > I apologize for the delay. I've been pretty busy lately, like all of
> > us, I guess. But I've finally done my homework.
> >
> > I've added some lines in the test.pro you've sent me a while ago, for
> > the nonlinear magnetic vector potential case, using the Newton-Raphson
> > method and the differential reluctivity tensor.  It seems to work, but
> > the pro-file is a bit messy.  I'm planning to make a well documented
> > 2D example and make it available to all GetDP users in the very near
> > (let's be optimistic!) future.
> >
> > If you have any questions or if you need more help, just give a ring.
> > I'll try to answer much quicker.
> >
> > Bye,
> >
>
> --
> Janvrot IVM
> janvrot at infolink.com.br
>
>   ------------------------------------------------------------------------
>                    Name: nonlin_a.pro
>    nonlin_a.pro    Type: unspecified type (application/octet-stream)
>                Encoding: base64
>
>                  Name: solver.par
>    solver.par    Type: unspecified type (application/octet-stream)
>              Encoding: base64

--
Dr. ir. Johan Gyselinck
University of Liege
Dept. of Electrical Engineering (ELAP), Institut Montefiore
Campus du Sart Tilman, Parking P32, Bâtiment B28, B-4000 LIEGE, Belgique
phone: +32-4-3663732 -- fax: +32-4-3662910 -- http://elap.montefiore.ulg.ac.be/


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20021115/6ba81571/attachment.html>