cyl.axis

[Christophe Geuzaine]

Samuel Kvasnica wrote:
> 
> Christophe Geuzaine wrote:
> 
> > Yes. It seems that there is a small problem with 2nd order interpolation
> > in a few elements in the mesh (near the symmetry axis). I'll investigate
> > this when I have some free time. By now, you should stick with 1st order
> > interpolation (and maybe apply some smoothing on the result if
> > smoothness is an issue).
> 
> ok, thank you very much for your reply. Anyway, I still have a problem how to
> correctly set a
> constraint on the z-axis. Maybe I'm wrong, but I think that only the r-component of B
> should be zero
> there. Any hint ?
> 

For symmetry, you need n . b = 0 on the vertical axis. Since b = curl a,
this means n x a = 0. You should thus have all circulations of 'a' along
the perpendicular edges to be zero. And this is what you do with 

 Name MagneticVectorPotential_2D ;
    Case {
      { Region SurfaceGh0  ; Value 0. ; }
    }
  }

My best,

-- 
Christophe Geuzaine

Tel: +32-(0)4-366.37.10    mailto:Christophe.Geuzaine at ulg.ac.be
Fax: +32-(0)4-366.29.10    http://www.geuz.org