I think I've found a bug in the whitneyharmonicsolver.
The current density is wrong in the harmonic solution.
Magnetic field with eddy currents 2d/3d
-
- Posts: 27
- Joined: 13 Aug 2013, 16:50
- Antispam: Yes
Re: Magnetic field with eddy currents 2d/3d
- Attachments
-
- input.zip
- (412.25 KiB) Downloaded 419 times
-
- whitney_current_density.png (51.27 KiB) Viewed 4861 times
-
- harmonic_current_density.png (28.94 KiB) Viewed 4861 times
Re: Magnetic field with eddy currents 2d/3d
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?
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?
Re: Magnetic field with eddy currents 2d/3d
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
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
-
- Posts: 27
- Joined: 13 Aug 2013, 16:50
- Antispam: Yes
Re: Magnetic field with eddy currents 2d/3d
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
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
-
- Posts: 27
- Joined: 13 Aug 2013, 16:50
- Antispam: Yes
Re: Magnetic field with eddy currents 2d/3d
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
In Elmer I tried
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.
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
Code: Select all
Potential = real 0
Code: Select all
P re {e} = real 0
P im {e} = real 0
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.
-
- Site Admin
- Posts: 4843
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Magnetic field with eddy currents 2d/3d
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.akirahinoshiro wrote: Could you explain why your definitions in linear system worked and mine not?
-Peter
Re: Magnetic field with eddy currents 2d/3d
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
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