2D axisymmetric Inductive heating example with periodic recomputation of EM fields. This example uses the monolithic solver to solve the FEM/BEM system which allows to use higher timesteps (See choice of numls).
References :
[1]Approximate analytical solution for induction heating of solid cylinders
Todd A. Jankowski∗, Norma H. Pawley, Lindsey M. Gonzales, Craig A. Ross, James D. Jurney, Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM 87545, USA, Annual Mathematical Modelling 40 (2016) 2770-2782
*EM_2DAXI *EM_BOUNDARY *EM_CIRCUIT *EM_CONTROL *EM_CONTROL_TIMESTEP *EM_EOS_TABULATED1 *EM_MAT_001 *EM_ROTATION_AXIS *EM_SOLVER_BEM *EM_SOLVER_BEMMAT *EM_SOLVER_FEM *EM_SOLVER_FEMBEM_MONOLITHIC *END
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LST-Ansys)
$
$
$X
$X------------------------------------------------------------------------------
$# UNITS: (S.I)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
2D Axisymmetric Inductive heating example
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls
2 &numls
*EM_CONTROL_TIMESTEP
1,&em_macrdt
*EM_SOLVER_FEMBEM_MONOLITHIC
0,1,1.e-12,1.e-3
*EM_ROTATION_AXIS
0,0,0,0,0,1,32
*EM_2DAXI
100,12,,,62,82
*EM_2DAXI
101,1,,,51,71
*EM_2DAXI
102,2,,,52,72
*EM_2DAXI
103,3,,,53,73
*EM_2DAXI
104,4,,,54,74
*EM_2DAXI
105,5,,,55,75
*EM_2DAXI
106,6,,,56,76
*EM_2DAXI
107,7,,,57,77
*EM_2DAXI
108,8,,,58,78
*EM_2DAXI
109,9,,,59,79
*EM_2DAXI
110,10,,,60,80
*EM_2DAXI
111,11,,,61,81
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CIRCUIT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
1 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
1 51 71
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
2 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
2 52 72
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
3 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
3 53 73
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
4 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
4 54 74
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
5 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
5 55 75
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
6 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
6 56 76
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
7 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
7 57 77
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
8 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
8 58 78
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
9 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
9 59 79
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
10 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
10 60 80
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
11 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
11 61 81
*EM_BOUNDARY
21,9
*EM_BOUNDARY
22,9
*EM_BOUNDARY
23,9
*EM_BOUNDARY
24,9
*EM_BOUNDARY
25,9
*EM_BOUNDARY
26,9
*EM_BOUNDARY
27,9
*EM_BOUNDARY
28,9
*EM_BOUNDARY
29,9
*EM_BOUNDARY
30,9
*EM_BOUNDARY
31,9
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM MAT and EOS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
1 2 &em_cond1
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
2 4 &em_cond2 100
$ 2 4 &em_cond2
*EM_EOS_TABULATED1
100,201
$---+----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-6
*EM_SOLVER_BEMMAT
$# matid reltol
2 1e-6
*EM_SOLVER_BEM
$# reltol maxit stype precon uselas ncyclbem
1e-6 5000 2 2 1&em_bemmtx
*EM_SOLVER_FEM
$# reltol maxit stype precon uselas ncyclbem
1e-4 5000 1 1 1&em_femmtx
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM OUTPUT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_OUTPUT
$ matS matF solS solF mesh
2 2 2 2 0
*END
