[Paraview] ParticleTracer filter makes Paraview crash on 'play'

Thu Mar 17 03:27:32 EDT 2016

Hi Andy,

I have tried with PV 5.0 and I got a lot further than previous (I didn’t test my findings below on earlier PV versions)

The issue with missing particles after 1st timestep was simply because the time-stepping was too coarse making all the particles leave the domain at next time-step.

Adding the temporal interpolator to get much smaller timesteps I managed to get some quite descent animations although I still have some issues with the path lines jumping around which I would assume is caused by the wrong “Id Channel Array” setting.

Regards,
Per

From: Andy Bauer [mailto:andy.bauer at kitware.com]
Sent: 23. februar 2016 16:12
To: Per Adamsen <padamsen at grundfos.com>
Cc: paraview at paraview.org
Subject: Re: [Paraview] ParticleTracer filter makes Paraview crash on 'play'

Hi Per,
Can you share either a Python script or state file which reproduces the problem for PV 5.0? I still have your OpenFOAM file to test with so I don't need that.
Thanks,
Andy

On Tue, Feb 23, 2016 at 2:10 AM, Per Adamsen <padamsen at grundfos.com<mailto:padamsen at grundfos.com>> wrote:
Hi Andy,

I just want to refresh this conversation – sorry for not passing all messages through this mailing list… my mistake.

After I shared the case you wrote that you had tried a few things but with no luck in PV 5.0. As mentioned below I have tried the following 64-bit versions
- 3.98.1 on Win7
- 4.4.0 on Win7
- 4.2.0 on Linux HPC cluster
Regards,
Per

From: Andy Bauer [mailto:andy.bauer at kitware.com<mailto:andy.bauer at kitware.com>]
Sent: 23. november 2015 16:18
To: Per Adamsen <padamsen at grundfos.com<mailto:padamsen at grundfos.com>>
Cc: paraview at paraview.org<mailto:paraview at paraview.org>
Subject: Re: [Paraview] ParticleTracer filter makes Paraview crash on 'play'

Hi,

My guess is that it has to do with OpenFOAM's way of saving out extra fields for the initial time step. Several temporal filters assume that the arrays are ordered the same for different time steps but the OpenFOAM reader doesn't follow this convention.

Can you share a data set that reproduces the issue?
Andy

On Mon, Nov 23, 2015 at 8:23 AM, Per Adamsen <padamsen at grundfos.com<mailto:padamsen at grundfos.com>> wrote:
Hi,

I have been trying to visualize flow using the particle tracers as described at http://openfoamwiki.net/index.php/HowTo_use_particle_tracer_in_paraFoam.

It works when using it on the cavity tutorial but when I try to transfer the setup to one of my own 3D OpenFOAM cases of a pump volute with e.g. 20 timesteps and 900,000 cells Paraview crashes with a segmentation fault when I hit the ‘Play’ button to step through the individual timesteps.

The pipeline setup is exactly as described in the HowTo (linked above) and I have tried both with and without the Temporal Interpolator and on Linux and Windows platform – all with the same crash result.

I have tried the following 64-bit version:
3.98.1 on Win7
4.4.0 on Win7
4.2.0 on Linux HPC cluster

At first I thought the reason for the crash was due to the use of interfaces in the model but I have tried a model without interfaces at all and the result is the same.

Are there any restriction to the data which can be used?

The following is the trace right up to pressing ‘Play’:
#### import the simple module from the paraview
from paraview.simple import *
#### disable automatic camera reset on 'Show'
paraview.simple._DisableFirstRenderCameraReset()

# get active source.
q97_volutefoam = GetActiveSource()

# Properties modified on q97_volutefoam
q97_volutefoam.CaseType = 'Decomposed Case'

# get active view
renderView1 = GetActiveViewOrCreate('RenderView')
# uncomment following to set a specific view size
# renderView1.ViewSize = [1025, 773]

# get color transfer function/color map for 'p'
pLUT = GetColorTransferFunction('p')

# show data in view
q97_volutefoamDisplay = Show(q97_volutefoam, renderView1)
# trace defaults for the display properties.
q97_volutefoamDisplay.ColorArrayName = ['POINTS', 'p']
q97_volutefoamDisplay.LookupTable = pLUT
q97_volutefoamDisplay.ScalarOpacityUnitDistance = 0.01089979797382614
q97_volutefoamDisplay.SelectInputVectors = ['POINTS', 'U']
q97_volutefoamDisplay.WriteLog = ''

# reset view to fit data
renderView1.ResetCamera()

# show color bar/color legend
q97_volutefoamDisplay.SetScalarBarVisibility(renderView1, True)

# get animation scene
animationScene1 = GetAnimationScene()

# update animation scene based on data timesteps
animationScene1.UpdateAnimationUsingDataTimeSteps()

# get opacity transfer function/opacity map for 'p'
pPWF = GetOpacityTransferFunction('p')

# create a new 'Plane'
plane1 = Plane()

# Properties modified on q97_volutefoam
q97_volutefoam.CellArrays = ['U', 'k', 'nut', 'omega', 'p', 'gradP', 'myPtot', 'myYPlus', 'uPlus', 'y', 'yPlus']

# Properties modified on plane1
plane1.XResolution = 100
plane1.YResolution = 100

# show data in view
plane1Display = Show(plane1, renderView1)
# trace defaults for the display properties.
plane1Display.ColorArrayName = [None, '']
plane1Display.SelectInputVectors = ['POINTS', 'Normals']
plane1Display.WriteLog = ''

# create a new 'ParticleTracer'
particleTracer1 = ParticleTracer(Input=q97_volutefoam,
    SeedSource=plane1)
particleTracer1.SelectInputVectors = ['POINTS', 'U']

# Properties modified on particleTracer1
particleTracer1.StaticSeeds = 1
particleTracer1.StaticMesh = 1

# show data in view
particleTracer1Display = Show(particleTracer1, renderView1)
# trace defaults for the display properties.
particleTracer1Display.ColorArrayName = ['POINTS', 'p']
particleTracer1Display.LookupTable = pLUT
particleTracer1Display.SelectInputVectors = ['POINTS', 'U']
particleTracer1Display.WriteLog = ''

# hide data in view
Hide(q97_volutefoam, renderView1)

# hide data in view
Hide(plane1, renderView1)

# show color bar/color legend
particleTracer1Display.SetScalarBarVisibility(renderView1, True)

Regards,
Per Adamsen

_______________________________________________
Powered by www.kitware.com<http://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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20160317/10440333/attachment.html>