Hi JC, It's not clear to me how surf_int relates to total energy - What is it that you wish to compute? Is it the volume integral of p + 1/2 rho U^2 ? ... or ... ? Regarding computation of WSS, you might also get some ideas by looking at comp_torq() in navier5.f in the source directory. Let us know if this clarifies at all and we'll try to help further. Paul On Mon, 19 Jul 2010, [email protected] wrote:

Hi Nek's,

I am about to run a 3D-stenotic flow DNS, and I would like to compute the total energy in the computational domain. So far, I have only a routine for 2D (surf_int). Does such a 3D routine exists? My second question concern the velocity gradient normal to the wall. Is there any way to output its value at various wall points for each time step? I was thinking ofadding a few lines to the hpts routine, but I am sure about what to write exactly. Could anyone give me a clue?

Sincerely, JC _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users