Elmer Discussion Forum

Hi, Dr. Raback, I wonder can I simulate the movement of the iron dust in a static or dynamic magnetic field in the ELMER? Which solver should I use? As we know the iron dust is much smaller, and the size is approximate 4-20 microns. In my case, the magnetic field is produced by a coil with direct cu...

Hi Jiahao, Thank you for your hints. Actually, all the scalers can be automatically saved, with or without the "Save Component Results". I do not know the detail reason. Currently, I am trying to combine the Elmer with MATLAB-Simulink. My Researchgate is as follows: https://www.researchgat...

Hi,
I have solved the problem now.

The reason is that the Elmer version with MPI can not be invoked by several MATLAB parallel threads. The MPI of MATLAB conflicts with the MPI of Elmer.
Hence, I re-install the Elmer version without MPI, it works.

BR
Jenwel

Hi Peter, Thank you again for your response. I have opened the matlabpool as: matlabpool 3 and "parfor" for my loop. The ElmerSolver is invoked in my m file of MATLAB as: [status, results]= system(strcat('ElmerSolver case_',num2str(i),'.sif'),'-echo'); I can find that additional 3 threads ...

Hi dear Elmer designers, I want to optimize the electrical machine with Elmer. Thus, many models with different shape and size should be calculated by Elmer. I have tried the Elmer and MPI for parallel computation, but it solves even slower than only using a single CPU core for my model. Now, I can ...

Hi, Peter

I want to ask that how can I obtain the area of each element in "CalcFields.F90"? Is there a variable that denotes the area of the element, in m^2 for the 2-D element?

Thank you in advance!

BR
Jenwel

Hi, Peter Thank you so much for the answers! I understand now! I revised the code as: DO i = 1, GetNOFActive() Element => GetActiveElement(i) Bnode1 = sqrt(EL_MFD % Values(EPerm(Element % DGIndexes(1)) * 3 - 2)**2 + & EL_MFD % Values(EPerm(Element % DGIndexes(1)) * 3 - 1)**2) Bnode2 = sqrt(EL_MF...

Hi ! OMG, after I read the data in the mesh document of my model, I realized that all the elements are divided into TWO main categories, the active elements and boundary elements. In the file "mesh.elements" and "mesh.boundary", I found the number of active elements is 24518, the...

Hi I found the codes in "CalcFields" DO i = 1, GetNOFActive() Element => GetActiveElement(i) n = GetElementNOFNodes() np = n*pSolver % Def_Dofs(GetElementFamily(Element),Element % BodyId,1) nd = GetElementNOFDOFs(uSolver=pSolver) CALL GetElementNodes( Nodes ) BodyId = GetBody() Material =>...

Hi Peter, That's the model I fitted into "CalcFields" now. Your codes basically agree with mine, except several variable names are different. :D Now, the elemental loss can be solved. But the elements in air gap are also calculated by corresponding B field value, which is not expected. Act...