TEAM 6 problem
TEAM (Testing Electromagnetic Analysis Methods) represents an open international working group aiming to compare electromagnetic analysis computer codes. TEAM Workshops are meetings of this group. A series of TEAM Workshops was started in 1986 and has been organized in two-year rounds, each comprising a series of ”Regional” workshops and a ”Final” Workshop, as a satellite event of the COMPUMAG Conference. The TEAM problems consist in a list of test-problems, with precisely defined dimensions, constitutive laws of materials, excitation, etc., and each backed by a real laboratory device, on which measurements can be made.
Problem 6 [1] of the TEAM workshops for eddy current code comparison is a hollow conducting sphere in a uniform magnetic field (in the Z direction). The field is spatially uniform but varies with a given frequency. The objective is to measure the real and imaginary components of the magnetic flux in and outside the sphere.
This example is the equivalent of the T.E.A.M 11 problem solved in the frequency domain.
[1]RESULTS FOR A HOLLOW SPHERE IN UNIFORM FIELD (BENCHMARK PROBLEM 6). C.R.I. EMSON. COMPEL—The International Journal for Computation and Mathematics in Electrical and Electronic Engineering Vol. 7, No. 1/2, 89-101 © 1990 James & James Science Publishers Ltd
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_EXTERNAL_FIELD
*EM_DATABASE_POINTOUT
*EM_MAT_001
*EM_OUTPUT
*EM_POINT_SET
*EM_SOLVER_FEM
*EM_SOLVER_FEMBEM_MONOLITHIC
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$ emsol
4 &Freq &em_femmtx&em_bemmtx
*EM_CONTROL_TIMESTEP
$# tstype dtcons lcid
1 &em_dt
$---+----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,&B0
1000,&B0
$---+----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-9
*EM_SOLVER_BEMMAT
$# matid reltol
2 1e-9
*EM_SOLVER_BEMMAT
$# matid reltol
3 1e-9
*EM_SOLVER_FEMBEM_MONOLITHIC
,,1.e-16,1.e-5,500
*EM_SOLVER_BEM
$# reltol maxit stype precon uselas ncyclbem
1e-6 1000 2 4 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.,0.
2,40.,0.,0.
3,45.,0.,0.
4,47.,0.,0.
5,50.,0.,0.
6,51.,0.,0.
7,52.,0.,0.
8,53.,0.,0.
9,54.,0.,0.
10,55.,0.,0.
11,57.,0.,0.
12,60.,0.,0.
13,65.,0.,0.
14,70.,0.,0.
15,75.,0.,0.
16,80.,0.,0.
17,85.,0.,0.
18,90.,0.,0.
19,95.,0.,0.
20,100.,0.,0.
22,0.,0.,40.
23,0.,0.,45.
24,0.,0.,47.
25,0.,0.,50.
26,0.,0.,51.
27,0.,0.,52.
28,0.,0.,53.
29,0.,0.,54.
30,0.,0.,55.
31,0.,0.,57.
32,0.,0.,60.
33,0.,0.,65.
34,0.,0.,70.
35,0.,0.,75.
36,0.,0.,80.
37,0.,0.,85.
38,0.,0.,90.
39,0.,0.,95.
40,0.,0.,100.
*EM_OUTPUT
$ matS matF solS solF mesh
3 3 3 3
*END