Errors when extracting parameters over time

Submitted by NTT1508 on Thu, 07/07/2016 - 04:13

Dear mates,

I am struggling with post-processing by Paraview. There are several matters I would like to carry out on Paraview but it seems to not work properly.

- How to extract the variation of a particle over time, e.g., particle coordinates such as X, Y, Z? I selected a particle (Select Points Through (g)) and used filter/Plot Selection Over Time, but this seems to show the variation of a "fixed" point in the domain over time. When the particle moves insignificantly, this selection shows reasonable result, but when the particle moves much, its parameters does not change accordingly.

- To export the variation of number of particles in a certain part of the whole domain over time, I used Clip/Box to extract a certain zone. I can see the number of particles changing over time in this zone on Paraview, but I dont see any solution to extract the number of particles in the zone effectively.

Please advise how to solve these problems by using Paraview. I believe Christian can give some hints on these issues.

Thank you,

Regards,

richti83's picture

richti83 | Tue, 07/12/2016 - 09:16

Hi,
yes, the "select points function" in paraview uses the internal ID in PV, not the particle ID, when the sequence of rows in dumpfile changes bc. a particle is on another proc (the dumpfiles are ordered by core more or less) the internal PV id renew.
You need to use a filter Threashold by Particle ID, select the ID(s) of interest and plot this over time.

2) Use a Python Programmable filter to access the number of points in the clip

input = self.GetInputDataObject(0, 0)
output = self.GetOutputDataObject(0)
COUNT=input.GetNumberOfPoints()
outputarray = vtk.vtkStringArray()
outputarray.SetName("Text")
outputarray.SetNumberOfTuples(1)
outputarray.SetValue(0,"N=%1.0f" % (COUNT))
output.GetRowData().AddArray(outputarray)

Set type to "vtkTable"

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

NTT1508 | Thu, 07/14/2016 - 08:07

Dear Christian,

That is always great to have your advice. The problem with the ID of particles and extract over time have been solved with respect to your explanation.

Your codes to extract the number of particles in a clipped regions showed an error, but it helps me understand how to approach it. I will improve it (my Python coding will restrain me for a time :(.

Thank you much,

Regards,