CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Paraview & paraFoam

Paraview slice type - cylinder?

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree7Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   July 11, 2013, 07:18
Default
  #21
New Member
 
Christian Wei
Join Date: Jul 2013
Posts: 4
Rep Power: 3
Nephi is on a distinguished road
hi
thanks for your reply. Im going to try this.
meanwhile i got a hint to use the calculator filter to add an radius column to each point an then use the contour filter for only showing data at a certain radius.
after i got that i flattened the "slice" with an caltulater doing something like this
Code:
0*iHat+(asin(coordsY/(coordsX^2+coordsY^2)^(1/2))*(coordsX^2+coordsY^2)^(1/2))*jHat+coordsZ*kHat
this works verry well for me while i will never analyse more than a half cylinder because i just look at one blade.
with best regards
Christian
wyldckat likes this.
Nephi is offline   Reply With Quote

Old   July 22, 2013, 23:56
Default
  #22
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 298
Rep Power: 8
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
p, li { white-space: pre-wrap; } Hi
it works for me here.
try this

Code:
from paraview import vtk 
import numpy as np  
 
pdi = self.GetPolyDataInput() 
 
c2p=vtkCellDataToPointData() 
c2p.SetInputData(pdi) 
c2p.Update() 
pdi.GetPointData().AddArray(c2p.GetOutput().GetPointData().GetArray(0)) 
 
pdo =  self.GetPolyDataOutput() 
newPoints = vtk.vtkPoints() 
numPoints = pdi.GetNumberOfPoints() 
 
for i in range(0, numPoints): 
    coord = pdi.GetPoint(i) 
    x, y, z = coord[:3] 
    if x > 0: 
      r = np.sqrt(x*x+y*y) 
      t = np.arctan(x/y) + 0.5 * np.pi 
      if y < 0: 
         t = t + 1* np.pi 
    else : 
      r = np.sqrt(x*x+y*y) 
      t = np.arctan(x/y) +  1.5* np.pi 
      if y > 0: 
        t = t - 1* np.pi 
    newPoints.InsertPoint(i, r, t, 12*z) 
pdo.SetPoints(newPoints) 
#pdo.GetCellData().AddArray(pdi.GetCellData().GetArray(0)) 
pdo.GetPointData().AddArray(pdi.GetPointData().GetArray(0))
here I export my result on a cyliner surface to a vtk file by "sample" utility. what I got is data in cell so I convert it to the point. so if you just load OF result directly, you just need to comment some of same
Attached Images
File Type: jpg CylindricalTransform.jpg (77.1 KB, 23 views)

Last edited by waynezw0618; July 23, 2013 at 04:56.
waynezw0618 is offline   Reply With Quote

Old   July 23, 2013, 04:57
Default this is the results
  #23
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 298
Rep Power: 8
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
here is the results from the code
Quote:
Originally Posted by waynezw0618 View Post
p, li { white-space: pre-wrap; } Hi
it works for me here.
try this

Code:
from paraview import vtk 
import numpy as np  
 
pdi = self.GetPolyDataInput() 
 
c2p=vtkCellDataToPointData() 
c2p.SetInputData(pdi) 
c2p.Update() 
pdi.GetPointData().AddArray(c2p.GetOutput().GetPointData().GetArray(0)) 
 
pdo =  self.GetPolyDataOutput() 
newPoints = vtk.vtkPoints() 
numPoints = pdi.GetNumberOfPoints() 
 
for i in range(0, numPoints): 
    coord = pdi.GetPoint(i) 
    x, y, z = coord[:3] 
    if x > 0: 
      r = np.sqrt(x*x+y*y) 
      t = np.arctan(x/y) + 0.5 * np.pi 
      if y < 0: 
         t = t + 1* np.pi 
    else : 
      r = np.sqrt(x*x+y*y) 
      t = np.arctan(x/y) +  1.5* np.pi 
      if y > 0: 
        t = t - 1* np.pi 
    newPoints.InsertPoint(i, r, t, 12*z) 
pdo.SetPoints(newPoints) 
#pdo.GetCellData().AddArray(pdi.GetCellData().GetArray(0)) 
pdo.GetPointData().AddArray(pdi.GetPointData().GetArray(0))
here I export my result on a cyliner surface to a vtk file by "sample" utility. what I got is data in cell so I convert it to the point. so if you just load OF result directly, you just need to comment some of same
Attached Images
File Type: jpg a.jpg (36.4 KB, 29 views)
waynezw0618 is offline   Reply With Quote

Old   October 23, 2013, 06:30
Default Slice with geometry
  #24
New Member
 
Abdalla
Join Date: Oct 2013
Posts: 4
Rep Power: 2
AS_allways is on a distinguished road
Hi all,
first of all, i am a beginner in using paraview.
here at work i have to do some postprocessing for a motor flow simulation. Actually i have to test how to slice with a geometry. I have seen and also successfully tried to cylinder-slice the flow using the programmable filter as shown at the beginning of this thread.
Since the programmable filter can not handle complex geometries i am wondering if i can use a CAD-Geometry (the bend cylinder is an .stl geometry file) to slice with.

best regards,
Abdalla
AS_allways is offline   Reply With Quote

Old   October 26, 2013, 10:42
Default
  #25
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 6,981
Blog Entries: 32
Rep Power: 68
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Greetings Abdalla and welcome to the forum!

The programmable filter is only used because the cylinder shape is not present in the "Slice" filter.

In your case, it's different and somewhat easier to do. Simply load up both data sets (the CFD results and the slice surface) and use the "Resample with DataSet" filter.

