Magnetic field with eddy currents 2d/3d

Numerical methods and mathematical models of Elmer
akirahinoshiro
Posts: 27
Joined: 13 Aug 2013, 16:50
Antispam: Yes

Re: Magnetic field with eddy currents 2d/3d

Post by akirahinoshiro »

I think I've found a bug in the whitneyharmonicsolver.
The current density is wrong in the harmonic solution.
Attachments
input.zip
(412.25 KiB) Downloaded 419 times
whitney_current_density.png
whitney_current_density.png (51.27 KiB) Viewed 4861 times
harmonic_current_density.png
harmonic_current_density.png (28.94 KiB) Viewed 4861 times
NickR7
Posts: 18
Joined: 29 Jul 2013, 18:03
Antispam: Yes

Re: Magnetic field with eddy currents 2d/3d

Post by NickR7 »

I've seen that sort of thing before and I think it was because I wasn't writing the vector fields to file correctly.

For one thing, I should think that your current density will be a complex quantity. So I would think you would want to write the imaginary and real parts of the current density to file separately. I doubt paraview can handle complex vector quantity (or maybe you need to do something specific to get it to recognize complex vectors and plot it as a time evolution or something).

Look and see how a test file for WhitneyAVHarmonic writes output to .vtu format (if there is an example of that).

Also, you still have

Simulation Type = Transient

Is this correct for the harmonic solver?
millim
Posts: 39
Joined: 04 Apr 2012, 14:58
Antispam: Yes

Re: Magnetic field with eddy currents 2d/3d

Post by millim »

Hi
How paraview can handle complex Vectors See:http://www.flickr.com/photos/ajvanne/79 ... otostream/

For transient you have to use something like cos(wt) and static,
for stady just the angular frequncy and harmonic

Br
akirahinoshiro
Posts: 27
Joined: 13 Aug 2013, 16:50
Antispam: Yes

Re: Magnetic field with eddy currents 2d/3d

Post by akirahinoshiro »

Now I tried using statcurrentsolver.
The solver runs fine but after adding the whitney parts the solver do not converge.
What did I wrong for the current?

I uploaded the unv file here:
http://www.file-upload.net/download-809 ... e.unv.html
Attachments
case.sif
(3.5 KiB) Downloaded 390 times
akirahinoshiro
Posts: 27
Joined: 13 Aug 2013, 16:50
Antispam: Yes

Re: Magnetic field with eddy currents 2d/3d

Post by akirahinoshiro »

I got it to work but only in full 3d and 2d.
If I cut the 3d model on the symmetry line I need to define a new boundary condition on the front face.
My problem is I do not know how to do that in Elmer.
In Comsol I set

Code: Select all

electric potential V0=0V
In Elmer I tried

Code: Select all

Potential = real 0

Code: Select all

P re {e} = real 0
P im {e} = real 0
and some others.
But the results won't converge or produces errors.
Is there anyone who knows what I did wrong?
The projects are uploaded here:
eddy_alloy.zip

Thank you again.
raback
Site Admin
Posts: 4843
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Magnetic field with eddy currents 2d/3d

Post by raback »

akirahinoshiro wrote: Could you explain why your definitions in linear system worked and mine not?
The solution of the potential is indefinite having a nontrivial null-space. Hence you cannot use an ILU preconditioner as indefinite problems don't have an inverse.

-Peter
jvelasco
Posts: 10
Joined: 03 Feb 2021, 16:45
Antispam: Yes

Re: Magnetic field with eddy currents 2d/3d

Post by jvelasco »

Hi guys,

This old thread contains TEAM Problem 7 discussion. Although it's been a while, I will leave this link here just in case anyone is interested in modelling eddy currents in 3D. The setup is quite straight forward and hopefully useful to the community.

https://github.com/ElmerCSC/elmer-elmag/tree/main/TEAM7

Best Regards,

Jonathan
Post Reply