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

[Andy Bauer]

Hi Per,

I'm glad you're able to move forward with this.

Did you look at https://cmake.org/pipermail/paraview/2014-May/031000.html
for information on the Id Channel Array? That may help.

Best,
Andy

On Thu, Mar 17, 2016 at 3:27 AM, Per Adamsen <padamsen at grundfos.com> wrote:

> 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>
> 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]
> *Sent:* 23. november 2015 16:18
> *To:* Per Adamsen <padamsen at grundfos.com>
> *Cc:* 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>
> 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
>
> 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/b09a6d2f/attachment.html>