Elmer Discussion Forum

Hi The "ViewFactors" binary is called by "ElmerSolver" automatically. However, if the problem is the view factors then just compute those. If the view factors exist and the nodes are the same they are not recomputed unless requested. If you decleare your boundaries open then the ...

This piece of debug code in ViewFactors.F90 let's this case pass. Just checks the normal direction when R>20. ! This fixes the case by "Lhoeltke" on forum. Not a generic one! ! In line ~530 of ViewFactors. #if 1 BLOCK REAL(KIND=dp) :: r0(3), n0(3), rad,dotpr j = 0 k = 0 DO i=1,n Element =>...

Hi This is a strange one. The computation of view factors is one of the older pieces of Elmer and still there seems to be some strange bug. This piece of code has been used alot! I used quite a bit of time today to find the bug but I could not really find the main culprint. What happens is that ~25 ...

Hi Ryan, You jumped quickly to the parallel simulations of Elmer! Maybe you can show some nice color pictures. Do you have a reference for this case? It would be nice addition in "elmer-elmag" repo! Your setting's might not be fully optimal. For example, I think you didn't provide any BC's...

Hi, Linear solver setting are needed when solving for linear system. Usually this relates only to PDE's but it also relates to some postprocessing solvers where, for example, B=Curl(A) is also solved using the finite element method. No nonlinear system setting are needed for this since the primary s...

Hi It was previously possible not to define any linear solver settings. Then the code used defaults which happened to be "direct" solver with "banded" method. We thought that it was a very poor default so the user is now obliged to give one (for the CalcFields) solver. Linear Sys...

Hi I changed the MPI communicator and the solver may now function better. Prior to this it did work in parallel but only if the Coil was active in all partitions. The mortar stuff would be much more difficult to implement fully parallel. The issue there is that all search algorithms should be parall...

Hi I tried to fix this. At one point the saving of point data was made more generic (e.g. it can handle Hcurl elements at the node) but it may have broken this parallel operation which was not guarded by tests. Thanx for reporting. It is good idea to use the "exact" mode that find the elem...

Hi Alberto, The velocity should be the velocity of the particles. If there is no electric field the particles will stay at rest (in this test). When I increase |B| from 1 to 2 the results change at least for me. Maybe you set the electric field to zero and particles never start. This is really the o...

Hi I don't think the "double usage" is anything special. This lies in the very heart of isoparametric elements. Here is one paper from 1968... https://www.sciencedirect.com/science/article/abs/pii/0020768368900310 Now there is a pragmatic reason not improve the geometry description using t...