[Getdp] issue with 2D axisymetric solenoid magnet resolution

Frederic Trillaud ftrillaudp at pumas.ii.unam.mx
Tue Feb 12 15:34:59 CET 2013


Hi Ruth,

One additional question, I am still surprised by the results on the
infinite boundary as the field has some non-negligible values away from
the magnet, do you know where this comes from?, see picture attached.

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: potentialVector.png
Type: image/png
Size: 125491 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130212/5ecf8faa/attachment.png>