[kandelabr]

Hello,

I have results from a simulation and I run pimpleFoam -postProcess on them, which creates a postProcessing folder with calculated forces and torques.

How do I load those data in paraview for further calculations (efficiency) and display?

Thanks.

[kandelabr]

Well, if anyone from the future wants to know how I did it, here's the python script for ParaView 'programmable filter'. I needed to know efficiency of a pump so I calculated total torque. You might need something else.

It's not very efficient but for a small number of timesteps it should do its job.

# current time
timestep = inputs[0].GetInformation().Get(vtk.vtkDataObject.DATA_TIME_STEP())
timestep = float(timestep)

# forces: load from forces.dat file
input_file = open('postProcessing/forces_impeller/0/forces.dat', 'r')

for i in range(3):
    # first few lines do not contain data
    input_file.readline()

# read lines until the right timestep has been reached
for line in input_file:
    # values in a single line:
    # 0           1 2 3    4 5 6   7 8 9          10 11 12 13 14 15 16 17 18
    # Time forces(pressure viscous porous) moment(pressure viscous porous)
    values = line.split()
    values = [float(v.strip('()')) for v in values]

    if values[0] == timestep:
        # add pressure and viscous forces;
        # axial
        F_a = values[1] + values[4]

        # maximum radial force
        F_r1 = values[2] + values[5]
        F_r2 = values[3] + values[6]

        F_r = (F_r1**2 + F_r2**2)**0.5


        # torque = pressure + viscous
        T = values[10] + values[13]

        break

# do something with the data