Elmer Discussion Forum

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

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

Hi, Peter Thank you for the instruction. I agree that the fitting is not the best thing. I found it difficult to programme a new solver, so I add the suboutine in the "CalcFields.F90". Now the elemtal DG field can be solved at first, as shown below: CoreLoss_ElementField.jpg But, it is wei...

Hi, Now the core loss of each node has been solved. But, I do not know how to get the relationship between the element and its nodes, for the triangular type elements. Is there any variable that denotes such mapping? I can find some data in "mesh.elements", but I do not know how to invoke ...

Hi Peter, Thank you! The DG type of elemental field is necessary. And I still encounter the problem that I cannot define the similar variable Phys of Peddy for elemetal field. I checked the source code "CalcFields .F90" and "FourierLoss.F90", I found I should define another dummy...