[MgDynHarmonic,SaveScaler]current in body force

[xiejinzhi]

Hello，
Recently i build a 2D submarine cable model in gmsh and use MgDynHarmonic to calculate the induced current and joule heating.The body force are the 3 cores inside the cable,phase A,phsse B and phsse C.In phase A,B and C,the current is 982+j0,-491+j850,-491-j850（Phase angle difference of 120 degrees）.The current density in Paraview shows well.



But when i use SaveScalars to calculate the current density body integral in 3 cores, i got some problems:
1. the solver SaveScalars could only calculate the Real part of current(see the row of data with 0 Angular frequency ), but in sif file,it is "Variable = current density e".If i change it to "Variable = current density re e" or "Variable = current density im e",the result given is 0.
2.the current in the body force is not fixed value,it changes with Angular frequency(from 0 to 314.16).I think it should always be fixed value because it is force.

Any help is appreciated,thank you.

[xiejinzhi]

the sif file

[zeljkoT]

Hi xiejinzhi

I just came across your post because I'm dealing with a similar question.
Does the B vector field look reasonable?

Cheers

[xiejinzhi]

Hi xiejinzhi

I just came across your post because I'm dealing with a similar question.
Does the B vector field look reasonable?

Cheers

Hi
Use "circuits_harmonic" as power source(see the case in "elmerfem-devel\fem\tests\circuits2D_harmonic_massive")，not "body force",the problem will be solved.
Cheers

[zeljkoT]

Hi

Yes, indeed.
Works.
Thank you and cheers