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 nonzero 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: 20200803) 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.ikzberlin.de/s/Piej3YjiLxtGHKP
Best regards
Arved
Induction heating using MagnetoDynamics2DHarmonic solver [SOLVED]
Induction heating using MagnetoDynamics2DHarmonic solver [SOLVED]
 Attachments

 casemgdyn.sif
 (2.83 KiB) Downloaded 12 times

 casestatmag.sif
 (2.57 KiB) Downloaded 8 times
Re: Induction heating using MagnetoDynamics2DHarmonic solver
Dear ElmerTeam,
did you already have a look at this problem?
Can I do anything to help you with the issue? Do you need more information, another test case, ...?
Best regards
Arved
did you already have a look at this problem?
Can I do anything to help you with the issue? Do you need more information, another test case, ...?
Best regards
Arved

 Site Admin
 Posts: 3750
 Joined: 22 Aug 2009, 11:57
 Antispam: Yes
 Location: Espoo, Finland
 Contact:
Re: Induction heating using MagnetoDynamics2DHarmonic solver
Hi Arved,
I just had a look. It turns out that the new module did not find frequency where the old one had it. I guess when we included more different modules it turned out to be more natural to define it in Simulation section, for example. Now it looks for several places for it, but not the Solver section. Unfortunately there was no complaint so this could go unnoticed.
Attached are some modified .sif files that give pretty much the same results. I now use discontinuous fields for the newer solver as this is one of the things where it adds value. The idea is that the postprocessing is common to 2D/3D even though implementationwise they are rather different.
Attached is also .m file that makes the comparison.
I modified the code so that it will complain if the frequency is not found. There seems to be good agreement on the results. However, whether the sign of the vector potential could be discussed.
Thanx for nicely made test case!
Peter
I just had a look. It turns out that the new module did not find frequency where the old one had it. I guess when we included more different modules it turned out to be more natural to define it in Simulation section, for example. Now it looks for several places for it, but not the Solver section. Unfortunately there was no complaint so this could go unnoticed.
Attached are some modified .sif files that give pretty much the same results. I now use discontinuous fields for the newer solver as this is one of the things where it adds value. The idea is that the postprocessing is common to 2D/3D even though implementationwise they are rather different.
Attached is also .m file that makes the comparison.
I modified the code so that it will complain if the frequency is not found. There seems to be good agreement on the results. However, whether the sign of the vector potential could be discussed.
Thanx for nicely made test case!
Peter
 Attachments

 casestatmag.sif
 Modified sif with old module
 (2.85 KiB) Downloaded 10 times

 casemgdyn.sif
 Modified sif with new module
 (3.12 KiB) Downloaded 8 times

 compare.txt
 Matlab/octave file (change suffix to .m)
 (397 Bytes) Downloaded 9 times
Re: Induction heating using MagnetoDynamics2DHarmonic solver
Hi Peter,
thank you for your support! I also ran the test case and now get the following result:
For anyone who is interested in the theoretical background and the analytical solution: Have a look at Chapter 6 in S. Lupi: Fundamentals of Electroheat. Springer, 2017.
Best regards
Arved
thank you for your support! I also ran the test case and now get the following result:
For anyone who is interested in the theoretical background and the analytical solution: Have a look at Chapter 6 in S. Lupi: Fundamentals of Electroheat. Springer, 2017.
Best regards
Arved