I'm performing a 2D cylindrical coordinates simulation of a metal detector using MagnetoDynamics2DHarmonic solver.
The simulation is basically working but there is a strange anomaly. bSolver is giving me non-zero currents in regions where the
conductivity is zero. According to the models document the conduction current should be given by the equation
Jc = -j*sigma*omega*Az.
According to this equation it would be impossible for the current to be nonzero in any region with a conductivity (sigma) of zero.
Simulation files and some paraview screenshots scan be found at the link below.
https://1drv.ms/u/s!AkpB2-DyzyN6hro3esI ... A?e=QfQe1o
The metal detector coils are surrounded by a metal case where there are very high induced current densities. The case currents are
2000 times bigger than the spurious currents in air regions of zero conductivity.
Issue with currents computed by bSolver and MagnetoDynamics2DHarmonic
Issue with currents computed by bSolver and MagnetoDynamics2DHarmonic
- Attachments
-
- case.sif
- (2.88 KiB) Downloaded 232 times
-
- Site Admin
- Posts: 4827
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Issue with currents computed by bSolver and MagnetoDynamics2DHarmonic
Hi
The problem is that whereas Az is continuous sigma*Az is not. So what you are seeing is the effort to map discontinuous field to continuous basis functions. The remedy is to use discontinuous Galerkin (DG) basis as may be done with MagnetoDynamicsCalcFields (instead of BSolver).
-Peter
The problem is that whereas Az is continuous sigma*Az is not. So what you are seeing is the effort to map discontinuous field to continuous basis functions. The remedy is to use discontinuous Galerkin (DG) basis as may be done with MagnetoDynamicsCalcFields (instead of BSolver).
-Peter
Re: Issue with currents computed by bSolver and MagnetoDynamics2DHarmonic
I assumed that Elmer would model the conductivity as continuous. To deal with this I had added a very thin zero conductivity layer next to the high conductivity case region. (See attached figure). I'm seeing non-zero currents at nodes far away from the high conductivity regions (nodes where all adjacent nodes have zero conductivity). Shouldn't this result in the conductivity being zero in all triangles surrounding the node. How can BSolver still compute a non-zero current for such a node?
- Attachments
-
- Spurious currents many cells away from non-zero conductivity region
- spurios_current_mid_range_zoom.png (295.76 KiB) Viewed 3258 times
-
- Site Admin
- Posts: 4827
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Issue with currents computed by bSolver and MagnetoDynamics2DHarmonic
Hi, Its difficult to comment without the whole case. However, it is clear that if you you have a jump in currents over elements this cannot be presented by continuous function. This is difficult to ovecome without a discontinuous visualization space. -Peter