To extract point values from 2D SPHysics simulations

[vivekfrancis139]

Dear all,

I'm a PhD Research Scholar working on the 2D version of SPHysics codes and I've been making modifications to it. And these modifications, I can visualize through animations obtained via ParaView software.

The trouble I'm having is, though the animations are appearing perfectly, I'm not able to extract point values of water surface elevations ParaView (for instance from Case 3: Waves generated by a paddle in a beach). I placed wave gauges in desired locations after having loaded the '.vtu' files into ParaView, but when I plot the evolution of values over time, I'm not able to get any data at all. Each point value is always a zero. This could be because of a very basic error and I'm trying to sort it out (a similar approach in ParaView when applied to REEF3D software was giving me results!).

I really look forward to your answers. Any suggestion would be really appreciated.

[sarfaraz_mohammad]

Hi
You can NOT compute the water elevation directly from the paraview.
Imagine you want to compute water elevation at x=x0.
Just write a simple code of Matlab or fortran, to gather all particle located at x0-a<x0<x0+a. These particle should have density grater that 0.5*water_density.
"a" can be about 2cm, depending on your case. You need to validate it.
Then the particle which has the maximum elevation denotes your water elevation.

Regards!

[vivekfrancis139]

Thank you for your suggestion.I'm very grateful.
I've implemented the method and it is working well.