TEAM 12 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, excitations, etc., and each backed by a real laboratory device, on which measurements can be made.
The TEAM 12 problem is a coupled problem with moving conductor. A clamped beam is placed in a uniform magnetic field. The magnetic field will have a first component exponentially decaying with time that will generate an induced current in the beam that will in turn, interact with the second constant component of the field and create a Lorentz force that will cause the beam’s movement. The motion of the beam causes the current and deflection to be very different from what they would be if coupling were not present. The experimental results are based on a FELIX experiment [1].
In addition to the traditional T.E.A.M 12 problem, a modified version is also provided where the source of the magnetic field and force applied on the beam is induced by a moving magnet located at the end of the beam.
[1] T. Nakata and K. Fujiwara. Summary of results for benchmark
problem 13 (3-d nonlinear magnetostatic model). COMPEL-The international journal for computation and mathematics in electrical and electronic engineering, 11(3):345–369, 1992.
*EM_SOLVER_BEMMAT
*EM_CONTROL
*EM_CONTROL_COUPLING
*EM_CONTROL_TIMESTEP
*EM_EXTERNAL_FIELD
*EM_MAT_001
*EM_MAT_002
*EM_OUTPUT
*EM_PERMANENT_MAGNET
*EM_POINT_SET
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_FEMBEM_MONOLITHIC
*KEYWORD
$
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
1 &em_femmtx&em_bemmtx
*EM_CONTROL_TIMESTEP
1,&emdt
*EM_CONTROL_COUPLING
,2
*EM_SOLVER_FEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ reltol maxiter stype precon uselast ncyclfem msglvlFem
1.e-4 10000 1 1 1 0
*EM_SOLVER_FEMBEM_MONOLITHIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$ COUPLING KSP_TYPE KSP_ATOL KSP_RTOL KSP_MAXIT GMRES_RST FEM_PC DEBUG
0 1.e-9 1.e-4 100000 1
$--------1---------2---------3---------4---------5---------6---------7---------8
$ FORCE NEWT_STOL NEWT_ATOL NEWT_RTOL NEWT_MAXI
$ 1.e-9 1.e-9 1.e-5
$--------1---------2---------3---------4---------5---------6---------7---------8
$ LS_ON LS_FTOL LS_GTOL LS_RTOL LS_SAMPL LS_NUM_SP LS_MAXFUN
$ 1 1.e-9 0.01 1.e-5 30000
$--------1---------2---------3---------4---------5---------6---------7--------
*EM_SOLVER_BEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ reltol maxiter solvetype precon uselast ncyclbem
1.e-8 10000 2 2 1
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
1 1.e-12
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
2 1.e-12
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
3 1.e-12
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$ ID Part ID - North South M=0.38/mu0
3 3 4 5 &mag
*EM_MAT_002
$# mid mtype sigma eosid
1 4 &sigma1 0
*EM_MAT_002
$ Magnetic material 2
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid murel eosmu
3 4 &sigma3 0 &murel1
*EM_OUTPUT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matS matF solS solF mesh memory timing d3plotAsc
2 2 2 2 0 0 0
*END