CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Paraview & paraFoam (http://www.cfd-online.com/Forums/openfoam-paraview/)
-   -   Paraview slice type - cylinder? (http://www.cfd-online.com/Forums/openfoam-paraview/117815-paraview-slice-type-cylinder.html)

acocic May 15, 2013 12:46

Paraview slice type - cylinder?
 
Hello to everybody,

I computed the flow in axial fan, in one blade passage and I want see the flow structure in meridian section - so I need to make a slice with cylindrical surface over the blades, and then to "flatten" that surface.

Is it possibly to do that in Paraview (paraFoam)? Any help would be more that welcome. Thanks in advance.

Alex.

wyldckat May 18, 2013 14:51

Greetings Alex and welcome to the forum!

It's possible, but it requires two things:
  1. That you are using ParaView 3.12 or 3.14, including the Python capabilities inside ParaView. The reason for this is that in ParaView 3.98, the necessary VTK filter is not present.
  2. That you use the "Programmable Filter" instead of the normal "Slice" filter. The reason for this is because the "Cylinder" is an implicit function that for some reason they forgot to include in ParaView :(
Reference: http://www.paraview.org/pipermail/pa...ry/016019.html

Instructions:
  1. Apply the "Programmable Filter" to the item you want to cut.
  2. Use the following code in the first text edit box:
    Code:

    input = self.GetInputDataObject(0, 0)
    inp_copy = input.NewInstance()
    inp_copy.ShallowCopy(input)
    inp_copy.UnRegister(None)
    cutter = vtk.vtkCutter()
    transf = vtk.vtkTransform()
    transf.RotateX(90)
    cyl = vtk.vtkCylinder()
    cyl.SetCenter(0,0,0)
    cyl.SetRadius(0.4)
    cyl.SetTransform(transf)
    cutter.SetCutFunction(cyl)
    cutter.SetInput(inp_copy)
    cutter.Update()
    self.GetOutputDataObject(0).ShallowCopy(cutter.GetOutputDataObject(0))

  3. Notes on important details:
    • The cylinder is infinite, because it's an implicit function.
    • "cyl.SetRadius(0.4)" - defines the radius.
    • "cyl.SetCenter(0,0,0)" - defines the centre of the cylinder.
    • The transformation is done with these lines of code:
      Code:

      transf = vtk.vtkTransform()
      transf.RotateX(90)

    • If more transformations are needed, see the list of methods here: http://www.vtk.org/doc/release/5.10/...Transform.html
Best regards,
Bruno

acocic May 20, 2013 05:10

Hey Bruno,

It works just perfect!! Thank you very much! By default, I didn't have "Programmable filter" and Python support (I'm using OpenFOAM-1.6-ext and ParaView 3.12.0 which I compiled without Python support, I guess).

I downloaded then Paraview 3.12 again with

wget http://www.paraview.org/files/v3.12/...-x86_64.tar.gz

extract the files, and tried your code with "new" paraview on VTK data which I previously generated with foamToVTK. And then I got perfect cylindrical slice :).

What is left for me is to connect this version of paraview with OpenFOAM, to be called with paraFoam.

Thanks to your other link

http://www.cfd-online.com/Forums/blo...-openfoam.html

it won't be a problem ;)!

Best regards,
Alex.

acocic May 21, 2013 06:31

Hello once again,

I just wanted to add some new insights on this topic. Actually, we don't have to modify paraFoam script, and OpenFOAM bashrc file. This paraView version I mentioned in previous thread supports OpenFOAM files, and you can open files with *.foam extension.

So I just created aliases, something like paraViewWithPython, which points to binary of new paraview. When you start this command in the root of the case, it won't automatically load the case data, but you can open *.foam file which is there, and then use Programmable filter.

Best regards,
Alex.

wyldckat May 21, 2013 17:40

Hi Alex,

You use the following alias:
Code:

alias paraViewWP='touch case.foam && /long/path/to/ParaView-3.12.0/bin/paraview --data=case.foam'
It was based on the one shown here:http://openfoamwiki.net/index.php/Ho...using_ParaView

Best regards,
Bruno

acocic May 22, 2013 02:31

Hi Bruno,

Once again, it works perfect! Now the data are automatically loaded, and Programmable filter & mighty Python support are there!

Thank you once again for all your support and answers!

Best regards,
Alex.

akashjangid June 16, 2013 16:02

Cylindrical Slice in Paraveiw
 
Hello,

I computed the flow through a rotor(Diameter -1m) which is placed inside a wind tunnel.
Now i want to analyze flow at different span length of blades. So for that i need a cylindrical slice.
I am using 3.12.0 ParaVeiw version
I used this code but it does not show any output

input = self.GetInputDataObject(0, 0)
inp_copy = input.NewInstance() inp_copy.ShallowCopy(input)
inp_copy.UnRegister(None)
cutter = vtk.vtkCutter()
transf = vtk.vtkTransform()
transf.RotateX(90)
cyl = vtk.vtkCylinder()
cyl.SetCenter(0,0,0)
cyl.SetRadius(0.4)
cyl.SetTransform(transf)
cutter.SetCutFunction(cyl)
cutter.SetInput(inp_copy)
cutter.Update() self.GetOutputDataObject(0).ShallowCopy(cutter.Get OutputDataObject(0))

Can someone tell me why?
Thanks a lot!!

wyldckat June 16, 2013 16:15

Hi Akash and welcome to the forum!

(I forgot to welcome you to the forum, when I replied to the message you had sent me!)

And thank you for posting your question here!
But there is one big question left unanswered: where exactly in space is your rotor?

Because this code assumes that the rotor is centred in "0,0,0" and its axis is aligned with Z.

If you don't know exactly where the rotor is in space, you can use the source "Cylinder" and move this cylinder to the location where your rotor is. This source is located in the menu "Sources -> Cylinder".

Then, knowing the position of this cylinder, it'll be easier to configure the script to do the section-cut in the right place!

Best regards,
Bruno

akashjangid June 16, 2013 23:36

Hi wyldckat,

Yes my rotor is located at (0,0,0) and is rotating along z axis.

acocic June 17, 2013 02:58

Hello Akash,

I think that you made some mistakes in the code, or it's just because of copying. For me Bruno's code works perfectly. So if you compare with it, there are slight differences. These are the paraview messages for your code, and possible solutions

Remark 1:
File "<string>", line 3
inp_copy = input.NewInstance() inp_copy.ShallowCopy(input)
^
SyntaxError: invalid syntax

This should be two lines,
Code:

inp_copy = input.NewInstance()
inp_copy.ShallowCopy(input)

Remark 2:
File "<string>", line 15
cutter.Update() self.GetOutputDataObject(0).ShallowCopy(cutter.Get OutputDataObject(0))
^
SyntaxError: invalid syntax

Also two lines:
Code:

cutter.Update()
self.GetOutputDataObject(0).ShallowCopy(cutter.Get OutputDataObject(0)

Remark 3: no space in last line
NOT
Code:

self.GetOutputDataObject(0).ShallowCopy(cutter.Get  OutputDataObject(0))
BUT
Code:

self.GetOutputDataObject(0).ShallowCopy(cutter.GetOutputDataObject(0))
If you used foamToVTK, and you analyze VTK data, you should choose vtkPolyData instead "Same as Input" in Output Data Set Type (just bellow Apply Reset Delete)

Best regards,
Alex.

akashjangid June 17, 2013 04:23

Hello acocic,

I did the changes in the script as told by you.
But it is still showing some errors

Traceback (most recent call last):
File "<string>", line 21, in <module>
File "<string>", line 14, in RequestData
AttributeError: SetInput

Actually I'm a beginner in scripting python for ParaView.
Please help me with this ? Do i have to add any vtk modules for this ?

akashjangid June 17, 2013 13:24

hey,

I installed paraveiw 3.14 and in that the script worked correctly
Thanks !!

akashjangid June 18, 2013 00:44

Hello,

I got the cylindrical slice, now i want to flatten that surface
any idea how to do that ?
Thanks

wyldckat June 20, 2013 16:57

Hi Akash,

The ParaView filter "Warp by vector" comes to mind, but I think I've never used it :(

Best regards,
Bruno

Nephi July 6, 2013 10:52

empty output
 
hello
i have to do almoust the same task as Aleksandar and was verry happy to find your post and so i tried out the programmable filter on a dummycylinder that ist orientated along the y-axis as default and copied in the code into the programmable filter. there are no errors comming up but i do not get any results. the output geometry has 0 points and 0 cells.
is there anythig special i have to be aware of?
i am using paraview 3.14 on linux and on windows
with best regards
Christian

acocic July 7, 2013 05:09

Sharing some experience :-)
 
4 Attachment(s)
Hello Christian,

For me, it was very straightforward. Main flow was in Z-direction, so the axis of the cylinder matched with Z-axis, and Bruno's original code worked perfectly. After it I tried the code with my other case (LES in pipe, where X-axis was the symmetry), and the code didn't work. The solution was to change to rotational axis in

Code:

transf.RotateX(90)
to

Code:

transf.RotateZ(90)
and to put (leave) the center of cylinder at (0, 0, 0), which is actually pipe inlet. The length of the pipe was 0.1m and if I put the center of cylinder in the middle of the domain, (0.05, 0, 0), I also got nothing, On the other hand, if I delete "transform.Rotate" line, I've got cylinder which axis is parallel to Y, and which is shown no matter where the center is located.

So basically, you just have to play around with two lines of the code (cylinder center and axis of rotation) and see what's the most suitable in your case. Maybe, first delete the line with transf.Rotate and see what you'll get (the axis of cylinder should be in Y-dir).

Best regards,
Alex.

Nephi July 7, 2013 12:29

Hi Alex
thank you for your fast reply. But unfortunatly it still does not work.
I tryed out all rotations and all cylinders are centered in 0,0,0
with best regards
Christian

wyldckat July 7, 2013 14:24

Greetings to all!

@Christian: could you please give some indications are to where your geometry/data is located in 3D space? As well as the limits of said geometry/data?

Best regards,
Bruno

Nephi July 8, 2013 03:45

Hi Bruno,
i just added da cylindersurce centered at 0,0,0 with high 1 and a radius of 0.5 , which is symmetric around the y-axis.
I tried cutting with a rotation around x-axis, z-axis and wothout rotarion. I also changed the radius to 0.5, were i expected to get the whole cylindersurface as result.
with best regards
Christian

wyldckat July 10, 2013 16:29

Hi Christian,

Oddly enough, it does not work with a cylinder source. I have no idea why :(

But it does work with a sphere source! You can use the default settings and it will work!

Best regards,
Bruno


All times are GMT -4. The time now is 07:21.