Elmer Discussion Forum

In Elmer github repository, in commit 4334 dated Aug 2, 2019 Mika replaced 'Relative Reluctivity' by 'Reluctivity', to avoid misinterpretation.

Thank you Mika for your efforts!

With very kind support of Mika, Juhani and Matthias I finally got anisotropic results. First important note : "Relative Reluctivity" should be considered as "Reluctivity". The results are: left side: calculations in Octave 4.2.2, right side: results. results.jpg Material definiti...

Hi! I tried also Reluctivity tensor (directly, not relative) Material 2 ! nonlinear material Name = "Iron" Reluctivity (3,3) = Real \ 1.5915e+01 7.9577e+14 7.9577e+14 \ 7.9577e+14 3.9789e+04 7.9577e+14 \ 7.9577e+14 7.9577e+14 3.9789e+04 !Relative permeabiliy(3,3) = real \ ! 50000.1 0 0 \ !...

Hi Matthias!

line 1880 of elmerfem/fem/src/modules/MagnetoDynamics/WhitneyAVSolver.F90 is:
As I understand, mu variable is magnetic permeability. It seems to be scalar.

With best regards,

Roman

Hi Matthias! I simplified the script simply declaring the relative permeability tensor as: Material 2 ! nonlinear material Name = "Iron" Relative permeabiliy(3,3) = real \ 50000.1 0 0 \ 0 2.1 0 \ 0 0 2.1 Relative Permittivity = 1 End Solver works, but unfortunately, the material still reac...

Hi!

The easiest way is to take the most similar test case (from tests subdirectory) and substitute step-by-step to your mesh.

It is difficult to say something on the base of .sif file only

With best regards,

Roman

Hi Matthias! I have made some experiments during the night. I used working example: mgdyn_anisotropic_rel , I substituted .F90 file by permfile.F90: !------------------------------------------------------------------------------- SUBROUTINE permf(Model, n, X, Y) !------------------------------------...

Hi! unfortunately it is not working. :( I prepared file permefile.F90: SUBROUTINE permf(Model, n, X, Y) !------------------------------------------------------------------------------- USE DefUtils IMPLICIT NONE TYPE(Model_t) :: Model INTEGER :: n REAL(KIND=dp) :: X(*) REAL(KIND=dp), POINTER CONTIG ...

Hi Matthias! thank you for your help! I work under Linux Mint. I tried to build it, and it works! Previously, I have got an error, probably I used wrong command. This time I tried: elmerf90 reluctivity.F90 -o reluctivity.so There is declaration in the code, which I don't understand: REAL(KIND=dp) ::...

There is also example with anisotropic reluctance: mgdyn_anisotropic_rel
It is promising, but I can't build it.