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

[Per Adamsen]

Hi Andy,

This puzzles me…


1.       I was following your steps but when I hit ’play’ particles are only shown for initial timestep.

2.       Particles disappear in the following timesteps and nothing is shown at all.

3.       Returning to initial timestep makes paraview crash.

Regards,
Per

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

Hi,
I was able to get around this issue by doing the following:

  1.  Load your data set with just the U field
  2.  Use the Pass Arrays filter to only pass the U point data field.
  3.  Create a seed source (I used the line source)
  4.  Create the Particle Tracer filter

Please let me know if this doesn't work for you.

Best,

Andy

On Mon, Nov 23, 2015 at 10:18 AM, Andy Bauer <andy.bauer at kitware.com<mailto:andy.bauer at kitware.com>> wrote:
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/20151202/4efb8cf8/attachment.html>