[cfdFreshman]

Hello CFD community!

After searching and reading many articles I decided to register in this forum. I hope for some help of specialists.
This is part of the work I am doing during my internship at a development company.

My task - case description
Based on a similar script, which is able to take animations of solved cases, I have to take screens/animation(s) of a steady-state case. This case is solved with OpenFOAM and multiple processors ("mpirun"). The solver creates folders with an initiation timestep or rather the last timestep of the previous folder. Inside the folders the mesh is getting stretched while simulating the piston movement in 0.1 degree of crank angle. The last timestep is copied in the next folder if a specific accuracy can not be reached anymore. All timesteps between first and last in each folder are deleted after copying the last timestep. These steps are to be recorded too.
For this reason the subordinated folders of the case are existent decomposed and the screenshots/animation have to be saved during the solver is operating.

My solution on the whole
Course of action so far is starting a bash shell script which prepares the animation, then it takes a list of all subordinated folders and in final step it excecutes a traced python file of ParaView including ("all") settings with pvbatch for all of the listed folders. In a former version of the script I generated the animation from all screenshots. Now it is integrated in the .py, but it if its necessary I can go back on this. In this new instance I think the part of listing all folders is not required, however this is not the problem at the moment.

The challenge
Like I said, I traced my preparation of the geometry with all properties and the screens/animation in ParaView for all available timesteps. When I now excecute the script it does not make pictures in the correct quantitiy and field according to number and range of timesteps.
My assumption is that ParaView truely does not trace all properties. So I have to set them for myself.

Your knowlege is in demand
I would be very pleased, if could tell me your experience on adding some additional commands for my problem. Supporting I'll provide the scripts in the following posts as far as possible.

Thank you very much
cfdFreshman

[cfdFreshman]

rangeMin=0
rangeMax=110
folderName=caseXYZ

if [ ! -d ../$folderName ]; then

echo "Error"
return

fi

workingDir=`pwd`

curCase=1
t=1

rm -rf pics
mkdir pics
rm -rf clips
mkdir clips

for c in subcase* ;

do

cd $c
currentFolder=`pwd`
/cax/sw/OPENFOAM/LINUX_x86_64/ParaView/ParaView-4.4.0-Qt4-Linux-64bit/bin/pvbatch --use-offscreen-rendering $workingDir/t_komplett3.py $currentFolder U $curCase $workingDir $rangeMin $rangeMax

cd $workingDir
curCase=$((curCase+1))

# for picture in pics/anim_*;
# do
# mv -f $picture pics/cut_$t.jpg
# done
# t=1$t

done

#convert -delay 50 -loop 0 pics/cut*.jpg pics/animation.gif
#rm -f cut*

[cfdFreshman]

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

# create a new 'OpenFOAMReader'
casefoam = OpenFOAMReader(FileName='/usr2/tmp/csnickli/Test_paraview_workflow/case0571/case.foam')
casefoam.CaseType = 'Reconstructed Case'
casefoam.Createcelltopointfiltereddata = 1
casefoam.Adddimensionalunitstoarraynames = 0
casefoam.MeshRegions = []
casefoam.CellArrays = []
casefoam.PointArrays = []
casefoam.LagrangianArrays = []
casefoam.Cachemesh = 1
casefoam.Decomposepolyhedra = 1
casefoam.ListtimestepsaccordingtocontrolDict = 0
casefoam.LagrangianpositionsareinOF13binaryformat = 0
casefoam.Readzones = 0

# Properties modified on casefoam
casefoam.CaseType = 'Decomposed Case'
casefoam.MeshRegions = ['internalMesh']
casefoam.CellArrays = ['CO', 'CO2', 'H', 'H2', 'H2O', 'N', 'N2', 'NO', 'O', 'O2', 'OH', 'RGA', 'T', 'U', 'Ytf', 'Yto2', 'Z2', 'alphat', 'b', 'cellMotionU', 'epsilon', 'k', 'mut', 'p', 'reactFuel']

# get active view
renderView1 = GetActiveViewOrCreate('RenderView')
# uncomment following to set a specific view size
# renderView1.ViewSize = [1460, 912]

