[clarinha]

Hi all,


I'm simulating a river with detailed bed topography in OpenFOAM, using interFoam on a 3D domain. My goal is to obtain the dynamic pressure at the bed.

In ParaView, my rationale is to calculate the vertical distance between the free surface and the bed to obtain the flow depth at each point. Then, I'd calculate the hydrostatic pressure as water density*gravity*flow depth (=rho*g*h). Then, I'd make p minus hydrostatic pressure at the bed to obtain the dynamic pressure at the bed.

As far as I tried using the filters in ParaView, I'm not being successful at obtaining the flow depth! I tried using 'contour' to extract the free surface, setting alpha.water = 0.5. Then I used 'calculator' twice, one to extract coordsZ from the free surface and another to extract coordsZ from the bed. Then I grouped the two calculators using 'group datasets'. Then I used 'calculator' to subtract coordsZ from the two datasets. However, I get the following error message: "Warning: In vtkSMPVRepresentationProxy.cxx, line 288
vtkSMPVRepresentationProxy (000002578E2D7400): Could not determine array range."

Another possible option is to clip the water using alpha.water = 0.5, which would give me coordsZ from the bed to the free surface. But two issues arise: one, the clipped volume does not have 'p' if I use clip by point data, and two, I don't know how to proceed to subtract the elevation between the free surface and the bed.


Any help is appreciated!
Cheers