[Getdp] Simple (?) magnetostatics

Jasper funkybob at gmail.com
Wed Feb 8 13:46:49 CET 2006


Hello,

I'm new to GetDP and I intend to use it for 3D magnostatic modelling. As a
start I took the demo/example from the wiki for 2D magnetostatics and put a
small square magnet in a square of air. This all worked OK as far as I can
tell.
My next attempt was to take this very simple example and extend it to three
dimensions, thus a small magnetic cube in a cube of air. Using examples from
the presentation/wiki/mailing list I tried to adapt the problem definition
file to three dimensions.
But with the 3D example GetDP does not converge to a solution. From what I
gather from the examples it might have something to do with a Gauge
condition, which I need in 3D but I wasn't able to find any links explaining
this...
I've attached the files for the 2D and 3D example. Any help would be very
much appreciated! Thanks in advance,

Jasper
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060208/9f65c384/attachment.html>
-------------- next part --------------
/* --------------------------------------------------------------------------
    This is a sample gmsh geometry file
   --------------------------------------------------------------------------*/

/* defining some parameters */

w = 1e-3;			// Width of square ( x )
h = w;				// Height of square ( y )
p = w / 5;			// Mesh size
r = 5e-3;			// Air ``radius''

/* defining the points */

Point(1) = {  w/2,  h/2,  0, p };
Point(2) = {  w/2, -h/2,  0, p };
Point(3) = { -w/2, -h/2,  0, p };
Point(4) = { -w/2,  h/2,  0, p };

Point(5) = {  r/2,  r/2,  0, p };
Point(6) = {  r/2, -r/2,  0, p };
Point(7) = { -r/2, -r/2,  0, p };
Point(8) = { -r/2,  r/2,  0, p };


/* defining the lines */

Line(1)  = { 1, 2 };
Line(2)  = { 2, 3 };
Line(3)  = { 3, 4 };
Line(4)  = { 4, 1 };

Line(5)  = { 5, 6 };
Line(6)  = { 6, 7 };
Line(7)  = { 7, 8 };
Line(8)  = { 8, 5 };

/* defining the surfaces */

Line Loop(10) = { 1, 2, 3 ,4 };
Plane Surface(1) = { 10 };
Line Loop(11) = { 5, 6, 7, 8 };
Plane Surface(2) = { 11, 10 };

/* defining the physical entities (for which elements will be saved) */
Physical Surface(100) = { 1 };					// Magnet
Physical Surface(101) = { 2 }; 					// Air
Physical Line(1000) = { 5, 6, 7, 8 } ;			// System boundaries


-------------- next part --------------
A non-text attachment was scrubbed...
Name: square.pro
Type: application/octet-stream
Size: 1881 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060208/9f65c384/attachment.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: square_ms.pro
Type: application/octet-stream
Size: 2516 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060208/9f65c384/attachment-0001.pro>
-------------- next part --------------
/* --------------------------------------------------------------------------
    This is a sample gmsh geometry file
   --------------------------------------------------------------------------*/

/* defining some parameters */

w = 1e-3;			// Width ( x-direction )
h = w;				// Height ( y-direction )
d = w;				// Depth ( z-direction )
p = w / 10;			// Characteristic length for meshing
r = 5e-3;			// ``Radius'' of air

/* defining the points */

Point(1) = {  w/2,  h/2,  d/2, p };
Point(2) = {  w/2, -h/2,  d/2, p };
Point(3) = {  w/2, -h/2, -d/2, p };
Point(4) = {  w/2,  h/2, -d/2, p };
Point(5) = { -w/2,  h/2,  d/2, p };
Point(6) = { -w/2, -h/2,  d/2, p };
Point(7) = { -w/2, -h/2, -d/2, p };
Point(8) = { -w/2,  h/2, -d/2, p };


/* defining the lines */

Line(1)  = { 1, 2 };
Line(2)  = { 2, 3 };
Line(3)  = { 3, 4 };
Line(4)  = { 4, 1 };
Line(5)  = { 5, 6 };
Line(6)  = { 6, 7 };
Line(7)  = { 7, 8 };
Line(8)  = { 8, 5 };
Line(9)  = { 1, 5 };
Line(10) = { 2, 6 };
Line(11) = { 3, 7 };
Line(12) = { 4, 8 };


/* defining the surfaces */

Line Loop(13) = { 1, 2, 3 ,4 };
Plane Surface(1) = { 13 };

Line Loop(14) = { 5, 6, 7, 8 };
Plane Surface(2) = { 14 };

Line Loop(15) = { 1, 10, -5, -9 };
Plane Surface(3) = { 15 };

Line Loop(16) = { -3, 11, 7, -12 };
Plane Surface(4) = { 16 };

Line Loop(17) = { 4, 9, -8, -12 };
Plane Surface(5) = { 17 };

Line Loop(18) = { -2, 10, 6, -11 };
Plane Surface(6) = {18};


/* defining the volumes */

Surface Loop(7) = { 1, 2, 3, 4, 5, 6 };
Volume(1) = { 7 };


/* define the air cube */

/* defining the points */

Point(101) = {  r/2,  r/2,  r/2, 4*p };
Point(102) = {  r/2, -r/2,  r/2, 4*p };
Point(103) = {  r/2, -r/2, -r/2, 4*p };
Point(104) = {  r/2,  r/2, -r/2, 4*p };
Point(105) = { -r/2,  r/2,  r/2, 4*p };
Point(106) = { -r/2, -r/2,  r/2, 4*p };
Point(107) = { -r/2, -r/2, -r/2, 4*p };
Point(108) = { -r/2,  r/2, -r/2, 4*p };


/* defining the lines */

Line(101) = { 101, 102 };
Line(102) = { 102, 103 };
Line(103) = { 103, 104 };
Line(104) = { 104, 101 };
Line(105) = { 105, 106 };
Line(106) = { 106, 107 };
Line(107) = { 107, 108 };
Line(108) = { 108, 105 };
Line(109) = { 101, 105 };
Line(110) = { 102, 106 };
Line(111) = { 103, 107 };
Line(112) = { 104, 108 };


/* defining the surfaces */

Line Loop(113) 	 	= { 101, 102, 103 ,104 };
Plane Surface(101)  = { 113 };
Line Loop(114) 	 	= { 105, 106, 107, 108 };
Plane Surface(102) 	= { 114 };
Line Loop(115) 	 	= { 101, 110, -105, -109 };
Plane Surface(103) 	= { 115 };
Line Loop(116) 	 	= { -103, 111, 107, -112 };
Plane Surface(104) 	= { 116 };
Line Loop(117) 	 	= { 104, 109, -108, -112 };
Plane Surface(105) 	= { 117 };
Line Loop(118) 	 	= { -102, 110, 106, -111 };
Plane Surface(106) 	= { 118 };


/* defining the volumes */

Surface Loop(107) = { 101, 102, 103, 104, 105, 106 };
Volume(2) = { 107, 7 };


/* defining the physical entities (for which elements will be saved) */

Physical Volume(101) = { 1 }; 			// Cube
Physical Volume(102) = { 2 }; 			// Air

Physical Surface(1000) = { 101, 102, 103, 104, 105, 106 };		// Boundary
Physical Surface(1001) = { 1, 2, 3, 4, 5, 6 };					// Gauge
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cube.pro
Type: application/octet-stream
Size: 1568 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060208/9f65c384/attachment-0002.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cube_ms.pro
Type: application/octet-stream
Size: 2645 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060208/9f65c384/attachment-0003.pro>