<!doctype html public "-//w3c//dtd html 4.0 transitional//en">

<html>

Hi Janvrot,

<blockquote TYPE=CITE>Thank you very much for the files. I would not be

able to do it by my

<br>own.</blockquote>

You're welcome.

<blockquote TYPE=CITE> 

<br>I took some time to figure out how to compute the values for another

<br>magnetization curve, and I guess I got it. I ran the pro putting a

<br>magnet on, taking the inductor off and results seems ok.</blockquote>

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

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

<blockquote TYPE=CITE> 

Nevertheless, I'm having problems running in 3D, my target problem.

Is

this ok for 3D?

<br>I get

<br>Warning   : Null right hand side in linear system

<br>Error     : Zero row encountered in ILU</blockquote>

The source term (due to the permanent magnet?) is not properly written,

I guess. This may produce a zero RHS.

<br> 

<blockquote TYPE=CITE> 

<br>Attached you may find a complete pro file (without includes) and

<br>solver.par</blockquote>

This is still for the 2D problem, I guess. Can you send the new geo file,

so that I can run your stuff?

<p>Johan

<br> 

<br> 

<blockquote TYPE=CITE> 

<br>Thanks again

<br>Janvrot

<p>Johan Gyselinck wrote:

<br>>

<br>> Hello Janvrot,

<br>>

<br>> I apologize for the delay. I've been pretty busy lately, like all

of

> us, I guess. But I've finally done my homework.

<br>>

<br>> I've added some lines in the test.pro you've sent me a while ago,

for

<br>> the nonlinear magnetic vector potential case, using the Newton-Raphson

<br>> method and the differential reluctivity tensor.  It seems to

work, but

<br>> the pro-file is a bit messy.  I'm planning to make a well documented

<br>> 2D example and make it available to all GetDP users in the very near

> (let's be optimistic!) future.

<br>>

> If you have any questions or if you need more help, just give a ring.

> I'll try to answer much quicker.

<br>>

<br>> Bye,

<br>>

<p>--

<br>Janvrot IVM

<br>janvrot@infolink.com.br

<p>  ------------------------------------------------------------------------

<br>                  

Name: nonlin_a.pro

<br>   nonlin_a.pro    Type: unspecified type

(application/octet-stream)

<br>              

Encoding: base64

<p>                

Name: solver.par

<br>   solver.par    Type: unspecified type (application/octet-stream)

<br>            

Encoding: base64</blockquote>

<pre>-- 

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 -- <A HREF="http://elap.montefiore.ulg.ac.be/">http://elap.montefiore.ulg.ac.be/</A></pre>

 </html>