# get color transfer function/color map for 'p'
pLUT = GetColorTransferFunction('p')
pLUT.LockDataRange = 0
pLUT.InterpretValuesAsCategories = 0
pLUT.ShowCategoricalColorsinDataRangeOnly = 0
pLUT.RescaleOnVisibilityChange = 0
pLUT.EnableOpacityMapping = 0
pLUT.RGBPoints = [88616.515625, 0.0, 0.0, 1.0, 104562.15625, 1.0, 0.0, 0.0]
pLUT.UseLogScale = 0
pLUT.ColorSpace = 'HSV'
pLUT.UseBelowRangeColor = 0
pLUT.BelowRangeColor = [0.0, 0.0, 0.0]
pLUT.UseAboveRangeColor = 0
pLUT.AboveRangeColor = [1.0, 1.0, 1.0]
pLUT.NanColor = [0.498039215686, 0.498039215686, 0.498039215686]
pLUT.Discretize = 1
pLUT.NumberOfTableValues = 256
pLUT.ScalarRangeInitialized = 1.0
pLUT.HSVWrap = 0
pLUT.VectorComponent = 0
pLUT.VectorMode = 'Magnitude'
pLUT.AllowDuplicateScalars = 1
pLUT.Annotations = []
pLUT.ActiveAnnotatedValues = []
pLUT.IndexedColors = []

# show data in view
casefoamDisplay = Show(casefoam, renderView1)
# trace defaults for the display properties.
casefoamDisplay.CubeAxesVisibility = 0
casefoamDisplay.Representation = 'Surface'
casefoamDisplay.AmbientColor = [1.0, 1.0, 1.0]
casefoamDisplay.ColorArrayName = ['POINTS', 'p']
casefoamDisplay.DiffuseColor = [1.0, 1.0, 1.0]
casefoamDisplay.LookupTable = pLUT
casefoamDisplay.MapScalars = 1
casefoamDisplay.InterpolateScalarsBeforeMapping = 1
casefoamDisplay.Opacity = 1.0
casefoamDisplay.PointSize = 2.0
casefoamDisplay.LineWidth = 1.0
casefoamDisplay.Interpolation = 'Gouraud'
casefoamDisplay.Specular = 0.0
casefoamDisplay.SpecularColor = [1.0, 1.0, 1.0]
casefoamDisplay.SpecularPower = 100.0
casefoamDisplay.Ambient = 0.0
casefoamDisplay.Diffuse = 1.0
casefoamDisplay.EdgeColor = [0.0, 0.0, 0.5]
casefoamDisplay.BackfaceRepresentation = 'Follow Frontface'
casefoamDisplay.BackfaceAmbientColor = [1.0, 1.0, 1.0]
casefoamDisplay.BackfaceDiffuseColor = [1.0, 1.0, 1.0]
casefoamDisplay.BackfaceOpacity = 1.0
casefoamDisplay.Position = [0.0, 0.0, 0.0]
casefoamDisplay.Scale = [1.0, 1.0, 1.0]
casefoamDisplay.Orientation = [0.0, 0.0, 0.0]
casefoamDisplay.Origin = [0.0, 0.0, 0.0]
casefoamDisplay.Pickable = 1
casefoamDisplay.Texture = None
casefoamDisplay.Triangulate = 0
casefoamDisplay.NonlinearSubdivisionLevel = 1
casefoamDisplay.CubeAxesColor = [1.0, 1.0, 1.0]
casefoamDisplay.CubeAxesCornerOffset = 0.0
casefoamDisplay.CubeAxesFlyMode = 'Closest Triad'
casefoamDisplay.CubeAxesInertia = 1
casefoamDisplay.CubeAxesTickLocation = 'Inside'
casefoamDisplay.CubeAxesXAxisMinorTickVisibility = 1
casefoamDisplay.CubeAxesXAxisTickVisibility = 1
casefoamDisplay.CubeAxesXAxisVisibility = 1
casefoamDisplay.CubeAxesXGridLines = 0
casefoamDisplay.CubeAxesXTitle = 'X-Axis'
casefoamDisplay.CubeAxesUseDefaultXTitle = 1
casefoamDisplay.CubeAxesYAxisMinorTickVisibility = 1
casefoamDisplay.CubeAxesYAxisTickVisibility = 1
casefoamDisplay.CubeAxesYAxisVisibility = 1
casefoamDisplay.CubeAxesYGridLines = 0
casefoamDisplay.CubeAxesYTitle = 'Y-Axis'
casefoamDisplay.CubeAxesUseDefaultYTitle = 1
casefoamDisplay.CubeAxesZAxisMinorTickVisibility = 1
casefoamDisplay.CubeAxesZAxisTickVisibility = 1
casefoamDisplay.CubeAxesZAxisVisibility = 1
casefoamDisplay.CubeAxesZGridLines = 0
casefoamDisplay.CubeAxesZTitle = 'Z-Axis'
casefoamDisplay.CubeAxesUseDefaultZTitle = 1
casefoamDisplay.CubeAxesGridLineLocation = 'All Faces'
casefoamDisplay.CustomBounds = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
casefoamDisplay.CustomBoundsActive = [0, 0, 0]
casefoamDisplay.OriginalBoundsRangeActive = [0, 0, 0]
casefoamDisplay.CustomRange = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
casefoamDisplay.CustomRangeActive = [0, 0, 0]
casefoamDisplay.UseAxesOrigin = 0
casefoamDisplay.AxesOrigin = [0.0, 0.0, 0.0]
casefoamDisplay.CubeAxesXLabelFormat = '%-#6.3g'
casefoamDisplay.CubeAxesYLabelFormat = '%-#6.3g'
casefoamDisplay.CubeAxesZLabelFormat = '%-#6.3g'
casefoamDisplay.StickyAxes = 0
casefoamDisplay.CenterStickyAxes = 0
casefoamDisplay.SelectionCellLabelBold = 0
casefoamDisplay.SelectionCellLabelColor = [0.0, 1.0, 0.0]
casefoamDisplay.SelectionCellLabelFontFamily = 'Arial'
casefoamDisplay.SelectionCellLabelFontSize = 18
casefoamDisplay.SelectionCellLabelItalic = 0
casefoamDisplay.SelectionCellLabelJustification = 'Left'
casefoamDisplay.SelectionCellLabelOpacity = 1.0
casefoamDisplay.SelectionCellLabelShadow = 0
casefoamDisplay.SelectionPointLabelBold = 0
casefoamDisplay.SelectionPointLabelColor = [1.0, 1.0, 0.0]
casefoamDisplay.SelectionPointLabelFontFamily = 'Arial'
casefoamDisplay.SelectionPointLabelFontSize = 18
casefoamDisplay.SelectionPointLabelItalic = 0
casefoamDisplay.SelectionPointLabelJustification = 'Left'
casefoamDisplay.SelectionPointLabelOpacity = 1.0
casefoamDisplay.SelectionPointLabelShadow = 0
casefoamDisplay.ScalarOpacityUnitDistance = 0.0018799337935833274
casefoamDisplay.SelectMapper = 'Projected tetra'

