Elmer Discussion Forum

Hi,

Attached is a figure of the nodal heating. Many fields look quite similar to the eye. The imaginary part of the vector potential and this though are clearly different. Here you see the heating with just the IBC condition on the left and with the zero-current constraint on the right.

-Peter

Hi Roland, I added a new type of linear constraint to the library part of the code. This is applicable to Robin type of BCs where we want to constrain to total flux. For Neumann kind of BCs controlling of total flux is easy since we know what we are putting in. For Robin type of conditions it is mor...

Hi, There was just a wrong sign. Now things look good. Attached is comparison between the two models. CalcFields also works. Ask for "Calculate Nodal Heating = True" when using IBC. That can be used as "Temperature Load = Nodal Joule Heating" when solving for temperature. This wa...

Hi, Have you varied the timestep? I would guess you need easily ~20 steps for full cycle to get reasonable results. How does the damping depend on the this? The damping should be eliminated when timestep goes to zero. Time integration schemes tend to have some damping. In fact it is difficult to hav...

Hi,

You need to recompile / fetch new nightly build. This is in the library.

This is the commit:
https://github.com/ElmerCSC/elmerfem/co ... 3c42128635

-Peter

Hi I think there is something fishy with the mesh. The "surface with edges" on Paraview shows like there would be double set of edges on body 1. I agree with Kevin that ~5 orders of magnitude is challenging aspect ratio but the convergence is lost a lot earlier. I would expect somewhat mor...

Hi I had now some time to look at it closer. Indeed, there was a bug! That should be fixed now. Thanx for reporting! The bug manifested itself when there was different elementtypes within the same material such that the one with smaller number of gaussian integration points was first. If also the va...

Hi Roland,

Really nice test case!

The error came in the postprocessing part. It cannot yet be considered. I added skipping of the not implemented part. There seems to be some effect particularly on "Potential Im". Maybe you can compare to reference.

-Peter

Hi, maybe you can provide a full test case where this happens. -Peter

Hi Roland, Yes, a version that has been compiled after yesterday. Then add these values at some BC: Layer Electric Conductivity = Real Layer Relative Permeability = Real You should see the effect in the solution of A_phi with and without them. The Joule heating for the boundary is not yet included. ...