Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Skin effect error

Please login with a confirmed email address before reporting spam

Hello all,

I have problem in simulating induced current in the cross section of a wire which is copper. The problem is that based on the skin effect theory the current must be in one direction in the outer cross section of the wire and when we go to inner cross section, its density must fall, but my simulations show that the current is in one direction in the outer cross section and in the opposite direction in the inner cross section. I mean when we go from outer cross section to inner cross section the direction of the current changes alternatively.

I have attached the 2D model which I have created to explain the issue. In this model I have exited just one infinite wire of copper by external current density and 13MHz frequency but when simulate this wire I get the alternating current from outer cross section to inner cross section.

I also see the same problem in simulating two loop problem in which I have excited the first loop by external current density boundary condition and wanted to see current density of the second loop.

Thank you


3 Replies Last Post Dec 31, 2011, 10:48 a.m. EST
Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 31, 2011, 6:13 a.m. EST
Hi

are you sure you model is complete ? you are in 2D (not 2D axi) solving for only out of plane vector component, so you have a infinite long wire of copper, no air/domain around so your field is not looping even in any vacuum outside of your copper. And as you have no "return in curent".

For me with that you can only get a cylindrical solution, as you see it, no ?

--
Good luck
Ivar
Hi are you sure you model is complete ? you are in 2D (not 2D axi) solving for only out of plane vector component, so you have a infinite long wire of copper, no air/domain around so your field is not looping even in any vacuum outside of your copper. And as you have no "return in curent". For me with that you can only get a cylindrical solution, as you see it, no ? -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 31, 2011, 9:30 a.m. EST
I really appreciate your time and consideration.

With previous model I just wanted to say that even in this simple example (just one infinite), I do not get the right answer. What can I do to make it right?

I have attached the original model in which I have forced the upper loop to have current by external current density boundary condition and set its material to air to avoid this problem (alternating current from outer cross section to inner cross section), but in the bottom loop in which I wanted to see the induced current and its material is set to copper, I see the alternating current as I said before. For your information I use magnetic fields interface, make the entire mesh extremely fine and make the bottom loop mesh too more fine than extremely fine. I also added boundary layer mesh (I do not know if this is right to add boundary layer mesh) but the result did not change.

Thank you.
I really appreciate your time and consideration. With previous model I just wanted to say that even in this simple example (just one infinite), I do not get the right answer. What can I do to make it right? I have attached the original model in which I have forced the upper loop to have current by external current density boundary condition and set its material to air to avoid this problem (alternating current from outer cross section to inner cross section), but in the bottom loop in which I wanted to see the induced current and its material is set to copper, I see the alternating current as I said before. For your information I use magnetic fields interface, make the entire mesh extremely fine and make the bottom loop mesh too more fine than extremely fine. I also added boundary layer mesh (I do not know if this is right to add boundary layer mesh) but the result did not change. Thank you.


Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 31, 2011, 10:48 a.m. EST
Hi

I have some issues with your 2nd model too, for me Cu should be domains #2+3, and air domain #1
the curent should be defined in #2+3 (perhaps with opposite signs)

The mesh should have a boundary element for the Cu domains and the sub node boundary mesh for the 2x4 segments of the two Cu domain boundaries (the 1/4 circle segments), then you can relax somewhat your extremely fine mesh elsewhere.

Still, once changed, I agree the result looks strange, but I'm afraid that could come from your old version. The v4.0 was not nicely ironed out, you should use a more recent version or the stable 3.5a, there are probably even several bugs left in the early 4.0


--
Good luck
Ivar
Hi I have some issues with your 2nd model too, for me Cu should be domains #2+3, and air domain #1 the curent should be defined in #2+3 (perhaps with opposite signs) The mesh should have a boundary element for the Cu domains and the sub node boundary mesh for the 2x4 segments of the two Cu domain boundaries (the 1/4 circle segments), then you can relax somewhat your extremely fine mesh elsewhere. Still, once changed, I agree the result looks strange, but I'm afraid that could come from your old version. The v4.0 was not nicely ironed out, you should use a more recent version or the stable 3.5a, there are probably even several bugs left in the early 4.0 -- Good luck Ivar

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.