Induction heating using MagnetoDynamics2DHarmonic solver [SOLVED]
Posted: 01 Sep 2020, 15:02
Dear Elmer community,
I did some simulations for verification of axisymmetric induction heating and found that the results of the (obsolete) StatMagSolver agree well with the expected analytical solution, whereas the induced heat computed with the newer MagnetoDynamics2DHarmonic solver together with MagnetoDynamicsCalcFields or BSolver shows big deviations.
The investigated case is a cylinder in a long coil, where I analyze the total heat induced in the central part of the cylinder. The analytical value for the total induced heat is 193.24 W. The curves of the induced heat over the radius of the cylinder are plotted here: Furthermore, regarding the mgdyn solver, I made the following observations:
- The imaginary part of the magnetic flux density / field strength is zero in the whole domain
- The real part of the magnetic flux density / field strength is constant in the whole cylinder
- Changing the mesh to quadratic elements (by recombination in gmsh) produces a non-zero imaginary magnetic flux density / field strength, even though it is very low and looks weird
- MagnetoDynamicsCalcFields or BSolver produce the output: SolveSystem: Solution trivially zero!
I am using the solver Version: 8.4 (Rev: Release, Compiled: 2020-08-03) downloaded from the link provided on this website.
Attached you can find the the sif files I used for my simulations, the mesh files (in two different resolutions) are uploaded here: https://cloud9.ikz-berlin.de/s/Piej3YjiLxtGHKP
Best regards
Arved
I did some simulations for verification of axisymmetric induction heating and found that the results of the (obsolete) StatMagSolver agree well with the expected analytical solution, whereas the induced heat computed with the newer MagnetoDynamics2DHarmonic solver together with MagnetoDynamicsCalcFields or BSolver shows big deviations.
The investigated case is a cylinder in a long coil, where I analyze the total heat induced in the central part of the cylinder. The analytical value for the total induced heat is 193.24 W. The curves of the induced heat over the radius of the cylinder are plotted here: Furthermore, regarding the mgdyn solver, I made the following observations:
- The imaginary part of the magnetic flux density / field strength is zero in the whole domain
- The real part of the magnetic flux density / field strength is constant in the whole cylinder
- Changing the mesh to quadratic elements (by recombination in gmsh) produces a non-zero imaginary magnetic flux density / field strength, even though it is very low and looks weird
- MagnetoDynamicsCalcFields or BSolver produce the output: SolveSystem: Solution trivially zero!
I am using the solver Version: 8.4 (Rev: Release, Compiled: 2020-08-03) downloaded from the link provided on this website.
Attached you can find the the sif files I used for my simulations, the mesh files (in two different resolutions) are uploaded here: https://cloud9.ikz-berlin.de/s/Piej3YjiLxtGHKP
Best regards
Arved