[Getdp] Re: (no subject)
Johan Gyselinck
Johan.Gyselinck at ulg.ac.be
Fri Nov 15 10:16:30 CET 2002
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>