Elmer Discussion Forum

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.researchgate.net/prof...

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 of MATLAB ...

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_MFD % ...

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 number of boundary ...

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 => GetMateri...

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. Actually, the...

Hi, Peter Now I have calculated the elemental B field. I think maybe I can only calculate the core loss of the active elements. In the solver "MagnetoDynamicsCalcFields", I use GetNOFActive() to know the number of all active elements is 24518. But the number of total elements is 76602. Can I get the...