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, excitations, etc., and each backed by a real laboratory device, on which measurements can be made.
The TEAM 1 problem consists of a conducting cylinder of finite length [1] placed in a uniform magnetic field. The magnetic field decays exponentially with time. The main objective of the problem is to calculate the total circulating current in the cylinder at various locations (Z=0 m and Z = 0.5 m) as well as the Total Power Loss.
[1] K. Davey, THE FELIX CYLINDER PROBLEM (INTERNATIONAL EDDY CURRENT WORKSHOP PROBLEM 1) , COMPEL, Vol 7. Num 1& 2 11-27 (1988).
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_EXTERNAL_FIELD
*EM_CIRCUIT_ROGO
*EM_DATABASE_PARTDATA
*EM_MAT_001
*EM_OUTPUT
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_SOLVER_FEMBEM/*EM_SOLVER_FEMBEM_MONOLITHIC
*KEYWORD
*TITLE
TEAM 1a
*PARAMETER
R Tend 0.24
R dtplot 2.e-3
R dt 1.e-6
$
$--- em
R em_cond 25.38e6
R em_dt 2.e-4
R B0 0.1
R Tshift 0.12
R Tau 0.0397
*DEFINE_CURVE_FUNCTION
20
DELAY(LC200,&Tshift, &B0)
*DEFINE_CURVE_FUNCTION
200
&B0*EXP(-(TIME)/&Tau)
*EM_EXTERNAL_FIELD
$--------1---------2---------3---------4---------5---------6---------7---------8
$ fieldId fieldType FieldDef lcidX lcidY lcidZ
1 1 1 21 20 21
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls emdtinit emdtmax emtinit emtend ncyclFem ncyclBem
1
*EM_CONTROL_COUPLING
,2
*EM_CONTROL_TIMESTEP
1,&em_dt
$ uncomment those lines to turn on monolithic solver
$ timestep (em_dt) can be risen for example to 1.e-3
$*EM_SOLVER_FEMBEM_MONOLITHIC
$,,1.e-16,1.e-4
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eos skinDepth
1 4 &em_cond
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eos skinDepth
2 4 &em_cond
*EM_OUTPUT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matS matF solS solF mesh
2 2 2 2
*EM_CIRCUIT_ROGO
1,1,1,1
*EM_CIRCUIT_ROGO
2,2,1,1
*EM_CIRCUIT_ROGO
3,3,1,1
*EM_CIRCUIT_ROGO
4,4,1,1
*EM_DATABASE_PARTDATA
1
*MAT_RIGID
1 7000.6.8900E+10 0.330000 0.000 0.000 0
1,7,7
*MAT_ELASTIC
2,7000,6.89e10,0.33
*CONTROL_HOURGLASS
1 0.100000
$ Warning : normally not recommended to prescribe dt like
$ that due to mass scaling but here we only care about EM solution
*CONTROL_TIMESTEP
,,,,&dt,16
*DEFINE_CURVE
16
0,1.e-4
1000,1.e-4
*CONTROL_TERMINATION
&Tend
*SECTION_SOLID
1 1 0
*DATABASE_BINARY_D3PLOT
&dtplot
*DEFINE_CURVE
21
0,0
0.1,0
*PART
1 1 1
*PART
2 1 2
*INCLUDE
mesh_team1.k
*END
