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. |
Greetings Alex and welcome to the forum!
It's possible, but it requires two things:
Instructions:
Bruno |
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. |
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. |
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' Best regards, Bruno |
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. |
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!! |
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 |
Hi wyldckat,
Yes my rotor is located at (0,0,0) and is rotating along z axis. |
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() File "<string>", line 15 cutter.Update() self.GetOutputDataObject(0).ShallowCopy(cutter.Get OutputDataObject(0)) ^ SyntaxError: invalid syntax Also two lines: Code:
cutter.Update() NOT Code:
self.GetOutputDataObject(0).ShallowCopy(cutter.Get OutputDataObject(0)) Code:
self.GetOutputDataObject(0).ShallowCopy(cutter.GetOutputDataObject(0)) Best regards, Alex. |
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 ? |
hey,
I installed paraveiw 3.14 and in that the script worked correctly Thanks !! |
Hello,
I got the cylindrical slice, now i want to flatten that surface any idea how to do that ? Thanks |
Hi Akash,
The ParaView filter "Warp by vector" comes to mind, but I think I've never used it :( Best regards, Bruno |
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 |
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) Code:
transf.RotateZ(90) 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. |
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 |
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 |
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 |
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 12:14. |