Hi Shuhao,
Thank you for your response. Yes, I found the TemporalStatistics filter as
well. Could you please let me know how I can save it back as a .vtu file?
Best Regards,
Parastou
From: ParaView on behalf of Shuhao Wu
Sent: Monday, October 23, 2017 4:46:24 PM
To: paraview@paraview.org
Subject: Re: [Paraview] Calculating time-averaged shear stress
I believe you're looking for the TemporalStatistics filter. Specifically
what you want maybe something along the lines of:
temporal_statistics = TemporalStatistics(Input=vtk_data)
temporal_statistics.ComputeAverage = 1
temporal_statistics.ComputeMinimum = 0
temporal_statistics.ComputeMaximum = 0
temporal_statistics.ComputeStandardDeviation = 0
This should make the computation slightly faster as you don't need to
compute and store the stddev, min, and max.
Shuhao
On 2017-10-23 02:36 PM, Eslami, Parastou wrote:
> Hi Everyone,
>
>
> I'm very new to Paraview and python scripting. I have a series of time
> dependent .vtu files and would like to
>
>1. Read the vtu files for each time step
>2. calculate the magnitude of shear stress
>3. calculate time-averaged shear stress for all the elements in my
> geometry
>4. output/write a vtu file and load it back to paraview.
>
> I think I have been able to successfully implement steps 1 and 2 (please see
> bellow) but am completely clueless on how to do 3 and 4. Could anyone help
> me with this? I also should mention that I have tried the Temporal Statistics
> Filter and it seems to be working fine, but would like to have a way to do it
> with pythons scripting as well.
>
>
> Thanks,
> Parastou
>
>
> import the simple module from the paraview
> from paraview.simple import *
> import numpy as np
> disable automatic camera reset on 'Show'
> paraview.simple._DisableFirstRenderCameraReset()
>
> dumm_vWSS = []
>
> for i in range(4500, 4550, 50):
> f = "\...\.../all_results_0" + str(i) + ".vtu"
> print i
> # create a new 'XML Unstructured Grid Reader'
> all_results = XMLUnstructuredGridReader(FileName=f)
> all_results.CellArrayStatus = ['GlobalElementID']
> all_results.PointArrayStatus = ['GlobalNodeID', 'vinplane_traction',
> 'vWSS', 'average_speed', 'average_pressure', 'pressure', 'velocity',
> 'timeDeriv']
> # create a new 'Calculator'
> calculator1 = Calculator(Input=all_results)
> calculator1.Function = ''
> # Properties modified on calculator1
> calculator1.ResultArrayName = 'WSS mag'
> calculator1.Function = 'mag(vWSS)'
> # get the data points
> SetActiveSource(calculator1)
> calculator1.UpdatePipeline()
> vWSSdata=servermanager.Fetch(calculator1)
> vWSSPoints =vWSSdata.GetPoints()
> vWSSPointData =vWSSdata.GetPointData()
> vWSSArrayData=vWSSPointData.GetArray('WSS mag')
> dumm_vWSS.append(vWSSArrayData)
> print dumm_vWSS
>
>
>
> Parastou Eslami, PhD
> Post-doctoral Fellow, Harvard University
> Department of Radiology, Massachusetts General Hospital
>
>
>
>
>
>
>
>
> The information in this e-mail is intended only for the person to whom it is
> addressed. If you believe this e-mail was sent to you in error and the e-mail
> contains patient information, please contact the Partners Compliance HelpLine
> at
> http://www.partners.org/complianceline . If the e-mail was sent to you in
> error
> but does not contain patient information, please contact the sender and
> properly
> dispose of the e-mail.
>
>
>
> ___
> Powered by www.kitware.com<http://www.kitware.com>
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at:
> http://paraview.org/Wiki/ParaView
>
> Search the list archives at: http://markmail.org/search/?q=ParaView
>
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/mailman/listinfo/paraview
>
___
Powered by www.kitware.com<http://www.kitware.com>
Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html
Please keep messages on-topic and check the ParaView Wiki at:
http://paraview.org/Wiki/ParaView
Search the list archives at: http://markmail.org/search/?q=ParaView
Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview
___
Powered by www.kitware.com
Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html
Please keep messages on-topic and check the ParaView Wiki at:
http://paraview.org/Wiki/ParaView
Search the list archives at: http://markmail.org/search/?q=ParaView
Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview