Query about the output pressure about incompact3d

[mathrack]

Hello,

I think you are correct. With the implicit time scheme for the y viscous term, one gets :
$p = \phi + \frac{\Delta t}{2 Re} \partial_{yy} \phi$ when using Crank-Nicolson
and $p = \phi + \frac{\Delta t}{Re} \partial_{yy} \phi$ when using implicit Euler

[mattfalcone1997]

Hi,

Thanks for your response. I was wondering whether this would subsequently require modification of the visualisation so that the pressure in the binary file corresponds to the physical pressure. This is also relevant for the statistics.

[mathrack]

This definition of the pressure is selected so that when summing prediction and correction and inserting $p$ instead of $\phi$, the predicted velocity does not appear in the viscous term.

Regarding the exact meaning of physical pressure in the framework of incompressible flows, ...

[mattfalcone1997]

I am only concerned with whether the pressure output in postprocessing should be modified to account for the implicit scheme. I do not believe that this is done at the moment (it seems that $\phi$ is output rather than $p$, although I may be wrong).

[mathrack]

The current code outputs $\phi$ and not $p$ when using the implicit scheme.

[mattfalcone1997]

Okay, thanks for the clarification. Is this something that should/will be updated at a future point? I am not sure how important this is as the difference is $O(\Delta t)$.

[mathrack]

I am not planning to open a PR for this modification. The estimation of $\partial_{yy} \phi$ could be inserted just before the call to transpose_y_to_x : https://github.com/xcompact3d/Incompact3d/blob/master/src/visu.f90#L245

[mattfalcone1997]

Thanks for your help. I may have a play around with this. I am happy to close this now.