Scripting for run same Python code 20-40 times (corresponding to different slices)

[proxy]

Hello all,

I am trying to average the data in the span-wise direction referring to the previous posts:

Averaging in spanwise direction and plotiing the data

Hi, I tried to simulate a 3D NACA aerofoil in OpenFOAM. I want the velocity profile at different section of this aerofoil. To do that I need to average the velocity field in the span wise direction and then plot the velocity at different locations of the aerofoil. In the given figures: figure1 I have shown the 3D section of the aerofoil and to explicitly show what I mean by the span-wise direction; figure 2 is the contour plot of the 2D slice of the aerofoil in which I want to plot the velocity…

Average of multiple slices in a spanwise direction

Hello All, Corresponding to the thread Averaging in spanwise direction and plotiing the data I got the idea that by creating multiple slices and averaging all those slices I will get the spanwise average value of the velocity. I managed to create multiple slices in the spanwise direction but don’t know how to average those slices to get the “averaged spanwise slice” so that after getting the “averaged spanwise slice” I would extract the line plot at different location of the aerofoil. Any sugge…

I am now slicing the span-wise area in 10-20 slices and then will average these slices to get the averaged span-wise slice. Using the help of @ paraview_ran I got a running script for taking the slice which is as follows:

from paraview.simple import *
#### disable automatic camera reset on 'Show'
paraview.simple._DisableFirstRenderCameraReset()

# create a new 'OpenFOAMReader'
demofoam = OpenFOAMReader(FileName='<FOAM_PATH>')
demofoam.MeshRegions = ['CYLINDER']
demofoam.CellArrays = ['U', 'p']

# Properties modified on demofoam
demofoam.CellArrays = ['U']

# get active view
renderView1 = GetActiveViewOrCreate('RenderView')

####Then slice data using the code similar to this:


# find source
e0foam = FindSource('demo.foam')

# create a new 'Plot On Intersection Curves'
plotOnIntersectionCurves1 = PlotOnIntersectionCurves(Input=e0foam)
plotOnIntersectionCurves1.SliceType = 'Plane'

# Properties modified on plotOnIntersectionCurves1.SliceType
plotOnIntersectionCurves1.SliceType.Origin = [0.0, 0.0, 0.9]
plotOnIntersectionCurves1.SliceType.Normal = [0.0, 0.0, 1.0]

# get layout
layout1 = GetLayout()

# save data
SaveData('SAVE_LOACATION/demo.csv', proxy=plotOnIntersectionCurves1)

To work with this I need to run this code (pvpython slice1.py) at different span-wise location for 10 times, instead I want to automate this slicing process by script/bash file, is it easily doable ? Any comment suggestions will be a great help.

[mwestphal (Mathieu Westphal (Kitware))]

@proxy

is it easily doable ? Any comment suggestions will be a great help.

Yes, use a for loop in your script ?

[proxy]

Thanks @mwestphal for your help. Could you please share some link or video to include the looping ?

[mwestphal (Mathieu Westphal (Kitware))]

Python For Loops

[proxy]

Thanks a lot @mwestphal for your prompt reply and sharing this.

[proxy]

@mwestphal I guess this is a general link for the Python, is there any specific video/link dealing with the loop for the slices (closely related to Paraview)?

[emonson (Eric Monson)]

@proxy I don’t know if this will be helpful, but here’s an example of a script I created last year to loop through both time and end point coordinates for a line, gathering the data in a Pandas DataFrame, and finally saving out as a CSV file. I ran it as a macro in Paraview after loading the data source I wanted.

emonson/neutron_scattering/blob/master/PlotOverRadialLines.py

Embedded in the comments is a reference to a video I watched that helped me set up the basics for the script. I know it’s not looping through slices, but hopefully you can get some ideas from how I loop through line end point coordinates.