[Getdp] Magnetodynamic 3D - Electromagnetic circuit element

[Christophe Geuzaine]

Florin CIUPRINA wrote:
> Hello Christophe,
> 
> 
>>You should try to increase Nb_Fill: this will make the incomplete LU
>>decomposition (used as a preconditioner) be closer and closer to the
>>complete LU. This should improve convergence a lot.
> 
> 
> Thanks. Indeed, this has improved the convergence a lot.
> 
> I tried several resolutions for different frequencies between 10 and 1.e5 Hz.
> Looking at the results I noticed that, when the frequency increases, the J 
> vectors
> "flow" closer to one L-shaped surface, so looking along X axis there is an 
> asymmetry,
> the J field lines preferring the L-shape surface at z=0. But from the 
> symmetry point
> of view, J should have the same behaviour near z=zmax.
> With other words I have a symmetrical problem with respect to the geometry
> and boundary conditions (a=0 on all the boundary of the conductor,
> and v=1 on one terminal and v=0 on the other terminal) and the distribution 
> of J vectors is
> asymmetrical when the frequency increases. I refined the mesh for having 
> element size lower
> than the skin depth, but the  results are still asymmetrical.
> I have no idea where this assymmetry comes from, therefore I attach again 
> the problem,
> hoping that you will find a second to take a look at it. Your opinion would 
> help me a lot.
> 

Hello Florin,

I'm starting a very busy work period right now, and I'll probably not
have the time to look into this in detail. I hope Patrick, David, or
someone else can check it out.

Christophe

-- 
Christophe Geuzaine

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