movingMeshGran example

Submitted by phsieh2005 on Thu, 02/25/2010 - 14:10

Dear LIGGGHTS group,

I managed to compile LIGGGHTS and pizza.py on OpenSUSE 11.2 64 OS with mpich2.

The first example I ran is the movingMeshGran example. This case created quite a few new stl files (dump###.stl). I will assume that these stl files corresponds to each time step. How to animate the stl files together with the particles? I can animate the particle motion in paraview, but, I do not under what the particles are doing. I guess I will need the stl geometries to understand.

Your assistance will be highly appreciated.

Pei-Ying

Forums:

ckloss's picture

ckloss | Fri, 02/26/2010 - 07:35

Hi Pei,

yes the STL files correspond to the time-steps, you can directly view them in paraview. (The same procedure as if you would open a VTK files, paraview regards the dump###.stl files as belonging together, one per time-step

Regards,
Christoph

alexander.polson | Sat, 04/24/2010 - 17:54

Hi

I also ran this example. I now also have the stl files, but am not having success in visualizing the results.

Any advice on the steps from pizza (granular) to Paraview?

In Paraview I could see only the bucket, but not the spheres if I open directly. Must I do some pizza.py steps before?

Thank you

Alexander

ckloss's picture

ckloss | Sun, 04/25/2010 - 21:21

Hi Alex,

this example outputs:

+ the STL files (the bucket)
+ the dump file

If you use the pizza.py from out website, you can, with the
v.manyGran()
routine,
write the data to VTK. For each time-step, the particle data and the bounding box are written to VTK

So you just open the STL file list (a a whole) and the two VTK file lists - you should then see both the particles and the bucket.

Cheers,
Christoph

alexander.polson | Thu, 04/29/2010 - 21:27

Hi Christoph

Thank you very much!

Your explanation above, as well as the detailed explanation in "readme" in the "gran_pizza-16Mar09.tar_.gz" file helped me to successfully complete the whole process up to ParaView.

Also, I managed to compile under Ubuntu 9.1 using "make fedora".

Kind regards

Alexander

ckloss's picture

ckloss | Fri, 04/30/2010 - 16:13

Hi Paul,

I added the pizza.py doc to the newer version also!

Christoph

Somnath | Wed, 05/05/2010 - 08:40

Hi Christoph,

I am having trouble trying to compile LIGGGHTS (liggghts-1p0) on Mac OS X 10.6+. The earlier LIGGGHTS package (lammps-24Jan10) had compiled smoothly using make myFedora11. I have tried copying Makefile.myFedora11 to the new ligggths-1p0 and compiling, but to no use. This is the error I get with myFedora11 and liggghts1p0:

fix_meshGran.cpp: In member function ‘void LAMMPS_NS::FixMeshGran::calcTriCharacteristics(int, double***, double**, double***, double***, double*)’:
fix_meshGran.cpp:157: error: ‘isnan’ was not declared in this scope

Any ideas?

Thanks,
Som

ckloss's picture

ckloss | Wed, 05/05/2010 - 12:57

Hi Som,

try adding

#include < cmath >
(without the blanks between cmath and the parantheses)

to the include list of these files (can be found at the beginning of the code):

fix_meshGran.cpp
and
fix_wall_gran_hooke_history.cpp

Let me know if that worked

Christoph

Somnath | Thu, 05/06/2010 - 08:49

Hi Christoph,

Including cmath didnot work. However, using std::isnan() instead of isnan() in:

fix_meshGran.cpp
and
fix_wall_gran_hooke_history.cpp

solved the problem. The rest compiled smoothly. Thanks!

Regards
Som

schteppe | Wed, 02/09/2011 - 15:03

Hi Christoph.

I've tried doing as you say but I still can't see any bucket. I've tried all example simulations that contains .stl meshes. Only the simple_chute.stl works to view. The granular pizza conversion works well and the particles becomes visible glyph spheres in Paraview after doing some settings, but no bucket.

To try to view the bucket, I import dump***.stl as dump..stl in Paraview (as I do with the particle data) and then I press Apply in the object inspector. Nothing appears. Under Information I see that the bucket mesh is a polygonal mesh with 6 points. I also see the interesting fact that the bounds for the mesh is only along the X axis (-6 to 6) and zero for the other axes. By using the "Show cube axes" option one can see that this is the case, a white line with ticks appear along the X axis.

I'm running Ubuntu 10.04 and Paraview 3.4.0. Can you help me?

Btw, thanks for all your awesome work.

Stefan

ckloss's picture

ckloss | Wed, 02/09/2011 - 15:54

Hi Stefan,

>>Btw, thanks for all your awesome work.
Thanks :-)

>>I've tried doing as you say but I still can't see any bucket
>>I've tried all example simulations that contains .stl meshes. Only the simple_chute.stl works to view.
I assume this is an issue with your system - try to google "paraview cannot see stl" - that might help?

Cheers,
Christoph

schteppe | Fri, 02/11/2011 - 09:20

A newer version of Paraview did it. I guess the one in the Ubuntu repos is a bit old and got some bugs. Thanks!

Stefan

alexander.polson | Sat, 04/24/2010 - 17:54

Hi

I also ran this example. I now also have the stl files, but am not having success in visualizing the results.

Any advice on the steps from pizza (granular) to Paraview?

Thank you

Alexander

aacusm | Thu, 03/03/2011 - 23:08

Dear Pei-Ying, i'm new user of LIGGGHTS and i trying to compile it on my opensuse11.1 64 OS with mpich2, but i can't. Can you tell me how you did it, please ....

Thanks in advance
Alejandro