[Getdp] issue with 2D axisymetric solenoid magnet resolution

[Frederic Trillaud]

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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Example.zip
Type: application/zip
Size: 348016 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130211/dd0c22c8/attachment.zip>