[Getdp] issue with 2D axisymetric solenoid magnet resolution

[Frederic Trillaud]

Hi Ruth,

I think that I found my answer. Here is a solution using VolAxiSphVolume
and VolAxi instead of VolSphVolume and Vol in the Jacobian method. The
picture in my previous email was obviously the module of B and not the
potential vector! (my bad). Here is a new picture with, I think, a
relatively good solution.

Frederic

On Tue, 2013-02-12 at 11:01 +0100, Ruth V. Sabariego wrote:
> Hi Frederic,
> The problem came from your geometry: the physical line for the
> symmetry and infinity where you define the dirichlet boundary
> condition was not properly defined.
> You have to use the geometrical lines and not the loops.
> See the small correction in the geo file. 
> 
> 
> Best regards,
> Ruth
> 
> 
> 
> 
> 
> On 11 Feb 2013, at 21:05, Frederic Trillaud
> <ftrillaudp at pumas.ii.unam.mx> wrote:
> 
> > Dear all,
> > 
> > I have an issue with the solution of a 2D axi-symmetric problem
> > consisting of a solenoid magnet. The magnet field is not as
> > expected!
> > 
> > Here are the files to run the model in the example.zip attached
> > folder.
> > 
> > solenoid2D.par -> parameters
> > solenoid2D.geo -> gmsh geometry
> > main.pro -> model in getdp
> > magnetostaticProblem.pro -> formulation
> > etc...
> > 
> > gmsh -2 solenoid2D.geo -msh solenoid2D.msh
> > getdp main -msh solenoid2D.msh -solve potentialVectorResolution
> > getdp main -msh solenoid2D.msh -pos fieldMap
> > 
> > Dirichlet is applied on the axis of symmetry and the on the external
> > shell. A fixed current density is given normal to the cross-section
> > plan.
> > 
> > I have tried with different Jacobian volume, VolAxi and
> > VolAxiSphVolume
> > but without success in addition to what is included in the given
> > test.
> > 
> > Best regards,
> > 
> > Frederic
> > <Example.zip>_______________________________________________
> > getdp mailing list
> > getdp at geuz.org
> > http://www.geuz.org/mailman/listinfo/getdp
> 
> —
> Dr. Ir. Ruth V. Sabariego
> University of Liege, Electrical Engineering and Computer Science
> http://ace.montefiore.ulg.ac.be/
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: moduleB.png
Type: image/png
Size: 141309 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130212/54b96b0f/attachment.png>