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