2D Axisymmetric Inductive heating — Welcome to LS-DYNA Examples

Home » Knowledge Base » EM / Electromagnetics » Tutorial Eddy Currents » 2D Axisymmetric Inductive heating

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).

Description

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

Animated Result

Temperature fringes

Watch animation on YouTube

Keywords

*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

Reduced Input

$-----------------------------------------------------------------------------
$
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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      &amp  
$ 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

Figures

Cylinder Flow

Download