TEAM 11 problem

Problem 11 of the TEAM workshops for eddy current code comparison is a hollow conducting sphere in a uniform magnetic field. The applied field varies with time as a step function, requiring a fully transient solution. Magnetic field in and outside the sphere at different instants is measured and compared to the analytical solution.

*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_EXTERNAL_FIELD
*EM_DATABASE_POINTOUT
*EM_MAT_001
*EM_OUTPUT
*EM_POINT_SET
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_SOLVER_FEMBEM

*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$ emsol
1 &em_femmtx&em_bemmtx
*EM_CONTROL_TIMESTEP
$# tstype dtcons lcid
1 &em_dt
$*EM_BOUNDARY
$9,1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM EXTERNAL FIELD $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_EXTERNAL_FIELD
1,1,1,300,300,200
*DEFINE_CURVE
300
0,0.
1000,0.
*DEFINE_CURVE
200
0,0.
1.e-5,1.
1000,1.
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM MAT and EOS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$ em_mid mtype sigma eos
2 2 &em_cond
*EM_MAT_001
$ em_mid mtype sigma eos
1 1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM SOLVER $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_SOLVER_BEMMAT
$# matid reltol
1 1e-8
*EM_SOLVER_BEMMAT
$# matid reltol
2 1e-8
*EM_SOLVER_BEM
$# reltol maxit stype precon uselas ncyclbem
1e-6 1000 2 2 1
*EM_SOLVER_FEM
$# reltol maxit stype precon uselas ncyclbem
1e-6 1000 1 1 1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM OUTPUT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_DATABASE_POINTOUT
1
1
*EM_POINT_SET
1
1,0.,0.,50.0
2,0.,0.,51.0
3,0.,0.,52.0
4,0.,0.,53.0
5,0.,0.,54.0
6,0.,0.,55.0
7,0.,0.,56.0
8,0.,0.,57.0
9,0.,0.,58.0
10,0.,0.,59.0
11,0.,0.,60.0
12,0.,0.,61.0
13,0.,0.,62.0
14,0.,0.,63.0
15,0.,0.,64.0
16,0.,0.,65.0
17,0.,0.,66.0
18,0.,0.,67.0
19,0.,0.,68.0
20,0.,0.,69.0
21,0.,0.,70.0
22,0.,0.,71.0
23,0.,0.,72.0
24,0.,0.,73.0
25,0.,0.,74.0
26,0.,0.,75.0
27,0.,0.,76.0
28,0.,0.,77.0
29,0.,0.,78.0
30,0.,0.,79.0
31,0.,0.,80.0
32,0.,0.,81.0
33,0.,0.,82.0
34,0.,0.,83.0
35,0.,0.,84.0
36,0.,0.,85.0
37,0.,0.,86.0
38,0.,0.,87.0
39,0.,0.,88.0
40,0.,0.,89.0
41,0.,0.,90.0
42,0.,0.,91.0
43,0.,0.,92.0
44,0.,0.,90.0
45,0.,0.,91.0
46,0.,0.,92.0
47,0.,0.,93.0
48,0.,0.,94.0
49,0.,0.,95.0
50,0.,0.,96.0
51,0.,0.,97.0
52,0.,0.,98.0
53,0.,0.,99.0
54,0.,0.,100.0
55,0.,0.,101.0
101,50,0.,0.
102,51,0.,0.
103,52,0.,0.
104,53,0.,0.
105,54,0.,0.
106,55,0.,0.
107,56,0.,0.
108,57,0.,0.
109,58,0.,0.
110,59,0.,0.
111,60,0.,0.
112,61,0.,0.
113,62,0.,0.
114,63,0.,0.
115,64,0.,0.
116,65,0.,0.
117,66,0.,0.
118,67,0.,0.
119,68,0.,0.
120,69,0.,0.
121,70,0.,0.
122,71,0.,0.
123,72,0.,0.
124,73,0.,0.
125,74,0.,0.
126,75,0.,0.
127,76,0.,0.
128,77,0.,0.
129,78,0.,0.
130,79,0.,0.
131,80,0.,0.
132,81,0.,0.
133,82,0.,0.
134,83,0.,0.
135,84,0.,0.
136,85,0.,0.
137,86,0.,0.
138,87,0.,0.
139,88,0.,0.
140,89,0.,0.
141,90,0.,0.
142,91,0.,0.
143,92,0.,0.
144,90,0.,0.
145,91,0.,0.
146,92,0.,0.
147,93,0.,0.
148,94,0.,0.
149,95,0.,0.
150,96,0.,0.
151,97,0.,0.
152,98,0.,0.
153,99,0.,0.
154,100,0.,0.
155,101,0.,0.
200,52.50,0,0.0
201,52.30,0,4.576
202,51.70,0,9.116
203,50.71,0,13.59
204,49.33,0,17.96
205,47.58,0,22.19
206,45.67,0,26.25
207,43.00,0,30.11
208,40.22,0,33.75
209,37.12,0,37.12
210,33.75,0,40.22
211,30.11,0,43.00
212,26.25,0,45.67
213,22.19,0,47.58
214,17.96,0,49.33
215,13.59,0,50.71
216,9.116,0,51.70
218,4.576,0,52.30
219,0.000,0,52.50
*EM_OUTPUT
$ matS matF solS solF mesh
2 2 2 2
*END