force stream lines to go up to end of simulation domain

ParaView Support
1

Hello,
I have created stream lines from a custom source,
the simulation is a tube with an internal mixer. (with axial direction be Z)
for the source of the stream lines,

  1. load the inlet patch
  2. use extract block
  3. use slice in the extract block in X direction by 10 elements
  4. use slice in the previous created slice in Y direction by 10 elements this creates a homogeneous grid

this works like charm, with the only issue that when I try to ‘show’ the distribution of the stream lines at the end of the tube (at 0.99 zMax), using slice, I am getting a smaller number of points than the ones from the source, some of them are stoped over the tube.
I have played with the limits of the stream line filter and went up to:

MaximumStepLength = 0.05
MaximumSteps = 200000
MaximumStreamlineLength = 1000000.0

when the tube only makes 0.05x0.05x1 in total dimensions.
the number of ‘seeds’ lost is quite small (~14% see photo down the post) but I would like to know if there is any other option/parameter to enable the seeds that are stopped in the walls to continue and exit the reactor.

image
image853×446 21.2 KB

thanks in advance,

2

They probably go out of domain, are you using the streamtracer ?

3

Hello Mathieu,
Sorry for the late response,
I am using stream lines from custom source.
I feel from the images that some of streams ‘stops’ over some walls I can imagine that is the 0 velocity from the wall, but not sure… As even if I increase the parameters that I post it this stream lines do not grow more.
Regards,

4

Please share data and a statefile

5

Hello Mathieu,
I cleaned the case to the bare bones.
here it is:
https://filesender.renater.fr/?s=download&token=a45e660e-3c34-4470-b2fc-9e3ee235b54e
there is a script “paraviewNew.py” inside that is the one I am running to get the streamlines&points. hope it helps.
(and as you can see if you check the length of pointsOfPatchData and pointsAtOutletData they are 7660 and 6592 respectively for a tube with one inlet and one outlet so, no particles should ‘go out’ of the system without going thought the outlet( where the pointsAtOutletData is taken))
sorry didnt saw that you asked also for a statefile, the python script will give the exact result i am facing, in any case you prefere here is the statefile
stateFile.pvsm (136.8 KB)
for example here are two streamlines that ‘stop’ at one moment indeferently of the values that I put in the stream line from custom source parameters:

image
image774×523 75.2 KB

6

Hello Mathieu,
any solution that you can imagine for this? is this a bug?
I played with every parameter in the filter (stream lines from custom source) that i could find and nothing solved it…
thanks in advance,

7

I tried your statefile with 5.11 and 5.12 and got segaults everytime. More investigation needed.

In any case, you may want to switch to point locator if not already.

8

Hello Mathieu,
what do you mean by segaults?
same question to the point locator, you mean the path in the python script? if that, yeah, I changed to:

foamFile='sim.foam'
casePath=str(__file__).replace('\\','/')
casePath=casePath.replace(casePath.split('/')[-1],'')
path=casePath+foamFile
9

I’m using the pvsm.

By segfault I mean a crash.

10

oh, my bad, maybe because not enought RAM? as it is creating a lot of streamlines, here is the same state with way less.
stateWithLessStreamLines.pvsm (236.6 KB) (this but even the previous pvsm file runned correctly for me over tuxedOS (ubuntu based) in paraview 5.12.0)
and even with less number of stream lines (less seed points) it is always the same issue of having lines that stop…

11

I dont see any streamlines stopping using 5.12.0:

image
image1171×895 29.8 KB

12

hello Mathieu,
have a look at the inlet and at the outlet you will see different number of streams,
also in the figure pretty sure it can be seen that is being cutted

image
image1485×1136 96.9 KB

image
image1183×804 56.6 KB

13

They go out of domain:

image
image2129×1373 118 KB

14

Hello mathieu,
thanks for the insight, but how can I resolve this? as the velocity on the wall (cylindrical surface) is equal to 0 in each direction, and this should not happen as it is not physically accurate. is there any way to stop the stream lines of ‘going out of the domain’?
the same happends over the mixer with one of the stream lines:
image

regards,

15

The particles has residual velocity so the 0 velocity at the surface means that it will keep going out slowly but surely. You need to have a small invward velocity at the surface to avoid that.

Alternatively, use a streamtracer that support boundary interaction like the LagrandianParticleTracker.

More complex to setup though: https://www.kitware.com/vtk-and-paraview-add-highly-configurable-particle-tracker/

1 Like