<div dir="ltr"><div><div><div>Hi,<br></div><br></div>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.<br><br>Can you share a data set that reproduces the issue?<br><br></div>Andy<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 23, 2015 at 8:23 AM, Per Adamsen <span dir="ltr"><<a href="mailto:padamsen@grundfos.com" target="_blank">padamsen@grundfos.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div link="#0563C1" vlink="#954F72" lang="DA">
<div>
<p class="MsoNormal">Hi,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span lang="EN-US">I have been trying to visualize flow using the particle tracers as described at
<a href="http://openfoamwiki.net/index.php/HowTo_use_particle_tracer_in_paraFoam" target="_blank">
http://openfoamwiki.net/index.php/HowTo_use_particle_tracer_in_paraFoam</a>. <u></u>
<u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">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. <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">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.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I have tried the following 64-bit version:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">3.98.1 on Win7<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">4.4.0 on Win7<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">4.2.0 on Linux HPC cluster<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">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.
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Are there any restriction to the data which can be used?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">The following is the trace right up to pressing ‘Play’:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">#### import the simple module from the paraview<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">from paraview.simple import *<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">#### disable automatic camera reset on 'Show'<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">paraview.simple._DisableFirstRenderCameraReset()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># get active source.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoam = GetActiveSource()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># Properties modified on q97_volutefoam<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoam.CaseType = 'Decomposed Case'<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># get active view<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">renderView1 = GetActiveViewOrCreate('RenderView')<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># uncomment following to set a specific view size<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># renderView1.ViewSize = [1025, 773]<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># get color transfer function/color map for 'p'<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">pLUT = GetColorTransferFunction('p')<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># show data in view<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay = Show(q97_volutefoam, renderView1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># trace defaults for the display properties.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay.ColorArrayName = ['POINTS', 'p']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay.LookupTable = pLUT<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay.ScalarOpacityUnitDistance = 0.01089979797382614<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay.SelectInputVectors = ['POINTS', 'U']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay.WriteLog = ''<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># reset view to fit data<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">renderView1.ResetCamera()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># show color bar/color legend<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoamDisplay.SetScalarBarVisibility(renderView1, True)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># get animation scene<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">animationScene1 = GetAnimationScene()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># update animation scene based on data timesteps<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">animationScene1.UpdateAnimationUsingDataTimeSteps()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># get opacity transfer function/opacity map for 'p'<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">pPWF = GetOpacityTransferFunction('p')<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># create a new 'Plane'<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1 = Plane()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># Properties modified on q97_volutefoam<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">q97_volutefoam.CellArrays = ['U', 'k', 'nut', 'omega', 'p', 'gradP', 'myPtot', 'myYPlus', 'uPlus', 'y', 'yPlus']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># Properties modified on plane1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1.XResolution = 100<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1.YResolution = 100<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># show data in view<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1Display = Show(plane1, renderView1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># trace defaults for the display properties.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1Display.ColorArrayName = [None, '']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1Display.SelectInputVectors = ['POINTS', 'Normals']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plane1Display.WriteLog = ''<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># create a new 'ParticleTracer'<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">particleTracer1 = ParticleTracer(Input=q97_volutefoam,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">    </span><span lang="EN-US">SeedSource=plane1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1.SelectInputVectors = ['POINTS', 'U']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># Properties modified on particleTracer1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1.StaticSeeds = 1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1.StaticMesh = 1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># show data in view<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1Display = Show(particleTracer1, renderView1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># trace defaults for the display properties.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">particleTracer1Display.ColorArrayName = ['POINTS', 'p']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">particleTracer1Display.LookupTable = pLUT<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1Display.SelectInputVectors = ['POINTS', 'U']<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1Display.WriteLog = ''<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># hide data in view<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Hide(q97_volutefoam, renderView1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># hide data in view<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Hide(plane1, renderView1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"># show color bar/color legend<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">particleTracer1Display.SetScalarBarVisibility(renderView1, True)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Regards,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Per Adamsen</span><u></u><u></u></p>
</div>
</div>

<br>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" rel="noreferrer" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" rel="noreferrer" target="_blank">http://markmail.org/search/?q=ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><br>
<br></blockquote></div><br></div>