# reset view to fit data
renderView1.ResetCamera()

# show color bar/color legend
casefoamDisplay.SetScalarBarVisibility(renderView1 , True)

# get opacity transfer function/opacity map for 'p'
pPWF = GetOpacityTransferFunction('p')
pPWF.Points = [88616.515625, 0.0, 0.5, 0.0, 104562.15625, 1.0, 0.5, 0.0]
pPWF.AllowDuplicateScalars = 1
pPWF.ScalarRangeInitialized = 1

# get animation scene
animationScene1 = GetAnimationScene()

# update animation scene based on data timesteps
animationScene1.UpdateAnimationUsingDataTimeSteps( )

# hide color bar/color legend
casefoamDisplay.SetScalarBarVisibility(renderView1 , False)

# create a new 'Slice'
slice1 = Slice(Input=casefoam)
slice1.SliceType = 'Plane'
slice1.Crinkleslice = 0
slice1.Triangulatetheslice = 1
slice1.SliceOffsetValues = [0.0]

# init the 'Plane' selected for 'SliceType'
slice1.SliceType.Origin = [0.018561832934210543, -0.05640140175819397, 0.010244861245155334]
slice1.SliceType.Normal = [1.0, 0.0, 0.0]
slice1.SliceType.Offset = 0.0

# toggle 3D widget visibility (only when running from the GUI)
Hide3DWidgets(proxy=slice1)