I've finished just now writing the following post, which explains in more detail how to do this: Mapping in ParaView post #4

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   February 11, 2014, 09:18
Default Get time in a Programmable Filter paraview 3.14
  #26
New Member
 
German Blesio
Join Date: Feb 2014
Posts: 2
Rep Power: 0
gerble is on a distinguished road
Hi to everyone. I'm new in Paraview (start learning this month for an intership) and I have some complications. I need to program a filter to change over time. I used
Code:
t = inputs[0].GetInformation().Get(vtk.vtkDataObject.DATA_TIME_STEP()
in paraview 3.98. But in paraview 3.14 doesn't work and I couldn't found a way of fixing this "mistake". The reason for using paraview 3.14 is that I need to clip a cylinder.

Thank you very much, and I'm sorry if I have any spelling/gramatical mistakes
gerble is offline   Reply With Quote

Old   February 11, 2014, 14:43
Default
  #27
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 6,981
Blog Entries: 32
Rep Power: 68
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Greetings gerble and welcome to the forum!

I don't have time to test this right now, but some time ago I wrote/updated this post: Annoying issue of automatic "Rescale to Data Range " with paraFoam/paraview 3.12 post #6

After some further developments on the topic, the resulting conclusion was that the only guaranteed way to get the time in all ParaView versions was to use this:
Code:
GetAnimationScene().TimeKeeper.Time
Problem is that this might depend on where it's being used. You might need to add this line to your code in the filter, before everything else:
Code:
from paraview.simple import *
Please let us know if either one of these worked.

Best regards,
Bruno
gerble likes this.
wyldckat is offline   Reply With Quote

Old   February 11, 2014, 22:50
Default
  #28
New Member
 
German Blesio
Join Date: Feb 2014
Posts: 2
Rep Power: 0
gerble is on a distinguished road
Unfortunatly it doesn't work out. I have tried in a Windows 7 system and in a Debian base computer. What happen is that when executing the paraview.simple import this messages appears: "<string>:1: SyntaxWarning: import * only allowed at module level" and the result of GetAnimationScene().TimeKeeper.Time is 0.0 . Before writing in this forum I tried another sentences that requieres this definition and this error always appears. If I don't write that sentece, if I only try to get the time an error appears informing that 'GetAnimationScene' is not defined.

On the other hand, I was thinking on define a variable in the python code so I can later modified through the "Animation View". It would consume more memory to make de animations (I think) but it will probably work equally to what I was trying. Is there a way of doing this? Sorry if I changed the issue of the Thread too much.

Thank you very much for your answer!

Last edited by gerble; February 12, 2014 at 11:47. Reason: spelling mistake
gerble is offline   Reply With Quote

Old   February 16, 2014, 12:10
Default
  #29
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 6,981
Blog Entries: 32
Rep Power: 68
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Hi German,

After doing some trial-and-error-and-research, this programmable filter seems to work in ParaView 3.12.0:
Code:
def GetUpdateTimestep(algorithm):
    """Returns the requested time value, or None if not present"""
    executive = algorithm.GetExecutive()
    outInfo = executive.GetOutputInformation(0)
    if not outInfo.Has(executive.UPDATE_TIME_STEPS()):
       return None
    return outInfo.Get(executive.UPDATE_TIME_STEPS(),0)
 

 # This is the requested time-step. This may not be exactly equal to the 
 # timesteps published in RequestInformation(). Your code must handle that
 # correctly
 req_time = GetUpdateTimestep(self)
 

 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.1*(0.5+req_time))
 cyl.SetTransform(transf)
 cutter.SetCutFunction(cyl)
 cutter.SetInput(inp_copy)
 cutter.Update()
 self.GetOutputDataObject(0).ShallowCopy(cutter.GetOutputDataObject(0))
It's based on these sources of information:
The "req_time" variable is the current time snapshot.


As for global variables, they simply don't work, because the Python workspaces are completely independent, namely the ones on each filter and on the main Python shell.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   May 26, 2014, 08:07
Default
  #30
New Member
 
Christoph Rosebrock
Join Date: Jun 2010
Posts: 1
Rep Power: 0
CROS is on a distinguished road
Hello everyone,
I was looking for a possibility to create radial cuts and the solution presented by Bruno seemed to be a good solution. Using Paraview 4.1.0 I got the same error like akashjangid

Quote:
Originally Posted by akashjangid View Post
Traceback (most recent call last):
File "<string>", line 21, in <module>
File "<string>", line 14, in RequestData
AttributeError: SetInput
Replacing
Code:
cutter.SetInput(inp_copy)
with

Code:
cutter.SetInputData(inp_copy)
keeps the code running in Paraview 4.1.0.

Thanks to all of you.

Christoph
wyldckat likes this.
CROS is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
multiphaseInterFoam: timestep error by simulating a co-extrusion nozzle Quatschinsky OpenFOAM Running, Solving & CFD 7 March 27, 2014 06:08
interFoam/kOmegaSST tank filling with printStackError/Mules simpomann OpenFOAM Running, Solving & CFD 3 February 17, 2014 17:06
Need help with boundary conditions: open to atmosphere Wolle OpenFOAM 2 April 11, 2011 07:32
cgnsToFoam problems with "QUAD_4" cells lentschi OpenFOAM Meshing & Mesh Conversion 1 March 9, 2011 04:49
buoyantSimpleRadiationFoam msarkar OpenFOAM 0 February 15, 2010 06:22


All times are GMT -4. The time now is 12:19.