Muhammad , You can try the Clip filter to generate box-shaped subdomains. Counting the number of particles in each cell is a little more complicated. Off the top of my head, I can't think of a filter that does that for you. But, you could write a Programmable Filter that does this. Briefly, the Programmable Filter would:
1. Set up a vtkCellLocator and give it the input to the Programmable Filter (which is assumed to be a subdomain you've created with the Clip filter) 2. Create a cell vtkDataArray (e.g., vtkIntArray) to store the counts of particles falling inside each cell. The array should have the same number of tuples as cells in the subdomain. 3. Iterate through all the points in a second input given to the Programmable Filter. For each point, find the cell containing the point using vtkCellLocator::FindCell(). Increment the count in the cell array you created in step 2 given the cell ID returned by the vtkCellLocator::FindCell() method. There is a lot of information about the Programmable Filter in the ParaView Guide and at https://www.paraview.org/Wiki/Python_Programmable_Filter. Hope that helps, Cory On Tue, Sep 5, 2017 at 2:04 AM, Muhammad Kashif Saeed <kashifsaee...@gmail.com> wrote: > > Hi everyone, > > I have posted this question some days ago but could not get any answer. Now > I am trying to explain it more. > > I am using open source software LIGGGHTS, which uses DEM (discrete element > method), My simulation contains a cylinder and 4 bladed impeller which > rotates inside the cylinder. I have inserted 3 types of particles in the > cylinder and rotate the impeller which causes mixing of particles that > process is done using LIGGGHTS software. > > Now for post processing, i am using ParaView. In Paraview, I can see the > particles mixed each other. but I want to find out the mixing status for > that purpose I want to use Subdomain based mixing index (SMI). > > For SMI calculation, it requires dividing the domain into subdomains like 3D > mesh and to find the total number of particles and number of particles of > each type in each cell of mesh/subdomain. > > Please let me know how can I create these subdomains in ParaView? > > > > _______________________________________________ > Powered by www.kitware.com > > Visit other Kitware open-source projects at > http://www.kitware.com/opensource/opensource.html > > Please keep messages on-topic and check the ParaView Wiki at: > http://paraview.org/Wiki/ParaView > > Search the list archives at: http://markmail.org/search/?q=ParaView > > Follow this link to subscribe/unsubscribe: > http://public.kitware.com/mailman/listinfo/paraview > -- Cory Quammen Staff R&D Engineer Kitware, Inc. _______________________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Please keep messages on-topic and check the ParaView Wiki at: http://paraview.org/Wiki/ParaView Search the list archives at: http://markmail.org/search/?q=ParaView Follow this link to subscribe/unsubscribe: http://public.kitware.com/mailman/listinfo/paraview