# show data in view
slice1Display = Show(slice1, renderView1)
# trace defaults for the display properties.
slice1Display.CubeAxesVisibility = 0
slice1Display.Representation = 'Surface'
slice1Display.AmbientColor = [1.0, 1.0, 1.0]
slice1Display.ColorArrayName = ['POINTS', 'p']
slice1Display.DiffuseColor = [1.0, 1.0, 1.0]
slice1Display.LookupTable = pLUT
slice1Display.MapScalars = 1
slice1Display.InterpolateScalarsBeforeMapping = 1
slice1Display.Opacity = 1.0
slice1Display.PointSize = 2.0
slice1Display.LineWidth = 1.0
slice1Display.Interpolation = 'Gouraud'
slice1Display.Specular = 0.0
slice1Display.SpecularColor = [1.0, 1.0, 1.0]
slice1Display.SpecularPower = 100.0
slice1Display.Ambient = 0.0
slice1Display.Diffuse = 1.0
slice1Display.EdgeColor = [0.0, 0.0, 0.5]
slice1Display.BackfaceRepresentation = 'Follow Frontface'
slice1Display.BackfaceAmbientColor = [1.0, 1.0, 1.0]
slice1Display.BackfaceDiffuseColor = [1.0, 1.0, 1.0]
slice1Display.BackfaceOpacity = 1.0
slice1Display.Position = [0.0, 0.0, 0.0]
slice1Display.Scale = [1.0, 1.0, 1.0]
slice1Display.Orientation = [0.0, 0.0, 0.0]
slice1Display.Origin = [0.0, 0.0, 0.0]
slice1Display.Pickable = 1
slice1Display.Texture = None
slice1Display.Triangulate = 0
slice1Display.NonlinearSubdivisionLevel = 1
slice1Display.CubeAxesColor = [1.0, 1.0, 1.0]
slice1Display.CubeAxesCornerOffset = 0.0
slice1Display.CubeAxesFlyMode = 'Closest Triad'
slice1Display.CubeAxesInertia = 1
slice1Display.CubeAxesTickLocation = 'Inside'
slice1Display.CubeAxesXAxisMinorTickVisibility = 1
slice1Display.CubeAxesXAxisTickVisibility = 1
slice1Display.CubeAxesXAxisVisibility = 1
slice1Display.CubeAxesXGridLines = 0
slice1Display.CubeAxesXTitle = 'X-Axis'
slice1Display.CubeAxesUseDefaultXTitle = 1
slice1Display.CubeAxesYAxisMinorTickVisibility = 1
slice1Display.CubeAxesYAxisTickVisibility = 1
slice1Display.CubeAxesYAxisVisibility = 1
slice1Display.CubeAxesYGridLines = 0
slice1Display.CubeAxesYTitle = 'Y-Axis'
slice1Display.CubeAxesUseDefaultYTitle = 1
slice1Display.CubeAxesZAxisMinorTickVisibility = 1
slice1Display.CubeAxesZAxisTickVisibility = 1
slice1Display.CubeAxesZAxisVisibility = 1
slice1Display.CubeAxesZGridLines = 0
slice1Display.CubeAxesZTitle = 'Z-Axis'
slice1Display.CubeAxesUseDefaultZTitle = 1
slice1Display.CubeAxesGridLineLocation = 'All Faces'
slice1Display.CustomBounds = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
slice1Display.CustomBoundsActive = [0, 0, 0]
slice1Display.OriginalBoundsRangeActive = [0, 0, 0]
slice1Display.CustomRange = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
slice1Display.CustomRangeActive = [0, 0, 0]
slice1Display.UseAxesOrigin = 0
slice1Display.AxesOrigin = [0.0, 0.0, 0.0]
slice1Display.CubeAxesXLabelFormat = '%-#6.3g'
slice1Display.CubeAxesYLabelFormat = '%-#6.3g'
slice1Display.CubeAxesZLabelFormat = '%-#6.3g'
slice1Display.StickyAxes = 0
slice1Display.CenterStickyAxes = 0
slice1Display.SelectionCellLabelBold = 0
slice1Display.SelectionCellLabelColor = [0.0, 1.0, 0.0]
slice1Display.SelectionCellLabelFontFamily = 'Arial'
slice1Display.SelectionCellLabelFontSize = 18
slice1Display.SelectionCellLabelItalic = 0
slice1Display.SelectionCellLabelJustification = 'Left'
slice1Display.SelectionCellLabelOpacity = 1.0
slice1Display.SelectionCellLabelShadow = 0
slice1Display.SelectionPointLabelBold = 0
slice1Display.SelectionPointLabelColor = [1.0, 1.0, 0.0]
slice1Display.SelectionPointLabelFontFamily = 'Arial'
slice1Display.SelectionPointLabelFontSize = 18
slice1Display.SelectionPointLabelItalic = 0
slice1Display.SelectionPointLabelJustification = 'Left'
slice1Display.SelectionPointLabelOpacity = 1.0
slice1Display.SelectionPointLabelShadow = 0

