Elmer Discussion Forum

Fantastic! I did not try loading the 827 with ElmerGUI yet, but I will put in the patch so I don't notice a problem when I do try. Edit: I did try, and it crashed. Looking at the code, it's not so simple to fix. For example there's a line int facenodes = degree * faceedges; Which gives 8 for 820, bu...

I'm using a C program to create the mesh. Maybe the order I call things has something to do with it. Sounds like you first create the Hexa8 (808) mesh, then tell gmsh to go 2nd order. I'm telling it to start with 2nd order at the beginning. There is a bit Mesh.SecondOrderIncomplete which is default ...

I saw that at some point and it didn't make much sense. Now that I read it multiple times, I'm starting to understand a bit. The files are different now (9 years later), but I can find the places that require changing. I think this will be a lot of fun to play with, thanks! Edit: I am definitely con...

Thanks! I get a bunch of errors ... GetElementDimension: unknown elementtype 0 GetElementDimension: unknown elementtype 0 GetElementDimension: unknown elementtype 0 ... But the final result seems to be form 820 which loads fine. The 827 is in elements.def, that's why I was curious how to get that on...

I can get gmsh to create a 2nd order mesh of hexahedra, but ElmerGrid does not seem to notice. Here's what gmsh outputs: .... Info : [ 90%] Meshing surface 62 order 2 Info : [ 90%] Meshing surface 70 order 2 Info : [100%] Meshing surface 75 order 2 Info : [100%] Meshing volume 1 order 2 Info : [100%...

Finally got what I wanted for the simple test. Now I have a lot more to learn, but I think I'll pause while I'm happy.

I dumped all the x, y, z points using a write statement, and copied that to a file. I then used those to feed into my formula, and realized I wasn't doing what I thought. I also found that paraview is "sticky" and I had to reset the amplitude scale every time I read in a new run. My playin...

Thanks, I will play with that. I dug into the StatElecSolve.F90 and it sure looks like this should have worked. Hopefully I'll figure out what I'm doing wrong soon!

Do you see a ring charge in the potential output? A blob of charge should change the potential distribution - a blob of dielectric not so much. I don't see any change from the charge present vs the charge not preset. So I'm doing something wrong! Edit: are there any examples that have charged partic...

Sure: Function initrho(model, n, dummyin) result(rho) USE DefUtils Implicit None ! function variables Type(Model_t) :: model Integer :: n Real(Kind=dp) :: dummyin, rho ! internal variables Integer :: dim Real(Kind=dp) :: x, y, z, r1, r2, r, rp Logical :: firsttime = .true. ! keep these variables Sav...