Hi all!
First of all I'd like to thank this forum, here I have found many answers I needed before. But now I couldn't find any at all.
Thing is, I have a quite complicated 2D model of a superconducting magnet. My model is a transversal section of the magnet where we can see all the conductors as "domains" and the electric current is perpendicular to the plane.
The electrical resistivity of the conductor is a function of the temperature and magnetic field. So I imported a table and COMSOL interpolates it to obtain the resistivity, that works fine.
Now I need to implement an electric current that follows
Where I0 and L are constants.
R is the resistivity (rho) integrated all over the conductor's domain, it reads:
I used "intop" to tackle this double integral over the conductor's domain, and the outer one is the length in the "z" direction.
COMSOL apparently does simulate with this setup, but the calculations are too too slow. It takes a lot of memory as well (plus than 80GB!) However, if I use I(t)=I0, it works nice and fast, so I'm sure there's a problem with the "intop" function.
So I ask you if there's something I am missing, or some suggestions to make the simulation faster.
Thanks in advance!
Lucas