# hide data in view
Hide(casefoam, renderView1)

# show color bar/color legend
slice1Display.SetScalarBarVisibility(renderView1, True)

# Properties modified on casefoam
casefoam.CellArrays = ['CO', 'CO2', 'H', 'H2', 'H2O', 'N', 'N2', 'NO', 'O', 'O2', 'OH', 'RGA', 'T', 'U', 'Ytf', 'Yto2', 'Z2', 'alphat', 'b', 'cellMotionU', 'epsilon', 'k', 'mut', 'p', 'reactFuel', 'Chi', 'EqT', 'T_inl', 'T_int', 'T_out', 'Ydefault', 'Z', 'dRoHR', 'fuel', 'liquidFrom0', 'p_inl', 'p_int', 'p_out', 'rho', 'voidFrac']

# hide color bar/color legend
slice1Display.SetScalarBarVisibility(renderView1, False)

# set scalar coloring
ColorBy(slice1Display, ('POINTS', 'U'))

# rescale color and/or opacity maps used to include current data range
slice1Display.RescaleTransferFunctionToDataRange(T rue)

# show color bar/color legend
slice1Display.SetScalarBarVisibility(renderView1, True)

# get color transfer function/color map for 'U'
uLUT = GetColorTransferFunction('U')
uLUT.LockDataRange = 0
uLUT.InterpretValuesAsCategories = 0
uLUT.ShowCategoricalColorsinDataRangeOnly = 0
uLUT.RescaleOnVisibilityChange = 0
uLUT.EnableOpacityMapping = 0
uLUT.RGBPoints = [0.0, 0.0, 0.0, 1.0, 103.3250007861863, 1.0, 0.0, 0.0]
uLUT.UseLogScale = 0
uLUT.ColorSpace = 'HSV'
uLUT.UseBelowRangeColor = 0
uLUT.BelowRangeColor = [0.0, 0.0, 0.0]
uLUT.UseAboveRangeColor = 0
uLUT.AboveRangeColor = [1.0, 1.0, 1.0]
uLUT.NanColor = [0.498039215686, 0.498039215686, 0.498039215686]
uLUT.Discretize = 1
uLUT.NumberOfTableValues = 256
uLUT.ScalarRangeInitialized = 1.0
uLUT.HSVWrap = 0
uLUT.VectorComponent = 0
uLUT.VectorMode = 'Magnitude'
uLUT.AllowDuplicateScalars = 1
uLUT.Annotations = []
uLUT.ActiveAnnotatedValues = []
uLUT.IndexedColors = []

# get opacity transfer function/opacity map for 'U'
uPWF = GetOpacityTransferFunction('U')
uPWF.Points = [0.0, 0.0, 0.5, 0.0, 103.3250007861863, 1.0, 0.5, 0.0]
uPWF.AllowDuplicateScalars = 1
uPWF.ScalarRangeInitialized = 1

# Rescale transfer function
uLUT.RescaleTransferFunction(0.0, 110.0)

# Rescale transfer function
uPWF.RescaleTransferFunction(0.0, 110.0)
# reset view to fit data
renderView1.ResetCamera()

# set active source
SetActiveSource(casefoam)

# Properties modified on animationScene1
animationScene1.PlayMode = 'Real Time'

# create a new 'Annotate Time Filter'
annotateTimeFilter1 = AnnotateTimeFilter(Input=casefoam)
annotateTimeFilter1.Format = 'Time: %f'
annotateTimeFilter1.Shift = 0.0
annotateTimeFilter1.Scale = 1.0

# update animation scene based on data timesteps
animationScene1.UpdateAnimationUsingDataTimeSteps( )

# Properties modified on annotateTimeFilter1
annotateTimeFilter1.Format = 'Time: %.1f'

# show data in view
annotateTimeFilter1Display = Show(annotateTimeFilter1, renderView1)
# trace defaults for the display properties.
annotateTimeFilter1Display.Interactivity = 1
annotateTimeFilter1Display.Color = [1.0, 1.0, 1.0]
annotateTimeFilter1Display.Opacity = 1.0
annotateTimeFilter1Display.FontFamily = 'Arial'
annotateTimeFilter1Display.Bold = 0
annotateTimeFilter1Display.Italic = 0
annotateTimeFilter1Display.Shadow = 0
annotateTimeFilter1Display.FontSize = 18
annotateTimeFilter1Display.Justification = 'Left'
annotateTimeFilter1Display.WindowLocation = 'AnyLocation'
annotateTimeFilter1Display.Position = [0.05, 0.05]

# Properties modified on annotateTimeFilter1Display
annotateTimeFilter1Display.FontSize = 12

# Properties modified on annotateTimeFilter1Display
annotateTimeFilter1Display.Shadow = 1

# Properties modified on annotateTimeFilter1Display
annotateTimeFilter1Display.WindowLocation = 'LowerLeftCorner'

# Properties modified on annotateTimeFilter1Display
annotateTimeFilter1Display.WindowLocation = 'UpperLeftCorner'

# set active source
SetActiveSource(casefoam)
animationScene1.GoToFirst()

# current camera placement for renderView1
renderView1.CameraPosition = [-0.2092945091248491, -0.00903714747185367, -0.02075525198175189]
renderView1.CameraFocalPoint = [0.01856183260679245, -0.00903714747185367, -0.02075525198175189]
renderView1.CameraViewUp = [0.0, 0.0, 1.0]
renderView1.CameraParallelScale = 0.12641506499001942

# current camera placement for renderView1
renderView1.CameraPosition = [-0.2092945091248491, -0.00903714747185367, -0.02075525198175189]
renderView1.CameraFocalPoint = [0.01856183260679245, -0.00903714747185367, -0.02075525198175189]
renderView1.CameraViewUp = [0.0, 0.0, 1.0]
renderView1.CameraParallelScale = 0.12641506499001942

# save animation images/movie
WriteAnimation('/usr2/tmp/csnickli/Test_paraview_workflow/pics/pic_komplett.tif', Magnification=2, FrameRate=1.0, Compression=True)

#### saving camera placements for all active views

# current camera placement for renderView1
renderView1.CameraPosition = [-0.2092945091248491, -0.00903714747185367, -0.02075525198175189]
renderView1.CameraFocalPoint = [0.01856183260679245, -0.00903714747185367, -0.02075525198175189]
renderView1.CameraViewUp = [0.0, 0.0, 1.0]
renderView1.CameraParallelScale = 0.12641506499001942

#### uncomment the following to render all views
# RenderAllViews()
# alternatively, if you want to write images, you can use SaveScreenshot(...).

[cfdFreshman]

Does really no one have any idea?
If you need more information or something is unclear just ask.

Thanks
cfdFreshman

[Flowkersma]

Hi,

Everyone probably thinks that you're asking a bit too much. You have provided 1000 lines of code (mostly from tracing tool) and there are many mistakes. It might be easier that you would start from scratch and try to really understand the process. Then if there is something specific that you don't understand, isolate that part and ask from here.

Regards,
Mikko

[cfdFreshman]

Hi,

thank you for your feedback. My intension was to give all available information to the community to get quick answers. I thought my question was as short as possible (q.v. Your knowlege is in demand).
But, somehow the problem ended in smoke a few days ago. The only thing I was doing was set some default values in ParaView for the start sequence of the program. And now it works. Nothing done differently from the previously explained procedure.

Regards,
Stephan